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Abstract. We study the Cauchy problem for the sine-Gordon equation in the semiclassical Umit with 
pure-impulse initial data of sufficient strength to generate both high-frequency rotational motion near the 
peak of the impulse profile and also high-frequency librational motion in the tails. We show that for small 
times independent of the semiclassical scaling parameter, both types of motion are accurately described by 
explicit formulae involving elliptic functions. These formulae demonstrate consistency with predictions of 
Whitham's formal modulation theory in both the hyperbolic (modulationally stable) and elliptic (modula- 
tionally unstable) cases. 
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1. Introduction 



The sine-Gordon equation in laboratory coordinates 

e^utt - e'^Uxx + sin{u) ^ 0, xeR, t > 0, (1.1) 

is a model for magnetic flux propagation in long superconducting Josephson junctions ^2], but perhaps may 
be most easily thought of as the equation describing mechanical vibrations of a "ribbon" pendulum (the 
continuum limit of a system of linearly arranged co-axial pendula with nearest-neighbor torsion coupling). 
These and other applications are discussed in detail in the review article [5- The correct Cauchy problem 
for this equation involves determining the solution consistent with the given initial data 

u{x,0) = F{x) and €Ut{x,0) ^ G{x). (1.2) 

Here F{x) and G{x) are independent of the fixed parameter e. This Cauchy problem is globally well-posed 
[4]: if p > 1 and F, F' , and G are functions in LP(IR) then there is a unique solution with u, u^, and Ut all 
lying in LJ^^(IR_|_; L^'(]R)). Moreover if the initial data have one more derivative, so that F" and G' are 
functions in LP{M.), then this further regularity is preserved as well so that u^x and Utx lie in LJ^^(K_|_; LP{M.)). 
These well-posedness results also hold in a slightly modified form when the initial displacement has nonzero 
asymptotic values: F(x) — >■ 27m±, n± g Z, as a; — >■ ±oo. In this case the topological charge n+ — rt_ is 
preserved for all time in the solution u. 

If suitable initial conditions F and G are fixed, one may therefore in principle construct the unique global 
solution u{x, t) of ( 1.1 ) subject to ( 1.2 1 for each positive e. Our interest is in the asymptotic behavior of this 
family of global solutions in the semiclassical limit e — >■ 0. The well-known elementary excitations of the sine- 
Gordon equation include solitons of kink (or antikink) and breather type; these have a width proportional to 
e while the length scales in the initial conditions (1.2) are independent of e. This suggests that when e ^ 1 
the initial conditions of the system can be viewed as preparing a "condensate" whose ultimate breakup will 
liberate approximately 1/e fundamental particles. 

The decay process will take some time to become complete, and during the intermediate stages one may 
expect that some solitons may partially emerge from the condensate moving with nearly identical velocities, 
thus forming a modulated wavetrain. The simplest models for these wavetrains are the periodic (modulo 
2tt) traveling wave exact solutions of (1.1 ) of the form 

'<pix,ty 



u{x, t) = U 



U{C + 2n) = U{0 (mod27r), <l>{x,t) ^ kx - ujt, 



(1.3) 



where k is the wavenumher and lu the frequency of the wavetrain. With this substitution, the sine-Gordon 
equation reduces to an ordinary differential equation that can be integrated once to 

2 



2^ 



dU 



cos([/) = E 



(1.4) 



where E is an integration constant having the interpretation of energy. There are four types of solutions 
subject to the periodicity condition: 

• Superluminal librational wavetrains correspond to > k^ and \£\ < 1. From a phase portrait it is 
evident that J7(C -|- 2tt) = U{(^) if the nonlinear dispersion relation 

-I -2 



-k^ = 27r^ 



H-£) \/cos{(l)) 



(1.5) 



is satisfied, and then U oscillates about a mean value of U 
less than tt. 



(mod 27r) with an amplitude strictly 



Suhluminal librational wavetrains correspond to < fc^ and \E\ < 1. From a phase portrait it is 
evident that U{C + 2tt) 



[/(C) if the nonlinear dispersion relation 

-I -2 

/•+ cos (^i- ) 

■2 - fc2 = -27r2 



s-i(£) Y^cos((/>) - £ 



(1.6) 



is satisfied, and then U oscillates about a mean value of U 
less than tt. 

2 



TT (mod 27r) with an amplitude strictly 



• Superluminal rotational wavetrains correspond to cj^ > fc^ and £ > 1. Here the phase portrait 



indicates that C/(C + 27r) 



f/(C) ± 27r if the nonhncar dispersion relation 

-I -2 



V/COS(0)~T^ 



(1.7) 



is satisfied, with [/'(C) strictly nonzero being largest in magnitude when C = (mod 27r). 
• Suhluminal rotational wavetrains correspond to < and £ < —1. Here the phase portrait 
indicates that U{C + 2tt) = U(() ± 27r if the nonlinear dispersion relation 

-2 



-TT ^cos((/)) - £ 



(1.8) 



is satisfied, with [/'(C) strictly nonzero being largest in magnitude when C = (mod 27r). 

In the classical mechanics literature |13j the term libration is used to characterize the kind of motion in 
which both position and momentum are periodic functions while the term rotation is used to characterize 
motions in which momentum is periodic but position is not because the momentum has a nonzero average 
value. Furthermore, the dichotomy of subluminal waves versus superluminal waves is important because the 
sine-Gordon equation is strictly hyperbolic with characteristic velocities Vp = ±1; thus subluminal waves 
have phase velocities bounded in magnitude by the characteristic velocity, while superluminal waves move 
faster than the (unit) characteristic speed. In the superluminal (respectively, subluminal) case, the energy 
value oi £ = 1 (respectively, £ = —1) corresponds to the separatrix in the phase portrait of the simple 
pendulum, at which point the period (respectively, wavelength) of the waves tends to infinity. In this limit, 
each of the four types of wavetrain degenerates into a train of well-separated kink-type solitons; for rotational 
waves the kinks all have the same topological charge, while for librational waves the pulses alternate from 
kink to antikink for zero net charge. 

In a body of work beginning with his seminal 1965 paper |24', Whitham developed a nonlinear theory of 
modulated wavetrains. The main idea in the current context is that one seeks solutions of the sine-Gordon 



equation ( 1 . 1 1 of the approximate form 



u{x, t) = U 



$(a;,i) 



0{e) 



(1.9) 



over space and time intervals of 0(1) length, where two essential changes are made in the leading term: 

(1) The parameters k, uj, and £ are no longer taken as constant, but are allowed to depend on {x,t) as 
long as the appropriate nonlinear dispersion relation is satisfied pointwise, and 

(2) The phase ^{x,t) is replaced with a general (nonlinear) function of {x,t) and the local wavenumber 
and frequency are derived therefrom by the relations 

- and ^.,t):=-_. 



k{x, t) 



(1.10) 



By consistency, the definition (1.10 1 imposes that k and uj are necessarily linked by 

dk duj 
dt dx 

an equation that expresses conservation of waves. It then follows that for the error term in (1.9) to remain 



= 0, 



(1.11) 



0{e) formally, one additional partial differential equation on k, uj, and £ , only two of which are independent 
due to the nonlinear dispersion relation, is required to hold. This equation may be derived by many different 
methods. Perhaps the most direct in this context is to appeal to an averaged variational principle [25j . The 



sine-Gordon equation (1.1) is the Euler-Lagrange equation for the variational principle: 

S 

du , 



J J L[u\ dxdt = 
where the Lagrangian density is 



+ sin(u) — 



1 



2„,2 

't 



2 2 



cos(u) 



(1.12) 



(1.13) 



having the interpretation of the difference between kinetic and potential energy densities. The procedure is 



to substitute the exact wavetrain into L, using the differential equation (1.4) satisfied by U to simplify the 
resulting expression: 



U 



kx — Lot 



\{u? - e)U\C? + cos((7(C)) - 2(5 + cos([/(C))) - 



(1.14) 



where for the exact solution — kx — ujt. This expression is periodic in with period 27r, so one may define 
its period average as 

(L) := - + cos(C/(C))) dC - £. (1.15) 



An exact expression for U{(^) is not necessary; using the differential equation (1.4), a shifted version of U 
may be used as the integration variable although the details are slightly different in the four cases; defining 
integrals 

/l(£) ^ / v/cos(^) - f (i</) > 0, -1< £ < I, 

^ J-cos-i(£) 

and 



1 /"^ 

Ir{£):=—^ \/cos(^) - f > 0, f < -1, 

and noting for future reference that 

/l (f ) > for -1< f < 1 while /r(£:) < for f < -1, 

the result is that 

{L) = J{£)^l.y\uJ^ - k^\ - £ 
where fi = sgn(ci;^ — fc^) distinguishes the superluminal and subluminal cases, and where 

Ji£) := hi-^i£) or J(f ) := /R(-Aif ) 



(1.16) 
(1.17) 



(1.18) 
(1.19) 
(1.20) 



depending on whether we are considering librational or rotational wavetrains, respectively. One then substi- 
tutes k = 9x and lo = —9t and formulates the averaged variational principle: 

jj (L) dxdt = and J J (L) dx dt ^ 0. (1.21) 

The first of these two equations reproduces in each case the corresponding nonlinear dispersion relation 

d{L) 



( 1.5 )-( 1.8 ) in the form 



d£ 



f{£)^i^/\uJ^~k^\-l = 0. 



The second is a first-order partial differential equation: 



d 


\ dm 


a 


\d{L)] 


dt 


du! 


dx 


dk 



= 0. 



(1.22) 
(1.23) 



This equation together with ( |1.11[ ) and the nonlinear dispersion relation (1.22) to eliminate one of the 
three variables yields a closed system of equations to determine these fields as functions of {x,t). From the 
exposition in |22j one learns to appreciate the utility of taking £ and the phase velocity 

Vp (1-24) 



as the two unknowns, and thus the calculations go as follows. Clearly, one has 



d{L) 
du! 



ujj{£) 



fc2| 



and 



d{L) 
dk 



kJ{£) 



fc2| 



(1.25) 



be eliminated in favor of Up and £ 



Using uj = fpfc together with the nonlinear dispersion relation in the form (1.22) shows that k and oj may 



k 



J'{£)^\vl-l\ 



and UJ — Vpk 



J'i£)JK~l\ 



(1.26) 



where cr = ±1 is an arbitrary sign whose role is to select different branches of the dispersion relation. 
Therefore, the variational modulation equation ( |1.23 1 becomes 



d 


VpJ{£) 


d 


J{E) 


di 




dx 


[jk-i\\ 



= 



(1.27) 



and the conservation of waves equation (1.11 ) becomes 



d_ 
di 



d_ 
dx 



(1.28) 



(the sign cr = ±1 drops out in each case). An application of the chain rule puts the system in the form 



d_ 
di 



1 



-J{e)J'{£) 



{I -vlf J' {£).]"{£) 
v^[J {£).]"{£) + J' {£f 



' d 


'Vp 


dx 


£ 



= 



(1.29) 



where 



V{vp,£) ■.= J{£)J"{£) + vlJ'{£f 



(1.30) 

Thus the dependence on the sign /x — ±1 also disappears except from within the definition (1.20) of J{£)- 
Actually, this system as written has an apparent singularity if Wp blows up as can happen in the superluminal 
cases when the wavenumber k vanishes. In these cases it is better to introduce the reciprocal phase velocity 

Up := - (1.31) 



and then in the overlap region where neither Vp nor rip vanishes, (1.29) takes the form 

np[J{£)J"{£) + J'{£)^ 
J{£)J'{£) 



d_ 
dt 



£ 



J\f{np,£) 



-{l-nlfj'{£)J"{£) 
np[J{£)J"{£) + J'{£f] 



' d 


rip 


dx 


£ 



where 



AA(np,f) := nlJ{£)J"{£) + J' {£f 



This latter form of the Whitham modulation equations has no apparent singularity when 
sponding to Up cx). 



(1.32) 

(1.33) 
corre- 



The characteristic velocities c = Cj, j = 0, 1, are the eigenvalues of the coefficient matrix of x-derivatives 



and therefore are the roots of the quadratic equation 



{vp[J[£)J"{£) + J'{£f] - V{vp,£)cf + (1 - vlfj{£)J'{£fj"{£) = 



(1.34) 



or, equivalently as the coefficient matrices in (1.29) and (1.32) are similar, 

{np[J{£)J"{£) + J'{£f] -M{np,£)cf + (1 - n^f J {£) J' {£Y J" {£) = 0. (1.35) 

Since V{vp,£) and Af{np,£) are real the Whitham modulation system in either form (1.29) or ( 1.32[ ) is 
hyperbolic (corresponding to real and distinct characteristic velocities) if and only if J{£)J"{£) < 0. Using 
(1.20) and ( 1.18[ ) then shows that the Whitham systems governing modulations of rotational waves (both 
types, superluminal and subluminal) are hyperbolic, while those governing librational waves (again, both 
types) are elliptic. The Whitham modulation theory has been generalized to handle modulated multiphase 
waves [121 Uni E] having any number of 27r-periodic phase variables; however the full implications of the 
resulting modulation equations generalizing (1.32) have apparently only been understood for waves with a 
maximum of two phases. 

The Whitham modulation theory makes predictions of so-called "modulational stability" of wavetrains on 
the basis of whether the quasilinear system of modulation equations is hyperbolic (modulationally stable) or 
elliptic (modulationally unstable) . One should think of modulational stability as linear (neutral) stability of 
perturbations to the wavetrain that have similar characteristic wavelengths and periods to the unperturbed 
exact wavetrain solution. Thus, hyperbolicity of the modulation equations suggests the absence of a "slow" 
sideband instability, but does not necessarily rule out instabilities to perturbations with wavenumbers far 
from the unperturbed wavenumber k or even sideband instabilities with exponential growth rates that are 

5 



far from the unperturbed frequency uj. Among the candidate wavetrain types for stabiHty of the hnearized 
equation 



u 







(in the sense of a L^(R) estimate on v and evt that depends on initial data but is independent of t) there 
are thus only the subluminal and superluminal rotational wavetrains. It is not difhcult to believe that the 
superluminal rotational wavetrains are linearly unstable, since in the limiting case of fc = and £ | 1 one 
obtains an orbit homoclinic to the exact constant solution u{x,t) = (2m + l)7r, to e Z which is obviously 
unstable to small spatially constant perturbations. Indeed, such perturbations cause all of the vertical 
pendula to "drop" simultaneously; the growth rate of the perturbation is clearly large compared to the zero 
frequency of the unperturbed solution explaining why the instability is not captured by Whitham theory. 
The definitive statement in the literature [21\ is that the subluminal rotational wavetrains are indeed linearly 
stable (and the only stable type among the four types of traveling waves), although we have not been able 
to verify the line of argument in all details (it is not clear to us from the proof that the stable solutions 
obtained in the case of subluminal rotational wavetrains form a basis of L^(M), or that the exponentially 
growing solutions obtained in the other three cases are relevant as they appear to be unbounded in x). 

In this paper, we will analyze the Cauchy problem for (1.1 1 with initial data F and G independent of e, 
in the semiclassical limit e — ^ 0. We will show that for a general class of "pure impulse" initial data, most of 
the real x-axis is occupied for small time (independent of e) by modulated superluminal wavetrains of either 
rotational or librational types. This result shows the relevance of the Whitham modulation theory even in 
some cases when it results in elliptic modulation equations. (In a forthcoming paper ^ we will show that 
the remaining part of the x-axis is occupied by more complicated oscillations that nonetheless have a certain 
universal form for t small.) 

1.1. Pure impulse initial data for the sine-Gordon equation. Connection to the Zakharov- 



Shabat scattering problem. Equation (1.1) is the compatibility condition for the Lax pair 

AE{w) (1 — cos(u)) sin(u) — ie{ux + Uf) 



: sin(u) + ie{ux + Ut) — 4_B(w) + 



4Diw) + (1 - cos(u)) 



: sm{u) + ie{ux + Ut) — 4Z?(w) — 



: sin(u) 



;(1 — cos(u)) 

- ie{u^ + Ut) 
-{1 — cos(m)) 



where 



Eiw) := - 



and 



D{w) := - 



(1.37) 



(1.38) 



(1.39) 



and w G C \ M-i- is the spectral parameter. Here and throughout this paper, the radical refers to the principal 
branch of the square root. Analysis of the Cauchy problem for (1.1) posed with initial data (1.2) may be 



A self-contained account of this analysis can 
i\J —w is used; the utility of w is related to an 



carried out in some detail by means of the inverse scattering transform based on the differential equation 
(1.37), the so-called Faddeev-Takhtajan eigenvalue problem, 
be found in our paper 1^, where the spectral parameter z = 
even symmetry of the spectrum in the z-plane. 

By pure impulse initial data we simply mean data for which the initial displacement F(x) vanishes 
identically. An elementary observation is that if F{x) = 0, then for < = the Faddeev-Takhtajan eigenvalue 
problem reduces to the Zakharov-Shabat eigenvalue problem: 



-i\ 

'l\}[x) 



il}(x) 
i\ 



1 



G{x), X:^E{w). 



(1.40) 



This is the eigenvalue problem in the Lax pair for the sine-Gordon equation in characteristic coordinates 
and the cubic focusing nonlinear Schrodinger equations. Although when viewed as an eigenvalue problem of 
the form £v — Av the operator C is nonselfadjoint, several useful facts are known about the spectrum of C 
in the case relevant here that ipi^) is real. It is easy to see that for real ip the discrete spectrum comes in 

6 



quartets symmetric with respect to the involutions A i— ?> A* and A i— > —A. Moreover, Klaus and Shaw [TS] 
have proved that ii tp £ L^{M.) D C^(M) is real, of one sign, and has a single critical point (so that the graph 
of is "bell-shaped"), then the discrete spectrum is purely imaginary and nondegenerate. 

In the context of pure impulse initial data for which G is a nonpositive (without loss of generality) function 
of Klaus-Shaw type, the necessarily purely imaginary eigenvalues A may be approximated by a formal WKB 
method applicable when e ^ 1. The result of this analysis is that with the Klaus-Shaw function G{x) one 
associates the WKB phase integral 

^'(A) - / ^ \/G{s)^ + 16A2 ds, < y := -iA < max f--G(x)V (1.41) 

where a;_ (A) < x+ (A) are the two roots of G{s)'^ + 16A^ when A is as indicated. Then one defines approximate 
eigenvalues A^ by the Bohr-Sommcrfeld quantization rule 

^{XD^TTefk+l) , k^0,l,2,...,N{e)-l, (1.42) 



(1.43) 



2^ 

where the asymptotic number of eigenvalues on the positive imaginary axis is 

1 1 



N{e) 



2 ' 4.."^"^ 



where || • ||i denotes the standard i^(M) norm. To each simple eigenvalue A of the Zakharov-Shabat eigen- 
value problem (1.40) there corresponds a proportionality constant 7 relating the solution having normalized 
decaying asymptotics as x — >■ —00 with the solution having normalized decaying asymptotics as a; — >■ -l-oo. If 
G is a real even function of x, then one can show by symmetry that 7 = ±1, and WKB theory for Klaus-Shaw 
potentials suggests that the proportionality constants alternate in sign along the imaginary axis. Thus for 
even G, to the approximate eigenvalue A^ we associate the approximate proportionality constant 7^ given 

by 

:= (-l)'=+i, fc = 0,l,2,...,iV(e)-l, for even G. (1.44) 
The final ingredient of the scattering data for G in the Zakharov-Shabat problem (1.40) is the reflection 
coefficient defined for real A. According to the WKB approximation, the reflection coefficient is small 
pointwise for A 7^ when e <C 1. 

In this paper, we study even, pure impulse initial data of Klaus-Shaw type for the sine-Gordon equation 
(1.1). Thus we assume 

Assumption 1.1. The initial conditions (1.2) for ( |1.1[ ) satisfy F{x) = 0. 



Assumption 1.2. In the initial condition ( |1.2[ ) for eut, the function G{x) is a nonpositive function of 
Klaus-Shaw type, that is, G G iy"'^(M) n G"'^(M) and G has a unique local (and global) minimum. 

Assumption 1.3. The function G is even in x: G{—x) = G{x), placing the unique minimum of G at x = 0. 



Note that under Assumptions 1.2 and 1.3 
(G(0), 0) — >■ M+, and in terms of it we may rewrite the WKB phase integral in the form 



G restricted to M+ has a unique inverse function G ^ 



*(A) = 



1 



VG(s)2 - w2 ds, A = 



0<v< -G(0). 



(1.45) 



Thus, ^{iv/A) defined on (0, — G(0)) is an Abel-type integral transform of the nondecreasing function G{x) < 
defined on IR+. 



Proposition 1.1. On its range, the inverse of the transform (1.45) is given by the formula 

-'^^(O) ip{v)dv 



G-^iw) = - 



V — w 



G(0)<w<0, 



where 



(1.46) 
(1.47) 



7 



The proof of this proposition is a rather straightforward appUcation of Fubini's Theorem and is given in 
Appendix |A] 

We wiU require that the WKB phase integral have certain analyticity properties to be outlined in Propo- 
sition |1.2| below. We now make an assumption on G that will be sufficient to establish Proposition |1.2| and 
that can easily be checked for a given G: 

Assumption 1.4. The function G is strictly increasing and real-analytic at each x > 0, and the positive 
and real-analytic function 

WGiOr - n 



^(m) 



< m < GiOy 



(1.4 



2G"(G-i(-Vr7i))' 

can be analytically continued to neighborhoods of m ~ and m — G{0)^, with ^^(0) > and ^^(G(0)^) > 0. 

We point out that the class of functions '^{m) satisfying Assumption |1 .4| obviously parametrizes a corre- 
sponding class of admissible functions G{x) by simply viewing (1.48 1 as an equation to be solved for x = G~^ 
^ . The solution is: 

■ (1.49) 



given 



'^{m) dm 



m 



VgW 



For example, the function ^{m) = C > clearly satisfies the analyticity and positivity conditions on ^ 
listed in Assumption 1.4 and in this case the integral in (1.49) can be evaluated in terms of elementary 
functions and the resulting function x = G~^{G) can be inverted to yield 



G = G(0) sech(^^x 
perhaps the simplest example of an admissible initial condition. 



(1.50) 



Proposition 1.2. Suppose that Assumptions 1.2-1.4 hold. Then the function ^'(A) defined by (1.45) for 

X ~ iv/A and < w < — G(0) is positive and strictly decreasing (to zero) in v. Furthermore, is real-analytic 
for < w < — G(0) and has an analytic continuation to neighborhoods of v = and v = — G(0), for which 

d 



and, for some S > 0, 



*(A) = 



*(A) 



and — *(A) < 0, for A 
dv 



-zG(0)/4, 



|G|| 



iA + ^/3„A2", \\\<6 



(1.51) 



(1.52) 



n=l 



where a > and f3n G M for all n. 



We provide the proof of this statement in Appendix |Aj In particular, Proposition |1.2| guarantees the 
existence of a simply-connected open set S C C containing the closed imaginary interval < —Ai\ < — G(0) 
in which 4* (A) may be considered as a holomorphic function of A whose restriction to that interval is a 
real- valued function given by (1.41 1. By Schwartz reflection we therefore will have 



5r(-A*) = *(A)*, AeS, -A*eS, 



(1.53) 



which also shows that without loss of generality we may simply take S to be symmetric with respect to 
reflection through the imaginary axis. By the strict monotonicity and reality of ^'(A) for < — 4iA < — G(0) 
it follows from the Cauchy-Riemann equations that there exists some positive number Si < S such that for 
A in the open rectangle D+ {A e S, < 5R{A} < Si, < 5{A} < -G(0)/4} the inequality S{^ (A)} > 
holds, with the strict inequality failing only as A approaches the imaginary axis. According to ( |1.53[ ), if 
A e £)_ := -D;, then 9{^'(A)} < 0. 

Our analysis will require an assumption about e > 0: 



Assumption 1.5. The small number e lies in the infinite sequence 

''^"^ iV = l,2,3,... 



47rA^ 



(1.54) 



For such e we have from (1.43) that N{e) — N. 



Clearly, according to (1.52), Assumption 1.5 implies that 

*(0) 



ttA^, iV==l,2,3, 



(1.55) 



Also, according to WKB theory, another implication of Assumption |1.5| is that the reflection coefficient is 
uniformly small for A G M, i.e. the choice (1.541 makes the reflection coefficient negligible in a neighborhood 
of A = 0. 

In fact, our strategy will be to replace the scattering data corresponding to initial data of the above 
type with its WKB approximation, admittedly an ad hoc step, and then to carry out rigorous analysis 
of the inverse-scattering problem in the limit e — > 0. The sequence of exact solutions of (1.1) generated 
by the spectral approximation procedure and indexed by TV = 1, 2, 3, . . . is an example of a scmiclassical 
soliton ensemble in the sense of (TS]. In the context of the sine-Gordon equation and its application to 
superconducting Josephson junctions we will call this sequence {uN{x,t)} of exact solutions the fluxon 
condensate associated with the impulse profile G{x). The accuracy of this procedure for studying the Cauchy 
problem is suggested by the fact that the fluxon condensate recovers the initial data at i = to within an 
error of 0{e) (see Corollary |1.1| below). Also, there exist special cases for which the fluxon condensate 
represents the exact solution of the Cauchy problem when e lies in the sequence e = e^v (see (1.54)), giving 
further justification to the procedure. 

1.2. Exact solutions. Impulse threshold for generation of rotational waves. If furthermore we 
suppose that the sine-Gordon system ( |1.1| on M is set into motion at i = by pure impulse initial data of 
the special form 

F{x) = and G{x) ^ -4Asech{x) (1.56) 
for some ^4 > 0, the Zakharov-Shabat problem reduces further to a special case that was studied by Satsuma 
and Yajima [30]. From their work it follows that ii e — eN ■— A/N for N = 1,2,3,..., the reflection 



coefficient for (1.401 vanishes identically as a function of A, and the eigenvalues in the upper half-plane are 
the purely imaginary numbers X — iA — i{k + \)^n^ for fc = 0, 1, . . . , A'^ — 1. The auxiliary scattering data 
consist of the proportionality constants linking eigenfunctions with prescribed decay as a; — > — oo with those 
having prescribed decay as a; — > ~\-co: these are simply alternating signs (—1)'''+^. It is easily confirmed that 
these scattering data correspond exactly to the WKB approximation described above when G is defined as 



above; in particular, the phase integral (1.41) evaluates to 

«'(A) = iirX + ttA, if G{x) = -AA sech(a;). (1.57) 
In such a reflectionless situation, it becomes possible to solve the inverse-scattering problem by finite- 



dimensional linear algebra, and thus we obtain an exact solution u{x,t) of (1.1) with e — ejq and e 



independent initial data (1.56) for each positive integer N. We will give more details about this procedure 
later (see the final paragraph of ^ , but for now we discuss the results of an empirical study of these exact 
solutions. 

Figure [O] shows plots of the exact solution of the Cauchy problem for (1.1) with initial conditions given 
by ( 1.56[ ) with A~l/A. These plots show that the initial impulse sets the pendula into nearly synchronous 
librational motion of frequency proportional to N (inversely proportional to e). However, an instability 
seems to appear in the modulational pattern, leading to a kind of focusing of wave energy near a; = 0. 
In a region of the (a;,t)-plane that seems to become more well-defined as N increases, the focused waves 
take on a different character; in particular, the synchrony of nearby pendula is lost as spatial structures 
with wavelengths inversely proportional to N spontaneously appear. Nonetheless, the oscillations present 
after the focusing event are still fairly regular and result in an approximately quasiperiodic spatiotemporal 
pattern. This type of dynamics is qualitatively similar to what is known to occur in the scmiclassical limit 
of the focusing nonlinear Schrodinger equation [T^ [TH [531 HZ] , another problem that has elliptic modulation 
equations as is expected here before the focusing due to the apparent librational motion. 

In Figure 1.2 we give plots analogous to those in Figure |1.1| but now we increase the amplitude of the 
impulse by choosing A = 3/4. The dynamics are evidently quite different: one clearly can identify three 
types of behavior near t — Q: 

• Nearly synchronous librational motion of the pendula is apparent for large where the initial 
impulse is weak. 
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-1 -0.5 0.5 1 

cos(u) 

Figure 1.1. Plots o/ cos(u) with A — 1/4 over the {x,t)-plane. The horizontal axis is 
—2.5 < a; < 2.5 and the vertical axis is < t < 5. Left: N = 4. Center: N = S. Right: 
N = 16. 
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-1 -0.5 0.5 1 



cos(w) 

Figure 1.2. Same as in Figure [T7| but now A = 3/4. 



• Nearly synchronous rotational motion of the pendula is apparent for small \x\, where the initial 
impulse is strongest. 

• Strongly asynchronous motion of undetermined type is apparent near two transitional values of x. 
These transitional points appear to shed kinks and antikinks. 

These plots therefore suggest that some sort of transition in the dynamics occurs when the maximum 
amplitude of the impulse, A, exceeds some threshold value between A = 1/A and A — 3/4. In our paper 
[3j we observed that in similar families of exact solutions an analogous transition occurs if the initial data 
{F, G) when viewed as a curve (parametrized by x) in the phase portrait of the simple pendulum equation 
e^Utt +sin(w) = crosses the separatrix (and the transition appears to occur near the specific x-values where 
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crossings occur). Thus one expects that the existence of a; G M where 

(l-cos(F(a;))) + iG(x)2 = 2 



(1.58) 



may lead to a different kind of dynamics in a small time interval near < = of length independent of e. 
This equation has real solutions when F{x) and G{x) are given by (1.561 ii A > 1/2. From the plots in 
Figure 1.2 it is clear that when A > 1/2 there is a symmetrical interval around a; = in which the impulse 
is sufficiently large to cause rotation of the angle u outside of the fundamental range — tt < u < tt, which 
leads to an emission of kinks carrying positive and negative topological charge in opposite directions from 
the ends of this (shrinking in time) interval, which bound a triangular region in the (x, t)-plane. This region 
therefore appears to contain a modulated superluminal rotational wavetrain, while for \x\ sufficiently large 
one observes what appear to be modulated superluminal librational waves. 

The goal of this paper is to show that this type of structure is universal for pure-impulse fluxon condensates 
with sufficient impulse present at t — 0. Thus we further impose 

Assumption 1.6. The function G{x) satisfies G{Q) < —2. 

By Assumptions |1.2[ |1.3[ and |1.6[ there exists a positive number Xait > defined by 

a^crit :=G-i(-2), (1.59) 

and based upon the above heuristic discussion we may expect the dynamics of pure-impulse fluxon conden- 
sates to be of a different character for \x\ < XcHt than for |x| > Xcrit- 

1.3. Statement of results. Our results concern the asymptotic behavior, in the limit iV t oo equivalent to 
Eat ^ 0, of the functions uiq{x,t) making up the fluxon condensate associated with the pure-impulse initial 
condition of impulse profile G(-). As mentioned above, for fixed N , un{x, t) is not exactly the solution of the 
Cauchy initial- value problem with the corresponding initial data (although it is an exact solution of (1.1)), 
and the proper definition will be given below in Definition |2.1[ The statements below concern two regions 
of the (2:,t)-plane depending on G(-) but not on e = eat. The region S\_ is specified in terms of a continuous 
time-horizon function T\_{x) > defined for \x\ > Xcrit with \iuY\^\t^^^^^^T\_[x) = 0; then {x,t) G S\_ if and 
only if \x\ > Xait and \t\ < Ti{x). The region S'r is specified in terms of a continuous time-horizon function 
Tf((x) > defined for < cccrit with lim|2,|^2.^^.j. Tfi{x) = 0, as well as two curves t = t±{x) with t±{0) = 0, 
t'+ix) > 0, and t'_{x) < 0. Then (x,t) G S'r if and only if \x\ < Xcrit, \t\ < Tr{x), and t ^ t±{x). See 
Figure [Ol 

t 

n 



Si 



■^crit 



t+{x) 



t+{x) 
Sr 



Sl 



-►X 



■^crit 



t.{x) 



Figure 1.3. The regions Sl and Sr. 



Let K{-) denote the complete elliptic integral of the first kind: 

ds 

lo v/(l-s')(l-^' 



K{m) 



< m < 1. 



(1.60) 



We will prove the following two results under Assumptions 1.1 - 1.6[ and also under the more technical 



Assumption 2.1 to be presented shortly. The asymptotic formulae for the fluxon condensate are in fact quite 
explicit in terms of the Jacobi elliptic functions sn(z;m), cn(z;TO), and dn(z;m), for which we cite the text 
by Akhiezer [T] as a reference. 

Theorem 1.1 (Small-Time Librational Asymptotics). There exist well-defined differentiable functions Hp : 
Si (—1, 1) and £ : Si (—1, 1) satisfying the initial conditions np{x, 0) = and £{x, 0) — ■^G{x)'^ — 1 
and the elliptic Whitham system (1.32) where J{S) := Ii{~£) and Ii{-) is given by (1.16), as well as the 
inequality 

dur. 



dt 

Defining an elliptic parameter m — m{x,t) by 



x^:r^{x,0) < 0, > Xcrit- 



m(a:, t) = mi{x, t) := ^ ^ ^J^' G (0, 1), {x, t) e Si, 



and a real phase t) by 



where 



<^{x,t) 



t') dt' , 



ijj[x, t) 



2K{m{x,t)) ^l-np(a;,t)2' 
the following asymptotic formulae hold pointwise for {x, t) G Si: 

, f 2^(x,t)K{m{x,t)) , \ , 



(1.61) 

(1.62) 

(1.63) 
(1.64) 



cos [ -UN{x,t) 



sm 



1 , ,A r^—, ( 2^{x,t)K{m{x,t)) , 
-UN{x,t) j = -y'm{x,t)sn I ■,m{x,t) 1 + 0[eN) 



(1.65) 



duN 

^N^{x,t) 



4:K{m{x,t)) a$ 
n Ik 



\/m{x,t) 



2^(x,t)K(m(x,t)) , , , , 



Moreover, the error terms are uniform for {x,t) in compact subsets of Si. The phase $(x,t) also satisfies 
d^/dx — k{x, t) := uj{x, t)np{x, t). 



The accuracy of the asymptotic formulae (1.65) is illustrated in Figure 1.4 



COs(ttAr) 



COs(wAr) 





Figure 1.4. The cosine of the exact solution UN{x,t) for the special initial data described 



in i'1.2 for A = 3/4 plotted with points, and the corresponding asymptotic formula plotted 
with curves, for fixed x = 1.5 and < t < 0.5 so that {x,t) € Si yielding librational motion 
as described by Theorem \l.l\ Left: N — 8, or equivalently — 0.09375. Right: N — 16, 
or equivalently = 0.046875. 
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Theorem 1.2 (Small-Time Rotational Asymptotics). There exist well-defined differentiable functions Up : 
Sr (—1, 1) and £ : Sr —>■ (1, +oo) satisfying the initial conditions np{x, 0) = and £{x, 0) = ■^G{x)^ — 1 
and the hyperbolic Whitham system (1.32) where J{£) = Ir{—£) and /r(-) is given by (1.17), as well as the 
inequality 



x^{x,0)>0, 



< |a;| < Xcrit- 



(1.66) 



Moreover, the functions np(x,t) and £{x,t) extend continuously along with their first partial derivatives to 
the curves t = t±{x), and they obey the {x^f)- independent inequalities 



and 



< 
1 

r 



1 - np{x,t) 
1 + np{x, t) 



{£{x,t) + ^£{x,t) 



1 < 



np{x,t) 
np(x,t) 



(£{x,t)-^£{x,tY 



1 > 



Defining an elliptic parameter m — m{x,t) by 

m{x, t) — mR{x, t) := 
and a real phase $(a;, i) by 



l+£{x,t) 



\ (g(o) - yG(oy 

\ (g(0) + v/G(0)2-4 
e(0,l), (x,i)e5R, 



> 0. 



uj{x,t')dt', 



where 



uj{x, t) 



^£{x, t) + ^£{x, t)2 - 1 + ^£{x, t) - v/?OMF 



2K{m{x,t)) 



y/1 - np{x,ty 



(1.67) 



(1.68) 



(1.69) 



(1.70) 



(1.71) 



the following asymptotic formulae hold pointwise for (x, t) G S'r 



cos 



UM{x,t) 
UN{x,t) 



/'2<S>{x,t)K{m{x,t)) 

\ TreAT 

/2^{x,t)K{m{x,t)) 
m : 



m{x,t)] +0{eN) 



m{x,t) +0(ejv) 



(1.72) 



duN , , 

^N^{x^t) 



4K(m(x,t))9$ (2^{x,t)K{m{x,t)) 



dt 



dn 



]m{x,t) +C'(ejv). 



Moreover, the error terms are uniform for {x,t) in compact subsets of Sr. The phase ^{x,t) also satisfies 
d^/dx = k{x, t) := uj{x, t)np{x, t). 



The accuracy of the asymptotic formulae (1.72) is illustrated in Figure 1.5 
We now make several observations about these results: 



The asymptotic formulae (1.651 correspond to a high-frequency superluminal librational wavetrain 
that is relatively slowly modulated through the {x, i)-dependence of the quantities rip and £. Indeed, 
Taylor expansion of these formulae about a fixed point (a;o, to) e shows that if one sets x = xq + 
cnx and t = to+e^t and takes the limit e^r J, holding x and t fixed, we recover an exact superluminal 
librational wavetrain solution of the sine-Gordon equation as a function of x and t characterized by 
the fixed energy £{xQ,tQ) and the linear phase <i>(a;o,<o) + k{xQ,tQ){x ~ xq) ~ w(xo,io)(^ ~ ^o)- It is 
clear in this case that the function U]\[{x,t) is confined to the range (— tt, vr) and that both un and 
its time derivative are periodic functions of the linear phase as is consistent with librational motion. 
Similarly, the asymptotic formulae (1.72) correspond to a high-frequency and slowly modulated 
superluminal rotational wavetrain, and in particular u^v is monotonic while its derivative is periodic, 
as is consistent with rotational motion. 

We have excluded the curves t = t± (x) from S'r as a matter of complete academic honesty, as our 
proof would require a technical modification to extend pointwise asymptotics to these curves, and to 
have uniformity for t w t±{x) is yet a further matter requiring a double-scaling limit. However, it is 
easy to see that the explicit terms in the asymptotic formulae (1.72) undergo no phase transition in 
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Figure 1.5. The cosine of the exact solution UN[x,t) for the special initial data described 
in l.'1.2 for A — i/A plotted with points, and the corresponding asymptotic formula plotted 
with curves, for fixed x = —0.15625 and Q <t < 0.5 so that (x,t) S Sr yielding rotational 
motion as described by Theorem Left: N = 8, or equivalently ejy = 0.09375. Right: 

N = 16, or equivalently epj — 0.046875. Note that for the given value of x, t_.{x) w 0.153, 
but that there is no evidence of any change in behavior near this value of t in either the 
exact solution or the asymptotic solution. 



the vicinity of the curves t = t±{x), and while we cannot honestly exclude the possibility that there 
is some phenomenon here to be captured by more detailed analysis, there is no indication of such in 
the plots shown in §1.2| 

As part of our proof, we show that in each case the asymptotic formula for e^duf^ / dt is consistent 
with those for cos(iujv) and sin(^MAr) in the sense that differentiation of the latter with respect to 
t assuming that the error terms remain subdominant after differentiation yields the former up to 
terms of order O(ejv). 

The inequalities (1.61) and (1.661 together with the initial condition np(a;,0) = give information 
about the direction of motion of the waves. For example, if a; > Xcvit and t > 0, the fluxon condensate 
behaves like a train of superluminal librational waves propagating rapidly to the left (since the phase 
velocity Up is large and negative), and if < a; < Xcrit and t > Q the condensate behaves like a train 
of superluminal rotational waves propagating rapidly to the right. 

The fields np{x, t) and £{x, t) have in each case exactly the interpretation of reciprocal phase velocity 
and energy as explained earlier in the context of formal modulation theory. Although we conclude 
that these quantities satisfy the Whitham system (1.32), we wish to emphasize that at no point 
in our proof do we apply techniques from the theory of partial differential equations {e.g., the 
Cauchy-Kovalevskaya method) to solve the Whitham system with the specified initial data. Instead, 
the fields np{x,t) and £(x,t) are constructed by means of the solution of a system of nonlinear 

We first obtain functions p — p{x,t) and 
4.7), or for (a;, t) G Sr the equations 



algebraic equations via the Implicit Function Theorem 
q = c\{x,t) by solving the equa tion s M — I = (see Proposition 
M — I — (see Proposition 4.9), and then recover np{x,t) and £{x,t) therefrom via (4.35) and 
(4.92) respectively. The fact that these functions satisfy the quasilinear partial differential equations 
(1.32) is, from the point of view of our methodology, more or less a coincidence. 
Interestingly, the asymptotics of the fiuxon condensate are accurately described by superluminal 
wavetrains, even though the fundamental particles of the condensate are all subluminal solitons of 
kink and breather type. 

The semiclassical asymptotics of the sine-Gordon equation are, at least in the case of a sufficiently 
strong initial impulse consistent with Assumption |1.6[ much more complicated even for small times 
t (independent of ejv) than in the case of more well-known integrable partial differential equations 
like Korteweg-de Vries (see [B]) and focusing nonlinear Schrodinger (see [HI [53]). For these latter 
problems, there is an initial stage of the evolution where the asymptotics are described for all a; G M 
in terms of elementary functions. Here, we require two different types of formulae both involving 
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higher transcendental functions, and even these do not suffice to describe the dynamics near the 
transitional points x — ixcrit- 



If, in contrast to Assumption 1.6 we instead take —2 < G(0) < 0, then it will be clear from our 
proof that the region S\_ extends to include a full neighborhood of (x, t) — (0, 0), and hence becomes 
connected, including the entire t = axis. Then Theorem [T]T] will suffice to describe the semiclassical 
asymptotics for small time uniformly for a; e M. This is consistent with the dynamics pictured in 
Figure 

Our results show that the two superluminal types of modulated single-phase wave behavior observed 
in Figure [L2| for small time and x bounded away from ixcrit are universal for the class of initial data 
we consider. It is reasonable to conjecture that the qualitatively different behavior observed in the 
space-time plane on the other side of nonlinear caustic curves might also be universal, and that it 
might correspond to various types of modulated multiphase waves. We note that, with appropriate 
modifications and additional work, the methods described in this paper are capable of handling these 
other cases. Specifically, for (x, t) in these regions the model Riemann-Hilbert problems 
refrhp:wOdotlibrational and |5.2| will need to be generalized to ones that are solved using the function 
theory of higher-genus hyperelliptic Riemann surfaces. The associated Whitham equations that 
generalize ( |1.32 1 will also have a correspondingly increased number of dependent variables (see 

[nmniin]). 



The relevance of the fluxon condensate {um{x, t)} to the Cauchy problem for ( 1.1 1 with e ~ ej^ and with 



pure- impulse initial data characterized by the even Klaus-Shaw function G is then the following result: 

Corollary 1.1. When t = Q, the fluxon condensate {ujv(2;,t)} associated with the pure-impulse initial 
condition of impulse profile G(-) satisfies 

UAr(x,0) = ©(ejv) (mod47r) and eN^^ix,0) = G{x) + 0{eN) (1-73) 

where the error estimates are valid pointwise for x ^ and \x\ ^ a^criti ^.''^d uniformly on compact subsets of 
the set of pointwise validity. 



In this sense, the fluxon condensate approximates the solution of the Cauchy problem for ( 1.1 ) with initial 



data (1.2 1 when e = and N is large. 



1.4. Outline of the rest of the paper. The remaining sections of the paper are devoted to the proofs 
of Theorems |1.1| (small-time librational asymptotics) and 1.2 (small-time rotational asymptotics). These 



are proven using the well-developed inverse-scattering method (see the brief discussion in |1.1| and the more 
detailed exposition in our paper [i]). The fluxon condensates {uN{x,t)} we study are deflned by their 
scattering data, effectively skipping the forward-scattering transform. Thus all of our analysis concerns the 
inverse-scattering transform. 

Our approach is to use the Riemann-Hilbert problem formulation of the inverse-scattering transform. 
Since the scattering data are reflectionless (that is, comprised of only eigenvalues and the corresponding 
proportionality constants), the associated Riemann-Hilbert problem for the 2x2 matrix- valued function 
J(u') has only poles and a completely trivial jump on the positive real axis (which could be removed at the 
cost of artificially doubling the number of poles through the transformation w = z^). This setup and parts of 
the subsequent analysis are similar to an analagous work on semiclassical soliton ensembles for the focusing 
nonlinear Schrodinger equation [Mj. 

Our first step is to make a local change of variables J{w) — > M(u;) in the Riemann-Hilbert problem in 



5 3.2 that removes the poles at the price of introducing further jump contours (which are more amenable to 



analysis). Depending on the value of x and t, different transformations are used in different parts of the 



spectral w-plane. These different transformations are illustrated in Figures 3.1 3.6 Note that only one of 



these cases (A = shown in Figure 3.1 1 is required to analyze solutions with even initial data for small times 
and X bounded away from the origin. 

In |3.4| we make another transformation M(z/;) — > N(it;) involving a g-function, a standard tool (first 
introduced in [6 ) for controlling the behavior of jump matrices in the Riemann-Hilbert problem. The details 
of the construction of the g-function are contained in Q Finding the (^-function involves identifying parts 
of the jump contours as a band /3 and parts as a gap 7. There are two possible topological configurations 
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for the band /3 (see Figures [4.1| and 4.2). In case L (leading to librational wavetrains), the band endpoints 
are a complex-conjugate pair (w = p ± i^/—q). In case R (leading to rotational wavetrains), the band 
endpoints are real [w — p ± y^). The band endpoints are chosen to enforce certain conditions making the 
jump matrices easy to analyze in the limit epf ^ 0. The key part of the analysis is proving there exists a 
(/-function independent of eN such that these conditions are satisfied for some fixed small time. This analysis 
is particularly delicate in a neighborhood of {x,t) — (0,0), as described in 5 4.3.2 In 5 4.4 we show that the 
band endpoints associated to an admissable g- function satisfy the Whitham modulation equations (1.32) in 
Riemann invariant (diagonal) form. 



In 5 5.1 we open lenses in the Riemann-Hilbert problem by making the transformation N(u;) — !■ 0{w) 
as illustrated in Figures |5.1H5.4[ This has the effect of making the jump matrices exponentially close to 
the identity in e/v except on the band /3, the positive real axis, and in small neighborhoods of the band 
endpoints. In 5 5.2 we construct the global parametrix 0{w), which is an explicit approximation of 0{w) 



that is valid in the whole complex w-plane. In 5 5.3 we prove rigorously that 0{'w) is an C'(ejv)-approximation 
of 0{w). Thus we can use 0{w) to compute the solution un{x, t) to the sine-Gordon equation modulo errors 
of size O(eAr). Finally, Appendix [b| contains the details of calculating the solutions UN{x,t) from the outer 
parametrix 0(w). A significant portion of the effort in Appendix [b] is devoted to translating more or less 
standard formulae involving Riemann theta functions of genus one into expressions involving Jacobi elliptic 
functions with suitable moduli. This is necessary to establish the simple form of the asymptotic formulae 
appearing in the statements of Theorems |1.1| and |1.2| 



1.5. Notation and terminology. With the exception of the identity matrix I and the three Pauli matrices, 

(1.74) 
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we will denote matrices by bold capital letters {e.g. M) and vectors by bold lowercase letters {e.g. v). We 
will use A for the closure of a set ^ C M^, and denote complex conjugation with an asterisk. 

In what follows, by a planar arc we will mean the image of a continuous and piecewise-smooth one-to-one 
map (0, 1) 9 1 1— ?> w{t) S C = with parameter t and nonvanishing derivative. By a planar contour we will 
mean the M^-closure of a finite union of arcs. Thus a contour is always a closed set in the topology of and 
may contain self-intersection points where various arcs meet. An oriented contour is a contour K written 
in a particular way as the M^-closure of a finite union, denoted K, of pairwise-disjoint arcs each of which is 
assigned an orientation in the obvious way according to its parametrization by t. An oriented contour may 
include at most a finite number of points at which the orientation is not properly defined. These may be 
self-intersection points, endpoints, or points dividing an arc into oppositely-oriented sub-arcs. 

If K is an oriented contour and / : C \ if — >■ C is an analytic function, we denote by /+(^) (respectively 
f -{£,)) the boundary value taken by f{w) as w ^ (, E K from the left (right) according to local orientation, 
if it exists. We use analogous notation for vector-valued functions {e.g. v±(^)) and matrix-valued functions 
{e.g. M-i-(^)). Given a contour K we define a metric on C \ K as follows: 

dK{wi,W2) inf length(P) (1-75) 

PCC\K 

where the infimum is taken over paths P d C \ K connecting the points w and z. If < a < 1, an analytic 
function / : C \ if — )• C is said to be uniformly Hdlder-a continuous if 

\f{wi)^f{w2)\ 



sup 



< oo. 



(1.76) 



This definition generalizes in the obvious way to vector-valued or matrix-valued functions, and also may be 
analogously defined relative to open domains U\ K, U C C, U open. 

2. Formulation of the Inverse Problem for Fluxon Condensates 

[ — > M satisfying the conditions set out in the introduction be given 



a 
on 



Let a function G : M - 
decreasing sequence {eat}?^ 
the positive imaginary axis by the Bohr-Sommerfeld quantization rule (1.42) with the WKB phase integral 
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and let 

be defined by (1.54). For each integer N > 0, define numbers A° — X% 



(1.41) and the conditions e — cm and N{e) = N. Regardless of the value of = 1,2,3, ... , the numbers 
{''^% k}k=o ^^^^ each have exactly two distinct preimages in the w-plane with | arg(— w)! < vr under the map 
A = E{w) if the following additional condition is satisfied. 

Assumption 2.1. The fraction 'i'{i/2)/\\G\\i is irrational. 

This simply guarantees that none of the numbers k}k=o coincide with A = i/2, the unique critical 
value of E{w) for | arg(— w)! < tt, which ensures the poles that will appear in the definition of Riemann- 



Hilbert Problem 2.1 below are all simple. Rationality of 4'(i/2)/||G'|| i can also be admitted with the same 
effect at the cost of passing to a subsequence of values of N . 
Let a function Q{w) be defined for | arg(— it;)| < tt as follows: 

Qiw) = Q{w] X, t) := E{w)x + D{w)t, (2.1) 

and define 

E{w) + X% . 

k=0 ^N,k 



According to Assumption 2.1 the denominator of IIn{w) has 27V simple zeros in the domain of definition, 
while the numerator of IIn{w) is analytic and nonvanishing. Let P/v denote the set of poles of IlNiw), '2N 
points consisting of 2Nb points in complex conjugate pairs on the unit circle and 2A"k points on the 
negative real axis in involution with respect to the map w — >■ 1/w. We may write Tlj^j^w) equivalently in the 
form 



U^{w) = TT ^-^^y (2.3) 




The basic Riemann-Hilbert problem of inverse scattering to construct the fluxon condensate corresponding 
to G is then the following. We are following the description of the inverse-scattering problem given in 
Appendix A of [4] but we are exploiting the symmetry 2 i— — z to formulate the problem in terms of the 
complex variable w = z'^, which introduces a jump on M+. 

Riemann-Hilbert Problem 2.1 (Basic Problem of Inverse Scattering). Find a 2 x 2 matrix function 
'H.{w) = Hn{w; x,t) of the complex variable w with the following properties: 

Analyticity: H(z«) is analytic for w e C \ {Pn U M+). 

Jump Condition: There is a neighborhood U — Un of M+ such that ll{w) is uniformly Holder-a 
continuous for all a G (0, 1] on U\R^ . Letting be oriented from left to right, the boundary values 
taken by H(w) on M_|_ are related by the jump condition 

H+(0-a2H_(0a2, (2.4) 

Singularities: Each of the points of Pn is a simple pole o/H(w). If y G Pn with E{y) = A^; ^ for 
k = 0,. . . ,N - 1, then 



Res H(w) = lim H(u;) 





(-l)'^+iRese2*'5("''^'*)/'"nAr(w) 



(2.5) 



These amount to one matrix-valued condition on the residue of H(w) at each of its poles. 
Normalization: The following normalization condition holds: 

lim H(w) = I, (2.6) 

where the limit is uniform with respect to angle for \ arg(— ti;)| < tt. 

It is an easy application of Liouville's Theorem that any solution of this problem must satisfy det(H(w)) = 
1, from which it follows that if Hi{w) and H2(t«) are any two solutions, the matrix ratio R(ti;) :— 
Hi(w)H2(w)~^ is analytic for |arg(— w)! < tt (has removable singularities at the points of Pn and ac- 
quires no additional singularities from inversion of H2(u')), is Holder-a continuous in U \ and tends to 
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satisfy R+(0 = cr2R-(C)cr2. If 



S(z) 



(2.7) 



the identity matrix as w — >■ oo. The boundary values taken by R(k7) on 
we set 

'R(z2), > 0, 

^a2R{z^)cT2, < 0, 

then it is clear that S{z) extends to an entire function of z e C with identity asymptotics as z — !• cxd. 
Hence Liouvillc's Theorem shows that S(z) = I and so by restriction to 3{z} > 0, Hi(ii;) = H2(iy) for 
I arg(— w)| < TT. Therefore, solutions to Riemann-Hilbert Pro blem 2.1 are necessarily unique if they exist. 

Given any solution H(w) of Riemann-Hilbert Problem 2.1 another is easily generated by setting 11^{'w) :— 
II{'w*)* . By uniqueness it follows that H''(w) = H(w), or equivalently, that the unique solution of Riemann- 
Hilbert Problem |2. 1| necessarily satisfies 

H(u;*) =H(u;)*. (2.8) 

If for some {x,t) e 



2.1 



for all integer N > No 

for sufRciently large Nq, then a construction similar to (2.7) involving the variable z such that w = shows 



H(ii') is the unique solution to Riemann-Hilbert Problem 



that if it exists the solution H(w) has convergent series expansions of the form 

oo 

Il{w)=J2il°N''i^,t)iV=wf, \w\<r 



fc=0 



and 



H(u;) = I + ^H^•'=(x,^)(^/=^)- 



> R 



(2.9) 



(2.10) 



fc=i 



for suitable numbers r and R independent of N. We will use the notation 

AM{x,t):^U%°{x,t), B%ix,t):=li'}f{x,tr'U%\x,t), B^{x,t) -.^ ll^'\x,t), (2.11) 

defining three matrices depending parametrically on {x,t) e and the integer N > Nq. These matrices 
necessarily satisfy the conditions 

det{AN{x,t)) = 1, Apfixjt) ^ a2AN{x,t)(T2, and AN{x,t) = AN{x,t)* , (2-12) 

tr(B^(x,i)) = 0, B^{x,t)^-a2B^{x,t)a2, and B^ (x^t) = B^ (x^t)* , (2.13) 

and 

tv{B%{x,t)) = 0, B%{x,t) = -a2B%{x,t)a2, and B%{x,t) = B%{x,t)* . (2.14) 

has a 



2.1 



Definition 2.1 (Fluxon condensates). Given {x,t) G K^, suppose that Riemann-Hilbert Problem 
solution for all integer N > Nq. Then, the fluxon condensate {uAr(x, t)}^^^^ associated with the impulse 
profile G is given (modulo addition of arbitrary integer multiples of Air) in terms of Aisi{x,t) as follows: 



cos I -UN{x,t) ) = AN^ii{x,t) and sin 



1 



UNix,t) = AN,2lix,t). 



(2.15) 



Note that the identities (2.121 ensure the the reality of these expressions as well as the Pythagorean identity 
sin{^UN{x, t))^ + cos{^U]y{x, t))^ 



1. 

Proposition 2.1. Suppose that C M'^ is open and an integer Nq > is given such that Riemann-Hilbert 
Problem 2.1 has a solution whenever (x, t) Cz Q and N > Nq, and that the fluxon condensate {upf{x, t)}^^j^ 
is defined as above for {x,t) € il. Then for each integer N > Nq, u = Uiq{x,t) is an exact real-valued 

(2.16) 



solution of the sine-Gordon equation (1.1) with e = €n- Moreover, 



du 



ejv^(x, t) = i?^,i2(x, t) + B^,i2(a^, t)- 



Proof. Consider the matrix F(w) defined in terms of the solution H(w) of Riemann-Hilbert Problem 2.1 by 

F(u;) := H(w)e-*^("''^'*)''='/"". (2.17) 

Aside from an introduced essential singularity at w = 0, the matrix Y{w) is analytic exactly where H(w) 
is, and has similar properties where analyticity is violated. Namely, F(iy) is Holder-a continuous inU \ M+ 
except at w = where the exponential factor has an essential singularity, and as a consequence of the 
identity Q+(C) = ^Q-(C) for ^ € M+ one has that the boundary values of F(w) satisfy the jump condition 
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F+(^) = (T2r_(^)(T2 for ^ G M+. Also, F{w) has simple poles at the points of Pjv, and if y G Pn with 
E{y) = >^°N,k then 

ol 

(-1)^+1 Res HjvH • (2.18) 



ResF(w;) = lim F(u;) 



Note that neither the jump nor residue conditions involve {x, t) G M^, from which it follows that the matrices 
UH := AieNF^{w)F{w)-^ = 'iieNllx{w)U{w)-^ + 4E{w)H{w)a3H.{w)-^ (2.19) 

and 

V{w) := 4ieNFt{w)F{w)-'^ = 4ieA,Ht(M;)H(w;)-i + 4D(w)H(w)ct3H(w)-i (2.20) 
are functions of 2; = i\/—w that are analytic for z G C \ {0}. These functions have the following asymptotic 



behavior as 



and u; — ^ 00: 



U(w) 



1 



-.iANix,t)a3AN{x,t) ^ + 0(1), w^O, 



-wias + i[B'^{x,t), as] + o{l), 
1 



w — >■ CX), 



(2.21) 



V(w) 



-w 



ziAN{x,t)a3AN{x,t) ^+0(1), w ^ 0, 



^wicr3 + i[B'^{x,t),<T3] + o{l), 



w 00. 



It then follows from Liouville's Theorem applied in the z-plane that in fact U(u') and V(u;) are Laurent 
polynomials in z of degree (1, 1): 



U(w) = y^ia3 + i[B'^{x,t),a3] + —^iAN{x,t)a3AN{x,t) ^ 

1 



Y{w) = V^i(^3 + i[BN{x,t),a3] 



ziANix,t)a3AN{x,t) ^ 



According to (2.12), we may write A^ixji) in the form 

A (-r A - |"cos(0) - sin(</>) 
AN(X,t) - ^^.^^^^ C0S((/.) 

where — 0Ar(a;, t) is a real angle. Likewise, according to ( 2.13[ ) we may write B^(a;, t) in the form 

B^(x,t)=[^^ d 

where c = CNix,t) and d — dAr(x,t) are real- valued fields. In terms of these we therefore have 

1 



U(it;) = y/—wia3 + 2ca2 H (i cos(2(/))o'3 + ism{2(/))ai) 

\/—w 



V(w) 



-wia3 + 2c(T2 



{i cos{2(f))a3 + isin(2(?!))cri) . 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



Since when considered as a function of x and t, the matrix F^w) is a simultaneous fundamental solution 
matrix of the first-order overdetermined system 

AieNF^ = UF and AieNFt = VF, (2.26) 

it follows that this system is compatible, implying that the zero- curvature condition 

4iejvUt - 4^e^^V, + [U, V] = (2.27) 

holds. Since {cti/ \/—w, (73/^— w, (T2} is a linearly independent set, the zero-curvature equation implies the 
following three equations: 

(eAr(/)t -I- eN^ix - c) cos{2(j)) = 
(ew^t + eN(f>x - c) sm{2(j>) = 

£NCt - eNCx + ^ sin(20) = 0. 
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(2.28) 



The first two are together equivalent to the single equation e^i't 
eliminate c from the third equation, yielding 



^NVxx 



- sin(2(/.) - 0, 



— c = 0, which may be used to 



(2.29) 



which is the sine-Gordon equation (1.1) for ujsi — 2(j) with e — e^. 
Obviously we have 



UN 



-UN 



sin(0) : 
cos((/)) 



ANai{x,t) 
■ AN.\\{x,t) 



(2.30) 



in accordance with (2.23) and (2.15). Finally, the identity (2.16) arises from the constant term in the 
Laurent expansion (in powers of \/ —w) of the differential equation 4iejvHi(i(;) -|-4£'(w)H(u')ct3 — V(w)H(w) 
equivalent to 4ieArFt(w) ~ V(ti;)F(w). □ 



Note that since Pn consists of points y with | arg(— ?;)| < tt, it is a consequence of the definitions (1.39) 



and (2.1 ) involving the principal branch of the square root that for each fixed iV and t we have 

0. 



lim Res 

>+oo w=y^Pt 



(2.31) 



From this it can be shown that for each u; G C \ {Pn U M+) and for each i G M, 'H.N{w\x,t) — > I as 
X — >■ +00. By contrast, the asymptotic behavior of H as x — > —00 is not at all clear, and this is more 
than merely a technical difhculty. The "preference" that H has for large positive x stems from a choice in 
defining the spectral theory for which Riemann-Hilbert Problem 2.1 is the inverse-spectral problem in terms 
of scattering "from the right" . Given this asymmetry, it should be no surprise if asymptotic analysis of H 
in the semiclassical limit of large must take a different route for — a; than it does for x. This asymmetry 
can be addressed in two different ways: 

• If the initial conditions F{x) := u{x, 0) and G{x) := eNUt{x, 0) for the sine-Gordon Cauchy problem 
are both cither even or odd functions then this symmetry is preseved in time and knowledge of the 
solution for, say, positive x is sufficient. The special initial conditions under consideration in this 
paper have even symmetry. Moreover, it can be shown directly from the conditions of Riemann- 
that the functions {uN{x^t)}'^^j^^^ 
functions of x. 



Hilbert Problem 



2.1 



approximating the solution to the Cauchy 

problem are all even 

• The inverse problem may be reformulated in terms of scattering "from the left" . 

In fact it will turn out that for certain values of {x, t) S of interest, neither of the above two approaches 
will be of much help. However, a modification of the second approach that can be thought of as formulating 
the inverse problem in terms of scattering theory "partly from the left and partly from the right" will indeed 
succeed. The idea is to select a subset A C Pn of the pole divisor of 'ii{w) and set V := Pn \ A. Then 



define from the solution Hjv(w; t) of Riemann-Hilbert Problem 2.1a new matrix function given by 



3 Niw^Xjt) :=Yi.N{w;x,t) 




(2.32) 



It is easy to see that if HAr(w;a;,t) satisfies Riemann-Hilbert Problem 2.1 then 3 Niw^x^t) satisfies the 
following equivalent problem. 

Riemann-Hilbert Problem 2.2 (Modified Problem of Inverse Scattering). Find a 2 x 2 matrix function 
3{w) — 3 N{w;x,t) of the complex variable w with the following properties: 
Analyticity: 3{w) is analytic for w G €. \ {Pn U M.^). 

Jump Condition: There is a neighborhood U — Un 0/ IR+ such that J(w) is uniformly Holder-a 
continuous for all a G (0, 1] on [/\ . Letting be oriented from left to right, the boundary values 
taken by 3{w) on IR+ are related by the jump condition 



J+(0 = ^2J-(e)fT2, 



(2.33) 
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Singularities: Each of the points of Pn is a simple pole of 3{w). //y G V C Pn with Eiij) = j. 
for k — 0, . . . , N — 1, then 



Res J(u') ~ lim 3{w) 



w=y 



0' 

Res e2^'3("''^'*)/""njvM 



w=y 



id if y e A C Pn with E{y) = X% f. for k ^ 0, . . . , N - 1 , then 



Res J(w) = lim J{w) 

w—y w^y 

where IIn{w) is re-defined as 



Res e~2*^('"'="'')/'"nArH-i' 

w=y 





TiNiw) := n 




(2.34) 



(2.35) 



(2.36) 



qeA 



(note that this definition reduces to the earlier one, see (2.2) and (2.3), in the special case when 
A = and hence V = Pn)- These amount to one matrix-valued condition on the residue of3{w) at 
each of its poles. 

Normalization: The following normalization condition holds: 



lim J{w) — I, 

■u;— J-oo 

where the limit is uniform with respect to angle for \ arg(— ti;)| < tt. 



(2.37) 



It is clear that in passing from H to J the nature of the residue conditions is changed near those points 
y & A C Pn and is left unchanged near those points y € V = Pn \ A. The special case of A = Pn and V = 
corresponds to reformulating the inverse problem in terms of scattering theory "from the left" . Indeed, we 
see that in this case the exponential e^*'^^'"'^'*)/'^™ has been completely replaced by its reciprocal, which 
makes the limit x — > — oo particularly transparent; from the conditions of this problem it is easy to prove 
that if A = Pn then 3 n {w; x , t) I as x ~oo whenever t G M. and w G C\ (P/v UK+) . More generally, this 
calculation suggests that if a; € M is a value for which semiclassical asymptotic analysis of ii{w) is difficult, 
the problem may be resolved by analyzing the equivalent matrix 3{w) instead for a particular choice of A. 

From now on we will take Riemann-Hilbert Problem |2.2| as the basic object of study in the limit iV — ^ oo, 
where the set A is to be chosen differently for different {x, t) e to facilitate the asymptotic analysis. 

Having formulated Riemann-Hilbert Problem |2.2[ we can now explain how the plots in §1.2| were made. 
For any choice of A C P/v it is easy to see that 3{w) is a rational function of z = i\J—w and it may therefore 
be written as a finite partial fraction expansion with simple denominators and constant term I to satisfy 
the normalization condition. The matrix coefficients in the expansion are then determined from the residue 
conditions ( 2.34 ) and (2.35 1, which imply a square system of linear equations to be solved for the coefficients. 



In this way, the construction of sin(|ujv(2;, t)) and cos(|uAr(x, t)) may be reduced to a finite-dimensional 
(of dimension proportional to iV) linear algebra problem parametrized explicitly by (x,t^ G M^. The linear 
algebra problem is ill-conditioned when is large (this difficulty can be partly ameliorated by judicious 
choice of A given x and t), but nonetheless implementing this approach numerically allows one to explore 
the phenomenology of the semiclassical limit while avoiding many traditional pitfalls (for example, stiffness 
and propagation of errors) of direct numerical simulation of the Cauchy problem for the sine-Gordon equation 



(1.1) when e is small. For more details about implementation of this direct approach to inverse scattering 



see our recent paper 0]. 

3. Elementary Transformations of 3(vj) 
We now embark upon a sequence of explicit invertible transformations with the aim of converting Riemann- 



Hilbert Problem 2.2 into an equivalent one that is better suited to asymptotic analysis in the limit TV — > oo. 



Throughout the rest of the paper we will use the following notation for the composition of the WKB phase 
integral ^' with the function E: 

B^{w) ^{E{w)). (3.1) 



3.1. Choice of A. The set P^v can be decomposed into a disjoint union Pat = U with P^ n P^ = 0. 
Here P^ consists of Nq nonreal complex-conjugate pairs of complex numbers on the unit circle in the 
w-plane, while P^ consists of TVk pairs of reciprocal negative real numbers {i.e. of the form {w, 1/w) with 
w < 0), none equal to —1. Since each conjugate pair of poles contributes to the fluxon condensate one 
breather soliton while each negative pole contributes one kink soliton, the fluxon condensate may be viewed 
as a nonlinear superposition of 2A^k kinks and Nq breathers. As a consequence of Assumption |1.6[ both A^k 
and A"b are proportional to N and as A^ ^ oo, P^ fills out the whole unit circle while P^ fills out a negative 
interval of the form [a, bl where: 



1 



:=-4(v/G(0)2-4-G(0) 



:=-J(v/G(0)2-4 + G(0))' = i, 



(3.2) 



Note that both a and b are independent of A^ and that a < — 1 < b < 0. If Assumption 1 1 .6| were not satisfied, 
that is, if G(0) > —2, then P^ would be empty and P^ would fill out a proper sub-arc of the unit circle as 
N ^oo. 

We will consider six different configurations for the subset A C Pn- To each fixed number Tqo G (o, — 1) U 
(—1, b) we associate a sequence {tn}n'^No '^^ ^^^^ numbers with limit Tqo by 



9o(roo) 



N > Nq. 



(3.3) 



(This equation can be solved for r^y near Tqc for A^ > Nq with Nq sufficiently large by the Implicit Function 
Theorem since the only critical point of 9q{w) in (a, b), namely w = —1, has been excluded.) According to 
the Bohr-Sonimerfeld quantization rule ( 1.42[ ), when A is large, r^r is approximately halfway between two 
neighboring points in P^. The six choices of A we consider are the following: 



A = . In this case J(w) = H(w) and Riemann-Hilbert Problems 2.1 and 2.2 coincide 
A = P^^. In this case we choose a < < 



near w 
A 



a. 



1 and set A = P^^ 



-I < Too < b and set A 



P^r^[a,TN] 



Thus A is localized 



P^ n [Tjv,b]. Thus A is 



-1 and set V = P^^ 



:=P^n[a,TN] 



Thus V is localized 



In this case we choose 
localized near w = b. 

V = . This is complementary to the case when A — 
. In this case we choose a < t^o < 
a. 

In this case we choose —l<T^<b and set 
b. 

oo as A^ — > cxd) as a transition point. One consequence of 



V - P^^ 
near w = 
V 



pK>- 



Pn n [tnM- Thus V 



localized near w 
We will refer to tm (and sometimes also its limit t< 



the definition ( |3.3[ ) in relation to the Bohr-Sommerfeld quantization rule ( 1.42 ) and the fact that Pn contains 
an even number of points is the identity 

00 (tat) 



= #A (mod 2) 



(3.4) 



where #A denotes the number of points in A. 

For most (a;,t), specifically outside of a small neighborhood of {x,t) — (0,0), the only cases we will use 
are A = and V = 0. In fact, A = is the only case needed to analyze the behavior for small times and 
positive X bounded away from the origin. It is then possible to use evenness of the Cauchy data (Assumption 
1.3) to obtain results for negative x. However, we include the case V = for completeness (and to admit 
future generalizations of our methodology to non-even Cauchy data). The analysis of the neighborhood of 
the origin that relies on the four additional cases in which neither A nor V is empty is somewhat more 
complicated and is presented in a self-contained fashion in 5 4.3.2 A reader not interested in these details 
can safely ignore all of the cases except for A = and V = along with references to any transition points, 
and may also skip most of §4.3.2| 

3.2. Interpolation of residues. Removal of poles. As a consequence of the Bohr-Sommerfeld quanti- 



zation rule (1.42), whenever y E Pn with E{y) = j, then 



{-If 



,N-l. 



(3.5) 



Thus the exponential function on the left-hand side analytically interpolates the signs (—1)'^ at the corre- 
sponding poles of J. This gives two different ways to combine the sign (—1)'^+'^ appearing in the nilpotent 
residue matrices in the conditions (2.34) and (2.35 1 with the residue factor, a fact we can use to formulate an 
invertible transformation of J into another matrix M that will not have any point singularities whatsoever. 




Figure 3.1. The case o/ A = 0. Here = 0. The set Poo := [a, b] U in which Pm 
accumulates for large N is shown in dark red for reference, and the branch cut M+ is shown 
with a black line. 




Figure 3.2. The case of A = ■ The regions and abut the real axis in the 
interval (a, tat) where a < Too < min(u)+, —1). 

We now introduce four disjoint open regions of the w-plane with | arg(— w)| < tt, denoted and fi^, 
such that E{^1) = D_|_ U where := fij U fiY U fi^ U O^, and where the open rectangles D± are as 
defined following the statement of Proposition [T]2j The details of the definitions of these regions are different 



depending on which of the six cases of choice of A we are considering (see Figures 3.1 3.6 1, but the following 

features are common to all cases: 

• f2 is independent of N , x, and t, and always contains the pole locus P/v for all N. Moreover, fi^ 
contains V C Pjv while fl^ contains A C P/v, where fl^ --QjuflY and fl'^ := flf U fl^ . This will 
allow us to remove the poles from J by making appropriate substitutions based on (3.5) in i7j and 
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Figure 3.3. The case of A ^ P^'^ . The regions and abut the real axis in the 
interval (r^r, b) where max(u'^, —1) < Tqo < b. 




Figure 3.4. The case ofW ^9. Here = 0. 



• Schwartz symmetry is present: D,Y — £7^* and V,^ — Q^* . 

• Either the common bomidary of fij with SlY or that of fl^ with fi^ (but not both) contains a 
Schwartz-symmetric closed curve that meets the real axis at w = 1 and exactly one other point, 
w = ui^ e (a, b). This closed curve is allowed to cross the unit circle at points other than w = 1, 
perhaps nontangentially. 

• Only the common boundary of fl^ or the common boundary of fl^ and flY, will depend 
on N (necessarily through the transition point tjv where these curves meet the real axis); the other 
curves must be independent of N. 

For convenience, we assume that the common boundary of and 51^, and also the common boundary of 
and 5lY, are arcs of level curves of Sy{i?(w)} as illustrated in Figures 
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3.1 



3.6 



Figure 3.5. The case of V ^ P^^ . The regions f2y and ilY. abut the real axis in the 
interval (OjTat) where a < Too < min(u'^, 1). 




Figure 3.6. The case of V = P^^ . The regions 
interval (tjv, b) where max(u'"'", 1) < Too < b. 



nY and ilY abut the real 



axis m the 



In each of the six cases we wiU now transform the matrix J(u') into an equivalent matrix M(w) 
Mjv(w; x, by the following explicit formula: 



m{w) 



3{w) 



3{w) 



1 



1 ±injv(w)~-^e~t'^*'5*^"'''^'*^'''*^"^'"^l/'^" 
1 



3{w), w e c\ (r2uE+). 

Note that by this transformation, M inherits from J the Schwartz symmetry M(u'*) — M(u)) 



(3.6) 



It is a consequence of the interpolation formula (3.5) and the residue conditions (2.34) and (2.351 that 



M(ti;) has only removable singularities in and up to the boundary of each of the regions of its definition, 
and thus it may be considered as a sectionally analytic function of w taking continuous boundary values 
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on but having jump discontinuities across an oriented contour E that is the positively-oriented boundary of 
U together with the negatively-oriented boundary of JlY U and two intervals of the positive real 
axis, which we take to be both oriented to the right. 

We introduce notation for various subcontours of E as follows: 

• The part of E consisting of the common boundaries of fi^ (respectively omitting only the 
segments on the positive real axis we denote by E^ (respectively E^). 

• The common boundary of and will be denoted E^^. The common boundary of fl^ and flY 
will be denoted E^"^. _ 

• The common boundary of (respectively O^) and O will be denoted E^ (respectively E^). These 
are just the non-real arcs of E \ (E^ U E^ U E"^^ U E^^). 

• The part of the common boundary of the sets fij (respectively $7^) on the positive real line will be 



denoted E^g 



^respectively E^q). 



• The contour M+ \ (E^q U E^q) will be denoted simply E>o. 



The contour E and its components are illustrated in Figure 3.7 for the case of A 



We also make 




oo 



Figure 3.7. The components of the contour E in the case A — ■ Compare with Fig- 



ure 



3.2 



Since in this case it is fJj and not that meet the positive real axis near w — 1, 
the contour component E^q does not exist. It is generally the case that, aside from the 
two disconnected components of E>o which are always taken to be oriented to the right, all 
contour arcs are oriented with regions labelled "+ " on the left and regions labelled "— " on 
the right. 



the accumulation set of Pjv, ^oo '■= [n, &] U , into an oriented contour by using homotopy to E^ U E"^ 
fixing the points w — a, w — tn , and = b to assign orientation. Note that Poo, taken without regard to 
orientation, is the inverse image under A — E[w) of the interval < — 4«A < — G(0). The contour E^ U E"^ 
followed by —Poo forms the boundary of a bounded region that we denote by Z C C. We write Z = Z+ U Z_ , 
where according to the orientation of E^ U E^ U (—Poo), dZj^ is positively oriented while dZ- is negatively 
oriented. 



3.3. Analysis of the jump conditions for M(ti;). The jump conditions satisfied by the continuous bound- 
ary values of M(i(;) across the contour E will involve the product T\n{w), and to formulate the jump condi- 
tions in a concise way it will be convenient to consider the asymptotic behavior of ^{^{w) as N ^ oo. But 
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first, we introduce some elementary logarithm-type functions. For y ^ M+, let 



l{w-y) ■.= \og[ ^^^^' ] (3.7) 



denote the principal branch, obtained by composing the principal branches of the logarithm and the square 
roots. For fixed y, this is a function of w that is single-valued and analytic except along a piecewise-linear 
branch cut consisting of a finite line segment from w — y to w — Q along which we have the logarithmic jump 
condition 1+ (w;y) — I- {w;y) = 2iTi and the semi- infinite ray along which we have l+{w;y) + l-(w\y) =0. 
As w — > oo in the domain of analyticity, l{w] y) — > 0. For y < 0, we define a new function by setting 

{— ZTT, |w| < 1 and > 0, 

in, \w\ < 1 and < 0, (3.8) 

0, \w\>l. 

If y < —1, then m^(z/;; y) has a unique analytic continuation to the interval w G (—1, 0) and has branch cuts 
consisting of the intervals [y, —1] and M+ as well as the upper and lower arcs of the unit circle. If — 1 < y < 0, 
then m^{w; y) has a unique analytic continuation to the interval w G (y, 0) and has branch cuts consisting of 
the intervals [— l,y] and M+ as well as the upper and lower arcs of the unit circle. We also define a function 
m^'^{w;y) for y < by a completely analogous formula in which the conditions < 1 and \w\ > 1 are 
respectively replaced by the conditions that w lie inside and outside of the region bounded by the nonreal 
arcs of U T,^. Finally, if y = — e*'^ with — tt < a; < tt, then we set 

{—iiT, \w\ < 1 and — tt < arg(— w) < oj, 
iiT, \w\ < 1 and w < arg(— w) < tt, (3.9) 

0, \w\>l. 

This function has a unique analytic continuation to the line segment connecting w = y with w = 0, and its 
branch cuts consist of the unit circle and the interval . 
Now, set 

L%{w) := ^ m^{w;y)eN - ^ m^{w;y)eN + ^ ■m^{w;y)eN - ^ m^{w;y)eN- (3.10) 

yev yeA yeV yeA 

Q{y}=0 S{y}=0 GfalT^O sfyj^o 

(For the various choices of A under consideration, only one of the latter two sums will be present in any one 
case.) This function is analytic for w S C\ (Poo UE+) because each summand is. It also satisfies the identity 
L'pf{w*) — L^(w)*; indeed this holds term-by-term in the first two sums, and in each of the last two sums 
(whichever one is present) we may group the indices y in complex-conjugate pairs and the desired Schwartz 
symmetry holds for each pair. The boundary values taken by L^(ti;) on M4. satisfy L^^(w) -I- L'^_{'w) ~ 
as this identity holds term- by-term in each of the sums. Finally, we have the identity 

Hat (lu) = e^" ii;eC\(PooUM+). (3.11) 

This identity would be completely obvious were it not for the ztiir contributions in the definitions of m^{w; y) 
and mP{w;y); that it holds in the presence of these contributions follows from the fact that = V U A 
consists of an even number of points. 

Along any arc of P^, the part of Poo that is the accumulation set of V, we may enumerate the points of 
V in order {. . . , yn-i, ym Vn+i, ■ ■ • } consistent with the local orientation of Poo. Note that this order is the 



same as that of increasing k in the Bohr- Sommcrf eld quantization rule (1.42 1. Thus, from (1.42 1 we have 

6'o(y„+i) - 6'o(y„-i) 2TTeN- (3-12) 
Expanding the left-hand side around y„ gives the spacing between consecutive points of V as 

Ay(y„) ^ ^ ^ ^ ^ ^ 
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Therefore, for each w € C \ (Poo U 



Hy}=o ?s{y}^o 



0'o{y)m''{w;y)Ay{y) + - ^ 9'o{y)m^{w;y) Ay{y) 

yeV ySV 

{y}=o a{y}5^o (3.14) 

B'^{y)m^{w, y)dy+- f 0'„{y)m^iw; y) dy + 0{el), 



where the second-order accuracy comes from the fact that the sum is a midpoint rule approximation to the 
integral away from y = ~1 (the extreme sample points are asymptotically half as far from the endpoints of 



as they are from their neighbors due to the 1/2 in the Bohr-Sommerfeld quantization rule (1.42) and 



the choice (3.31 of the transition point w — tn) and the fact that E'{—1) = 0. In a similar way, but taking 



into account that = Prx, \ P^ is oriented oppositely to the direction of increasing k in (1.42) 



- ^ m^{w;y)eN- ^ m^{w;y)tN 

yeA y6A 
S{y}=0 s{y}#o 

= -/ e',{y)m^{w;y)dy+- [ 9'o{y)m^{w-y)dy + 0{ej,). (3.15) 

The error terms are in fact uniform in w for w bounded away from P^o, and in particular, for such w, 

L%{w) L°{w) + 0{elj) and UN{w)e-^°^"''>/''' ^ 1 + 0{eN), N ^ oo, (3.16) 



L\w)~-l 9'oiy)m''iw;y)dy+ - I e'^{y)m^[w-y) dy. (3.17) 



where 

L\w) :^ i 

This "continuum limit" of L^^w) is analytic for w G C \ (P00UM+) and shares the same symmetry properties 
as its discretization: L°{w*) — L^{w)*, and L° (w) + L^{w) = for it; e M+. 

Another useful form of L'^{w) is easily obtained by integrating by parts. In the integral over Poo H M, one 
expects endpoint contributions from the two points y — a and y = b, as well as from the points y — —1 and 
y — tn (the latter only if neither A nor V is empty so that a transition point exists) where the orientation 
of Poo changes. In the integral over Poo H (C \ M) one expects endpoint contributions from the two points 
at y = 1 on opposite sides of the branch cut and from the two points at y — —1. However, since 
m^{w] —1) = mP{w\ —1) the contributions from y = —1 will cancel between the two integrals. Also, in the 
integral over Poo n (C \ M) the sum of the contributions from the two endpoints at y = 1 vanishes because 
E{y) = for both points and m^{w; 1+) + m^{w] 1—) — 0. Furthermore, in the integral over Poo n M the 
contributions from the endpoints y = a and y — h both vanish individually because both correspond to 
E{y) = — iG(0)/4, the point at which 0o(y) vanishes. On the other hand, the contribution from y = tj\j 
generally survives, but it takes a particularly simple form: zL29Q{T]\[)m^{w;T]y)/Tr, where the "+" sign 
(respectively "— " sign) corresponds to the case when Poo is oriented toward (respectively away from) the 



transition point y = tn- But according to the condition (3.3) characterizing the transition point t^t, this 
contribution can be written as 2ejvnTO'^(w; tjv) where n G Z, and according to (3.4), 7i — (mod 2). 
Explicitly differentiating m^(ti;; y) and m^{w; y) with respect to y to integrate by parts then yields 

tQ, \ f ^o(y) dy N ^oiQ\ 

L iw) = / , \-2€Nnm (w^tn). (3.18) 

Jp^ y/-yy-w 

This form of L*'(w) allows us to exploit analyticity of 6q{-) on Poo to apply Cauchy's Theorem, showing that 
L'{w)^Liw) + r-"''^^ (mod2.ze^), w e _ 

w e C \ z, 




-w f 0oiy) dy ^ 



where 



L(w):=- / — — +2eNnm'^'^(w;TN). (3.20) 

In this definition of L{'w), the function 9o{y) in the integrand denotes the analytic continuation of the function 
of the same name from Poo- The domain of analyticity of is therefore C\ (S^ U UK+). Like L^{w), 
L{w) has Schwartz symmetry {L{w*) = L{w)*) and satisfies L+{w) + L-{w) = for w G M+. Furthermore, 



\0 (mod 47rzejv), fg(S^UE^)nr ^ ' 

and moreover the exact relation i+(0 ^ ^-{0 ~ 2i^?o(0 holds for those ^ e (E^ U E^) n M that do not lie 
on the branch cut of m^'^{w; tjv)- We define for future reference 

Y{w) ■.^UNiw)e-^^'"^/'"', (3.22) 

a quantity that by (3.16) and (3.201 is uniformly 1 + ©(e^y) for w bounded away from Z and Poo- 
Let L(^) denote the average of boundary values taken by L{w) on E^ U E'^: 

m-.^liL+iO + L-iO), eeS^UE^. (3.23) 
We now introduce two quantities closely related to this average: 

^^{0 :=i(0 + 

and 

^^iO :=I(e) + 







3{a > 0, 




0, 






(3.24) 


— iTrnejv, 




3{a < 0, 








> 0, 




0, 


^eE^, 


= 0, 


(3.25) 




^eE^, 


< 0. 





It follows from Proposition 1.2 that (p^ and are analytic functions on each arc of E^ and E^ respectively. 



The self-intersection point ^ = lies either in E^ or E^, and the additive constants in the definitions 
(3.24) and (3.251 ensure that for the relevant function analyticity extends to a full complex neighborhood of 
^ — ti;+. Now set 

T^iO ■■= 2n^(e) COS {e],'9o{0) ^ € E^ (3.26) 

and 

r^(e) := 211^(0-' cos {e],'9oiO) e^''^^^^- , ^ e E^. (3.27) 
According to the Bohr-Sommerfeld quantization rule ( |1.42[ ), these functions are analytic where defined (all 
singularities are removable), and we shall denote the analytic continuations from the contours of definition 
by the same symbols. Moreover, it can be proved that both {w) ~ 1 and T'^{w) ~ 1, a fact whose proof 
is only slightly more subtle than the analysis of Y{w) presented above (the proof given in jj] relates the 
product IIn(w) to the Gamma function and applies the reflection identity r(z)r(l — z) sin(7rz) — n and 
Stirling's formula). We formalize our results concerning the functions Y, T^, and T'^ with this proposition: 

Proposition 3.1 (Baik et. ai, 2 ). The function Y{w) is analytic for w £ C \ (E^ U E'^ U Poo U M+). 
Uniformly on compact sets in the domain of analyticity disjoint from the region Z in between P^o and 
E^UE^, 

Y{w) ^l + 0{eN). (3.28) 
This estimate fails if w approaches either a or b from within the above region of definition, but nonetheless 
both Y(w) and Y(w)~^ remain bounded if these points are approached nontangentially to the real axis. On 
the other hand, uniformly on compact subsets of Z± , 

Y{w) = e=F2»«o(»)/.iv (1 OieN)) ■ (3.29) 



The function T^ (w) is analytic for w G f^Y U flY \ (M+ U {a, b, t^}), and uniformly on each compact set 
in the interior of this region, 

T^{w) = l + OieN). (3.30) 
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This estimate fails if w approaches either a or b from within the above region of definition, but nonetheless 
{w) remains uniformly bounded near these points as ejv i 0. 

The function T'^{w) is analytic for w G D,^ U fJ"^ \ (M-i- U {a, b, tj\[}), and uniformly on each compact set 
in the interior of this region, 

T^iw)^l + OieN). (3.31) 
As above, this estimate fails if w approaches either a or b from within the region where T^{w) is defined but 
still T^{w) remains bounded. 

Finally, we have the algebraic relations 

rV(y^) ^y(y^) |^l_^g±2»eoWA„^ ^ (3.32) 

and 

T'^(u;) =y(u;)^i (l + e^^*^''^"')/'^") , w efl^\R. (3.33) 



The proof of (3.28) and (3.29) has essentially been given above, and the algebraic relations (3.32 1 and 
(3.33) follow directly from the definitions of Y[w), T^ [w), and T^{w) with the use of the jump condition 
(3.21). The asymptotic relations (3.30) and (3.31) are easily proved from the other results as long as w is 



kept bounded away from Poo- For example, suppose w € Q^, in which case T^{w) = Y{w)~^{l + e'^^^°^'"'>^'^") 
according to (3.33), and w lies on the right of so w ^ Z^. If also w G E^^{D^), then w lies on the right 
of Poo so w ^ in which case (3.28) applies and we obtain (3.31) from the inequality 5{6'o(u')} > 0. On 



the other hand if also w e E^^{D-), then w lies on the left of Poo so w e in which case (3.29) applies 
and we obtain (3.31) from the inequality 3{6'o(t«)} < 0. Note, however, that the case when w is near Poo 
must be considered separately; see [l] for details. 

With the notation of Y{w), T^{w), and T^{w) estabhshed, we may write down the jump condition 
satisfied by M(w) across the various arcs of S in a simple form. Indeed, we have 



M+(e) = M_(0 



-iT^(C)e[2»Q(?)+'^'(«)lAiv 



M+(0-M_(0 



M+(0=M_(0 
M+(e)=M_(e) 



ie-[2iQ(«)+i(?)-«eo(«)]/eiv' 



1 irA(^)g-[2iQ(C)+¥.^(5)]/ejv 

1 

iy(^)e[2»Q(?)+i(?)+ieo(«)]/£iv i 

I jy(^)-lg-[2iQ(?)+L(C)+ieo(«)]/eiv 

1 -|_ e-2ieo(?)/ejv 

1 0' 

ir(Oe[2''3(C)+-£'(«)±*eo(C)]/£« i 

1 -ir(^)-ie"l2*^(^)+-^(^)=F*eo(?)]/£jv 



VA 



_jy(^)e[2«Q(?)+-E'(«)-««o(?)]/ejv 

M+(0 = M_(0 



AV 



M+(0=M_(0 







1 



(3.34) 

(3.35) 
(3.36) 
(3.37) 
(3.38) 
(3.39) 



M+(0 =C72M_(0(T2 



M+(0 = a2M_(0CT2 



1 + g»[eo+(«)-e(,-(0]Aiv jY_(^)e[2^'5-(^)+-^-(«)-'^°~(^^]/'^"" 
-iy+(^)el2*'5+(^)+^+(?)+'^o+(^)]/'" 1 

eeE^o, (3-40) 

1 iy^(^)-ie-[2iQ+(?)+i+(«)+»«o+(?)]/ejv' 
_iy_(^)-le-[2iQ_(e)+L_(C)-ieo-(?)]/ejv i _,_ g-^l0o+ii)-Oo-m/^N 

^eE^o, (3.41) 

and finally, 

M+iO = a2M_(0^72, e e S>o. (3.42) 
In deriving the jump conditions (3.40) and ( |3.41 1 we used the fact that for ^ e M+, Q+(0 + Q-iO ^ 
and nAr_|_(^)njv_ (^) = 1. Written this way, the jump relations display the key importance of the exponents 
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2iQ{£) + L('C) ± i6'o(0; 2?(3(^) + and 2iQ{^) + ip^{£^). It remains to introduce a mechanism to control 

the corresponding exponentials, and this is the purpose of the next transformation. 



3.4. Control of exponentials. The so-called (7-function. Let g{w) be a scalar function analytic for 
w e C \ (S^ U U In terms of this to-be-determined function, a new unknown can be obtained from 
M(ti;) as follows: 

It follows that 

g-ie(?)/eN Q 
g-ie(?)/eN Q 



N+(0 = N_(0 
N+(0=N_(0 



(3.43) 
(3.44) 
(3.45) 



and 



N+(0=N_(0 



pie(«)/ejv 



N+(e) = N_(e) 



g-i9(5)/cjv jyA|-^-)g-[0(C)+i77e„#A]/e„ 







pie(«)/ejv 



^{a ^ 0, 



where 



(3.46) 
(3.47) 

9(0 := -lig+iO - g-iO) and 0(0 2zQ(C) 1(0 - 5+(e) - 5- (0, ^eS^US^. (3.48) 

These definitions assume that the indicated boundary values of g exist unambiguously (independent of 
direction of approach). The key to the Deift-Zhou steepest-descent method in this context is to choose 5(^1) 
and the non-real arcs of the contour or so that the functions 9{£^) and 0(0 have properties that can 
be exploited to make the matrix N(z/;) easy to approximate in the limit ejv — >■ 0. 
Consider the analytic matrix functions defined as follows: 



-0-3/2 



^[2iQ{w)+L{w}-ieoiw)-2g(w)]/tN 



L^(u;) := < 



I j^^-[2iQ(w)+L(w)+iea(w)-2g(w)]/tN 
1 

1 0' 

j^^[2iQ(w)+L(w)-ieo(w)-2g{w)]/tN ^ 



w e nf, 



(3.49) 



(3.50) 



1 

_^^[2iQ{w) + L(w)+ieo(w)-2g(w)]/eN I ' 

where the square roots T^(u')^/^ a nd T '^(w)^/^ are defined to be the principal branches, so that in view of 
( [330l ) and ( [33l| ) from Proposition [st] we have r^(w)^/2 i ^nd T'^lw)^/'^ « 1 throughout the domain of 
definition of [w) and L'^(w) respectively. 

Proposition 3.2. The jump conditions for N(i(;) on the contours S^'^ and E'^^ may be written in the 
form 

N+(OLy(o = N_(OL^(0, eeE^^ (3.51) 

and 

N+(OL^(0=N-(OL^(0, eeS^^. (3.52) 



Proof Note that on S^^ and E"^^ we hav e 



(0 - Q-(0, L+iO = L_(0, g+(0 = 5-(0, >+(0 = i^-(0, 



and 6*0+ (0 = ^o-(0- The formulae (3.511 and (3.52) then follow from taking principal branch square roots 
factor-by-factor in the algebraic identities ( |3.32 1 and (3.331, a meaningful step since T,^^ and S'^^ are 
disjoint from Z. □ 
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Proposition 3.3. Suppose that 



(3.53) 



Then, the jump conditions for N(w) on the contours E^q and E^q may be written in the form 



N+(0L^(e) = fT2N_(0LY(0^2 



5V('^)g[2jQ+(4)+L+(4)-2g+(4)]/£„ 1 + {^) 

^gE^o, (3-54) 



where 



\2 

A^{£,)^0{eN) and 5^ = O ( eat^c 



X^E+iO>0, 



(3.55) 



and 

N+(0L$(0 = a2N_(0L^(e)fT2 



l+yl^(^) 5A(-^^g-[2jQ+(C)+L+(e)-2g+(C)]/eN 
5A(^^^g[2jQ+(C)+L+(5)-2ff+(^)]/ejv ^ + {£,) 



where 



A^(0-O(ejv) «nd S^(e)-ofeAr^e-"^/^"V A 



= i?_(O>0. 



/n both (3.551 and (3.57) the parameter a > zs defined in Proposition 1.2 Also, 

N+(0 = a2N_(0^2, ^eS>o. 



^gE^o, (3.56) 



(3.57) 
(3.58) 



Proof. The relation (3.58) follows directly from (3.42) and (3.43) taking into account the given condition 



(3.53 ) on the boundary values oi g. To prove (3.55 ) and (3.57 ), note firstly that from ( f .39 ), (2.f ), and (3.20 ) 



one has both 



Q+(O + Q-(O=0 and e2 



[L+(?) + L_(C)]Aiv 



1, 



(3.59) 



From the latter it follows also that y+(^)y_(^) = 1 for ^ e 1R+. Since M.^ is disjoint from we may define 
as the p rinci pal branch of the square root for w just above and below and it follows from (3.28) 



from Proposition 3.1 that F+(^)^/^y_ (^)^/^ = 1 for ^ G M+ also. Now from taking principal branch square 



roots factor-by-factor in (3.32) and (3.33) one can express the factors T'^ {w)^/'^ and T^(i(;)^/^ appearing 
in ( |3.49[ ) and ( |3.50p in terms of Y{wy^ and the principal branches (1 + g±eo(u')/eN)i/2^ which are also 
well-defined for w G K+. 



One may now combine these facts and definitions with the jump conditions (3.40) and (3.41) to see that 

1 + e«[eo+(«)-eo-(0]Aiv 



(3.54) holds where for ^ G S>o, 



A^(0 := 
5^(0 



{I ^ g2ieo+(C)/eiv)l/2(l g-2i0o-(e)/e]v)l/2 
g«fo+(C)/ejv _ g-ieo-(C)/eN 



- 1 



(3.60) 



'(1 + g2ieo+ (?)/£«) 1/2(1 -f e-2»eo_(4)/eN)l/2 ' 



and that (3.56) holds where for ^ e S^q, 

A^(0 : 
B^(e) : 



1 + e-*[^°+(«)-^°-(«)l/'^" 

(1 e"2ieo+(C)/ejv)l/2(X _|_ g2ieo- )l/2 



(3.61) 



(1 _^ g-2i0o+(«)/ejv)l/2(l ^_ g2jeo_(C)/ejv)l/2 " 

The desired estimates on A^{^), B^{^), A'^{^), and B'^{0 now follow from Proposition 1.2 in particular 
from the Taylor series of 5" (A) about A = 0. Indeed, given ^ e M+, i?+(^) + E-{£,) = 0, and moreover £^+(^) 
is real and positive for ^ G E^q and real and negative for ^ G E^q (the change of sign comes from the reversal 
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of orientation). Therefore, recalling (3.1), using (1.52 1 from Proposition 1.2 together with Assumption 1.5 
and writing 



we find 
and 



-2a\/e]s 



B^{0 = 2{-1) 



N 



^-aX/tN sin 7,1 



-2a\/ct 



Eat 



(3.62) 



0(ew), X^E+{0>0 (3.63) 



sin(e^V(A^)) 



I _|_ g-2aA/ejvg2ii/(A2)/e„ | 



O I — e 



, A-i?+(O>0, 



(3.64) 



therefore proving (3.551. In exactly the same way one obtains the estimates (3.57), thereby completing the 
proof. □ 



Note that in the special case in which ^'(A) is given by the formula (1.57), the error terms in (3.55) and 



(3.57) vanish identically, a fact that was exploited to simplify the analysis of the semiclassical limit of the 



focusing nonlinear Schrodinger equation with corresponding special initial data in |17j . 

4. Construction of g{w) 

We proceed under the assumption that the Riemann-Hilbert problem for N will reduce, under appropriate 
changes of variables, to a problem solved using a genus-1 Riemann surface [i.e., an elliptic curve), at least 
for small t and x bounded away from ixcrit- This assumption is consistent with the qualitative behavior 



observed in Figure 1.2 in which the solution appears to have one oscillatory phase for these values of x 
and t. In other parts of the space-time plane, the solution appears to have more than one oscillatory phase, 
suggesting that the corresponding Riemann-Hilbert problem will be solved using a Riemann surface of genus 
greater than one, which will require modifying the following calculations. 

4.1. Two types of "genus-1" ansatz for g{w). Let p and q be real parameters and consider the quadratic 
polynomial p, q)^ given by 

R{w-p,qf ■.= {w-pf -q. (4.1) 

When working in a region of the space-time plane where more than one nonlinear phase is expected, it is 
necessary to choose R{w)^ to be a higher-degree polynomial (with a corresponding increase in the number 
of parameters pj, qj specifying the roots of R{w)^). 

We must distinguish two cases, which we label as "L" and "R" as these will correspond ultimately to local 
asymptotics for the fluxon condensate U]y{x,t) in terms of periodic librational and rotational wavetrains, 
respectively. 

|T] This case is defined by the inequality q < 0. The quadratic i?(z«;p,q)^ has distinct roots forming 
a complex-conjugate pair w — p ± iy/~q. The roots are assumed to lie on the nonreal arcs of the 
contour U E'^ (or, rather, given p and q with q < the regions ilj and il^ are assumed to be 
positioned so that this holds). We define a subcontour (3 C T,^ U T,^ consisting of the closure of 
the arc of S connecting the two roots of i?(u;;p,q)^ via w — 1. Thus /3 is a simple contour (with 
no self- intersection points). Whenever we are in case L, we will assume that there is no transition 
point, i.e., either A = or V = 0. 
This case is defined by the inequalities q > and 



R 



a<p-^/q<p + ^/q<b. (4.2) 

The distinct real roots w — p± y/q therefore lie in the real interval [a, b] of U S'^ . We then assume 
further that the point w ~ u)+ where the nonreal arcs of US'^ meet the interval [a, b] lies between 
the two roots. We define a subcontour /3 C U T,^ as the union of the closure of the nonreal arcs 
of Ty U with the real interval [p — ^/q, p + y/q]. Thus /3 is a non-simple (self-intersecting) contour 
having a single self-intersection point (a simple crossing) at w = . Whenever we are in case R, 
we will assume that if there is a transition point w ^ tn, it lies in (3. 
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In both cases, we define 7 := (E'^ U E^) \ f3, and we assume the contours /3 and 7 inherit the orientation of 
S. See Figures 4.1 and 4.2 for illustrations of (3 and 7 in cases L and R respectively. 





Figure 4.2. The subcontours /? (red) and 7 (blue) for a co nfigu ration of type R. The ori- 
entations here correspond to the case A = ^^v*^ (^^^ Figure 3.2). This illustrates the fact 
that if a transition point tn is present it is assumed to lie in the subcontour /3. 



We define an analytic branch R{w) = R{w; p, q) of the square root of the quadratic R{w; p, q)^ by taking 
the branch cut to coincide with (3 and choosing the sign so that R{w) = w + 0(1) as w 00. Note that in 
case R as well as the borderline case f = 0, R{w) is a sectionally analytic function of w because the branch 
cut locus (3 has a self- intersection point at w = 

It will be useful below to have available some compact notation for certain sums of contour integrals. We 
therefore define 

[ FiOd^:^ [ Fi^)d^+ [ FiOd^+ [ FiOd^+ [ FiOd^ (4.3) 

Jc Janj\sv JaoY\sv Jdnf\s^ Jafi^\s^ 
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and we will use this formula in some situations where F(^) is to be understood along the four contours on the 
right-hand side in the sense of taking a boundary value from within the region whose boundary is involved 
in the integration. For example, if the regions fij meet the positive real axis in the contour E^q ^-^^d if F 
is a function having a jump discontinuity across this contour, then the two terms on the right-hand side of 



(4.3) involving integration over dfl^ \ require use of two different boundary values -F±(C) taken by F on 
E^Q. If F{^) depends on a parameter w through a Cauchy factor {£, — w)^^ then F{^) will be analytic 
in w in a bounded domain we denote as ft" := {ft U E^ U E^) \ {o, Tjy, b, 1}. 

In order to construct suitable functions g{w) with which to transform the matrix M.{w) into N(u'), we will 
need to impose certain relations between the real parameters p and q and the independent variables x and 
t. We will define two functions M{p,q,x,t) and I{p,q,x,t). Then the "moment condition" M{p,q,x,t) = 
and the "integral condition" /(p, q, x,t) = will be used to study the dependence of the roots of R{w; p, q)^ 
on X and t. For fixed x and t, the solvability of the equations M — and / = for p and q will be 
a necessary condition for asymptotic reduction to a model Riemann-Hilbert problem solvable in terms of 
elliptic functions. However, this will not be a sufficient condition. It may happen that for some {x, t) this 
system is solvable but certain necessary inequalities (see Proposition 4.9) fail along arcs of the subcontour 
7. In this case, it becomes necessary to introduce new arcs of the subcontour (3 near the points where the 
inequalities have failed in 7. Thus, /? becomes disconnected, and there are more endpoints (i.e., roots of 
-R(u))^). To determine these additional parameters it then becomes necessary to also include further equations 
Mj — and Ij — Q. Ultimately this will effect an asymptotic reduction to a model Riemann-Hilbert problem 
solvable in terms of higher-genus hyperelliptic functions. 

For the moment, we think of (p,q,x,i)^ S as a parameter vector, and begin by defining a function 
M = M(p,q,x,t) as 

q Jc 
Now let H{w) = H{w; p, q, x, t) be the function defined by 



M 



i?(e;P,q) 



(4.4) 



Hlw) 



1 



t 



Tp^ ' n Jc R{^;p,q){i~w) 



(4.5) 



This function is analytic where it is defined (the singularity at w = and corresponding branch cut along M_|_ 
are excluded from Q°), but it does have jump discontinuities across the the contours E^q oi' ^>o (depending 
on which of the six cases of A we are considering), as well as across E^'^ and E^^ if there exists a transition 
point. Note that H{w) satisfies H{'w*)* = H{w), so that in particular the zeros of H{w) in its domain of 
definition either lie on the negative real axis or come in complex-conjugate pairs. We may now define a 
second function / = /(p, q, x, t) as 




R+{i)H{£,)d^ 




(4.6) 



R-iOHiO dt 



/3nC_ 



Finally, let /(w) 



f{w; p, q, X, t) be given by the formula 
dQ , _ j_\ , x~t R{w]p,q) _/?(?«; p,q) 



^^{w;x,t) 



8v^ 



%{OV^id^ , IdL 



R{^-p,q){^-w) 2dw 



iw). 



(4.7) 



Proposition 4.1. Let parameters p, q, x, andt be given so that the quadratic i?(w;p,q)^ is in case L, case 
R, or the borderline case 0/ q = 0, and assume that the regions fij and $7^ are chosen so that the contour 
E^ U E^ is consistent with the roots of the quadratic with well-defined subcontours (3 and 7. Suppose also 
that the moment condition M = and the integral condition I — Q both hold. Then an analytic function 



g(w) is defined for w € C\{/3U . 
[l\ In this case we set 



as follows. 



f{w') dw' 



(4.8) 
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where the path of integration is arbitrary m C \ (/? U . 
In this case we set 



g{w) < 



fiw') dw' , G To 



f{w')dw', we Too, 



(4.9) 



where Tq and Too o-fe respectively the bounded and unbounded connected components o/C\ (/3UM+), 
and in each case the path of integration is arbitrary in the given domain. 



(In the borderline case of q = we also use the definition (4.9).j The function g{w) so-defined satisfies in 
all cases the Schwartz-symmetry condition 

and is a uniformly Holder-^ continuous map C\{I3U R+) C. The function 6* : U E'^ 
(3.48) satisfies 



S{6»(0} = o, Ce(s^ui]'^)n: 



{0 = ^R+{t,P,c{)H{0, e€/3- 



Also, we have 

The function (j) : U S'^ 

for some real number 
and 



g+(0 + .g-(e) = 0, eeM+. 



defined by (3.48) satisfies 



Finally, g satisfies the decay condition 



''''(e) = i?(C;p,q)i?(0, eer 

lim g{w) = 0. 



d^ 



(4.10) 

C defined by 



4.11) 
4.12) 

4.13) 

4.14) 

4.15) 
4.16) 

4.17) 
4.18) 
4.19) 



Proof. Even without the conditions AI = and / = 0, it is obvious from the definition (4.7) that the function 
f{w) is analytic at least for w e C \ {Yj^ U S"^ U M+). From the Plemelj formula and the definition (3.20) of 
L{w), one finds that 

/+(e)-/-(0 = 0, eex (4.20) 
which, upon taking into account the continuity of the boundary values taken by / along 7, shows that / is 
analytic in the larger domain C \ (/3 U M_|_). 



Now f{w) is automatically integrable at it; = 0. Indeed, (3.20) shows that L(w) has a well-defined limiting 
value as w — > with | arg(— w)| < tt (and is in fact uniformly Holder-^ for such w), the term involving the 
integral over 7 is clearly 0{{-w)-^^^), and the sum of the remaining two terms is as well (after canceling a 
term proportional to {—w)~^^'^ between them). This fact gives sense to the formulae defining g{w) in case L 
and also in the domain Tq in case R. But the moment condition M = also makes / integrable at w = 00: 



x-t 4 
, +x+t+ - 



i?(C;P,q) 



+ 0((-Z«)-3/2) 



w — > 00. 



(4.21) 



By a simple contour deformation argument (using the fact that i?(^;p,q) changes sign across /?), we may 



write the integral over 7 as minus one-half of the corresponding integral over C (recall (4.3 )) and thus identify 
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the term in brackets as the "moment" M{p,q,x,t). Therefore in case R, g{w) is well-defined and analytic 
for w € Too, and in case L we observe that g has a well-defined limiting value as w — >■ oo. 



The Schwartz symmetry (4.10 ) of g now follows immediately from the definition of g and the corresponding 
symmetry /(z*) = f{z)* obvious from (4.7), and the Holder- ^ continuity of 5 can be read off from the formula 
for / = g' . Also, since by definition g has no jump across the contour 7 we see from the definition (3.481 of 
that ( |4.12| ) holds. 

Next, observe that 



/+(0-/-(0 = i?+(?;p,q) 



-C7^i?(e;P,q)(e-0 



(4.22) 



= -R+{OHiO, 

where the second line follows from a residue calculation and a contour deformation like that used in identifying 



the leading term in (4.21) with a multiple of M{p, q,x,t). Similarly, 

f+io-f-io = m)Hio, 



„ . dQ , ^ , dL 



Simpler calculations show that 



Together with the definitions (3.48) of 9 and 



integration of which yields (4.14), and 



, these two relations establish (4.13) and (4.18). 

/_(0 = 0, eeM+, 



d^ 



(4.23) 



(4.24) 



(4.25) 



The latter shows that (j) defined by (3.481 is constant in each arc of (3. Since it follows from the Schwartz 



symmetry (4.101 of g and the jump condition (4.14 1 that g takes purely imaginary boundary values along 
M+, an examination of the functions 9 and along /3 near its intersection point w — 1 with shows 
that the limiting values of these functions taken along /3 (either of the two arcs) as ^ € /3 tends to <^ = 1 
are, respectively, real and imaginary. Since <j>{^*) — 0(0* > the identities (4.16) follow immediately. Now, 



Schwartz symmetry of g (4.101 also shows that both 9{^) and 4>{0 are real-valued for ^ £ (E^ U S ) n 



proving (4.11) and (4.15). If we are in case R, so that /3 contains two real arcs, we can now show that = 



on these two arcs, proving ( |4.17 l. Indeed, for ^ e /? n M, we may use constancy of 4> along the arc to obtain 



JR{0(O} - Hm^ ^m')} 

{'G/3nR 

= ^{2iQ + L}{w+) 



hm JR{5+(0+ff-(n}, 



(4.26) 



Note that while 2iQ{^) + L{^) is not continuous on f3 in a neighborhood of the self-intersection point ^ = 



due to jump discontinuities in L (see ( |3.24 | and (3.25) and the discussion just below these definitions), its 
real part is, which makes ^{2iQ + L} well-defined at ^ = regardless of the arc along which this point is 
approached, and explains the notation di{2iQ + L}{w~^). But g+{£,') + g- (^') — 2g+{^') — i9{^') and 9{^') is 
real for ^' e /3 n M, so 



^m)} = + L}{w+) - 2 hm n{g+{a}, 



(4.27) 



Now if [/ is a small neighborhood of the self-intersection point w ~ of /3, then g : t/ \ /3 — > C is uniformly 
Holder- i continuous, so the latter limit may be taken along a different arc of /3 with the same result: 



^{(f>{^)} = ^{2iQ + L}{w^ 



2 hm ^{g+iO}, ee/3n] 
C'e/3nC± 
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(4.28) 



where the upper or lower half-plane is used depending on whether the original arc of /3 H M was oriented to 
the right or left respectively. Now we can write this in terms of the limiting values of (j) and 6 along f3 n C±: 



= lim n{i0{C)}, ^e/3n] 

«'e/3nc± 



(4.29) 



where the second line follows because is a purely imaginary constant ±i<i> along f3 D C±. Finally, since 
9{^) has a real limiting value as ^ — > 1 with ^ E /3 C±, we obtain that = for e /3 H E as a 

consequence of the integral condition / = 0. But from (4.15), 0(^) is real for ^ g M, so (4.17) follows. 



The decay condition (4.19) is obvious from the definition in case R, and in case L it follows from the 
integral condition / = 0. Indeed, writing g{oo) in case L as the the integral 



3(00) 



/(C) dC, 



(4.30) 



we may split the integral into two equal parts, and in each part we deform the contour into the right half- 
plane in opposite directions, bringing the two contours against /3UM+. Using (4.24) cancels the contributions 
to 5(00) coming from integrals along the upper and lower edges of ]R+, leaving only 

5(oo) = ±^(/+(e)-/-(0) (4.31) 

where the sign is different depending on whether /3 C or (3 C T,'^. But in either case, we now can use 
(4.22) and the symmetry /(C*) = /(C)* to see that the condition / = guarantees that 5(00) = in case 
L. □ 

If p and q can be eliminated by means of the equations M(p, q, x,t) = and /(p, q, x, t) — 0, then we may 
view g as being a function of w depending parametrically only on x and t. Now, according to Proposition 4.1 
there will be a real constant $ associated with g. In the situation that p and q have been eliminated, we will 
have $ = ^{x, t), and it will be useful to characterize the dependence of ^{x, t) on the remaining parameters 
X and t. 

Proposition 4.2. Suppose thatp = p{x,t) andq — q(a;,t) constitute a differentiable solution of the equations 
Af(p,q,x,t) = and I{p,c\,x,t) — 0, so that g{w,x,t) — g{w;p{x,t),c{{x,t),x,t), and let the constant 
$ = $(a;,i) be obtained therefrom as described in Proposition 4-1 Then ^(x,t) is jointly differentiable in x 
and t with real-valued first- order partial derivatives, depending on x and t only via p and q, given by 



dx 



TT 

41? 



1 - 



1 



where in case L, 



V 



and in case R, 



V = 



2K{mn) 



TT 

42? 



1 



1 



e (0,1), 



mp := 



4v/p2 - q 



G (0,1), 



(4.32) 



(4.33) 



(4.34) 



where K{-) denotes the complete elliptic integral of the first kind defined by ( |1.60[ ). In particular, defining a 
quantity rtp by 

„ _ 1 - Vp' - q 

and noting that p^ — q > m both cases L and R, we see that Up is an algebraic function of p and q satisfying 
\n-p\ < 1. Also, < p^ — q < 1 implies Hp > while p'^ — q > 1 implies Hp < 0. 
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i + Vp^ 



(4.35) 



Proof. Consider the functions X{w) and T(w) defined in terms of g{w) — g{'w;x,t) as follows: 

TT dg{w) , . TT dg(w) 



X{w) 



dx 



T{w) 



-wR[w) dt 



(4.36) 



These functions are both analytic for w e C \ /3. Indeed, from (4.14) in Proposition 4.1 we see that in spite 
of the explicit presence of the square root \/~w in the definitions, X{w) and T(w) may be considered to be 
analytic in a neighborhood of the positive real axis, and then analyticity for w G C\/3 follows from elementary 
properties of the remaining factors. It is also a consequence of (4.19) and the relation g'{w) = ,f{w) with 
/(w) given by Kjh that g{w) = 0{{-w)-^/^) 



T(w) = 0{w ^) as w oo. 



for large w, and this implies that X{w) = 0{w ^) and 
±i$ for ^ e /3 n C± and 0(0 = for 



Differentiation (with respect to x and t) of the identities 0(0 
^ S /3nM using the definition (3.48) of in terms of giw) yields the following jump conditions for X{w) and 
T{w) along /3: 



x+(0-x_(0 = 



T+(0-r-(0 = 



^i?+(0 



2*-B(0t« 

_2zi?(0, 

2zI?(0 T « 
2zi^(0, 



9$ 

9a: ' 



Ce/3n( 



(4.37) 



(Note that /3 n M = in case L.) 

Since we know directly from their definitions that the functions X{w) and T{w) must be 0{w^'^) for large 
w, they are necessarily given by Cauchy integrals via the Plemelj formula in terms of the jump data (4.37): 

m)d^ 

//3nC+ 

1 



X{w) = 



T{w) 



^R+iOi^-w) 
^p^R+iOiC-w) 



dx 
~dt 



^i?+(0(O 

d^ 



w) 



fSnC- 



d^ 
d^ 



w) 



(4.38) 



But while these formulae indeed exhibit decay for large w, without further conditions we will have only 
X{w) — 0{w~^) and T{w) — O^w'^) as w — > oo. Imposing on these formulae the more rapid required rate 
of decay of 0{w^'^) as w — > oo reveals conditions determining the partial derivatives of $ with respect to x 
and t: 



dx 



1 

V 



13 



where 



Now the fractions E{w)/y 



V : 

'/3nC+ 
and D{w)/^ 

E{w) _ i 
~ 4 



-w 



dS, and 



d^ 



1h 



1 

V 



^R+iO 



dO 



V-^R+iO 



d^ 



_ V-^R+iO' 



-w are simple rational functions of w. 



1 



and 



D{w) 



-w 



and by elementary contour deformations. 



l±w-i 1 
— —7 — ^ aw — - 
R{w) 2 



1± w- 
R{w) 



dw. 



(4.39) 



(4.40) 



(4.41) 



(4.42) 



where the first integral is over a small positively-oriented circle surrounding only w = and the second 
integral is over a large positively-oriented circle outside of which R is analytic. (This result holds regardless 
of whether the branch point configuration is of type L or of type R.) Evaluating these integrals by residues 
at w = and w = 00 respectively yields 



R+i.0 



RiO) 



(4.43) 
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Therefore, using i?(0) = — — q establishes the forniulae (4.32). 

It remains to characterize the denominator T) in cases L and R. In case L, elementary contour deformations 
show that 



1 



V = 



By the substitution w ^ s given by 



dw 
-z~ 1 



1 



dw 



v/-w((w-p)2 -q) 



z + 1 



followed by z = ±\/ 1 — 



(4.44) 



(4.45) 



and comparing with the definition ( 1.60 ) of the complete elliptic integral of the first kind, we establish ( 4.33| ). 
On the other hand, in case R, by simple contour deformations, 

_ dw f° dw 

V= / ^ ^ ^ = / ^ (4.46) 

Jo ^-wR{w) Ju+^ a/-w)((w) - p)2 - q) 

By the substitution w i-^ s given by 

2 



w 



2s 



a bijection mapping the real path from w 



Wl 



s = 0, we then obtain (4.34) by comparing with (1.60) 



p + y/q < to w = onto the real path from s 



(4.47) 

= 1 to 
□ 



4.2. Finding g when t = Q. 

Proposition 4.3. Suppose t — Q, and that A = for x>Q while V = for x <Q (for x — we may choose 
either A = or V = 0j. Then the equations Af (p, q, x,t) = and /(p, c\,x,t) =0 are satisfied identically if 



p{x) = l-^G{x) 



X e 



(4.48) 



and q — q{x) — p(a;)^ 



1. 



Note that the proof will show that the equation A/(p, p^ — 1, a;, 0) = has no real solution p if a; < and 
A = 0. This explains the need of introducing in general the set A C Pn and shows that scattering theory 
"from the right" is insufficient to capture the semiclassical asymptotics for all {x,t). 



Proof. Consistently with the condition q = p^ — 1 we choose the branch cuts /3 of R{w;p,p^ — 1) to lie on 
the unit circle (and some intervals of the negative real axis in case R). Then it follows that 

R{w; p, p2 - 1) = V^RiEiw); p), (4.49) 

where for p < 1, R{X;p)^ = —16X^ — 2(1 — p) with R{\,p) being analytic away from the branch cut in the 
A-plane lying along the imaginary axis between the two imaginary roots of -R(A; p)^, and with branch chosen 
so that R{X,p) = AiX + 0{X~^) as A — ?► oo. Therefore, with q = p'^ — 1 and t = the moment condition 
M — can be written in the form 

1 f 0'jndf 

x = - . "^^^ ^ . 4.50) 



Recalling (3.1), this suggests taking A = E{^) as the variable of integration, yielding 

x='-f (4.51) 
TTjEic) R{X;p) 

where the factor of 2 appears because the contours of integration making up the scheme denoted C contain 
pairs of distinct points (^, 1/^) that have the same image in the A-plane. Here E{C) denotes the images 
of these contours in the A-plane counted with the same multiplicities as in the definition of C; this simply 



means that in (4.51 1 

• if A = we are integrating over two contours from A = — «G'(0)/4 on the positive imaginary axis to 
A = 0, with the two contours lying on opposite sides of the branch cut of _R(A;p), while 
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• if V = we are integrating over the same two contours as in the case that A = but with the 
orientation reversed. 

By deforming these contours to the imaginary axis taking into account the change of sign of i?(A;p) across 
its branch cut, and using the substitution A = m)/4 we therefore obtain 

' ^ ' ip{v) — *(zi;/4), (4.52) 



^ V^^aZp) ^«2_2(l-p)' ■ dv 

where cr = 1 for A = and cr = — 1 for V = 0. It then foUows from Proposition 1 1 . 1| that 

X ^ aG-\-y/2{l-p)) (4.53) 

which can be solved for p — p{x) if either a; = or a; has the same sign as cr, yielding in all of these cases 
the formula ( |4.48[ ). 

If the condition M(p,q,a::, 0) = is used to eliminate the parameter x from H{w;p,q,x,0), then subject 
also to the condition q = — 1 we obtain 

Riw;p,p'-l)Hiw;p,p'-l,x,0) = -^RiEiwy,p) [ f ^ + 1-) ^ ■ (4.54) 

where p — p{x) is given by (4.48). But since C is mapped onto itself, with orientation preserved, by the 
involution ^ i— > 1/^, we have 

^-wR{E{0;p) Jc^-'~wR{E{0;py 
so averaging these two formulae we obtain for p — p{x) given by ( 4.48[ ) that 

Ri.;py - l)i/(.;p,p^ - 1,.,0) = lRiEH;P)EH'^ L^^^W^Wm^Y ^'-''^ 
With the substitution A = E{^) this becomes 

Riw;p,p'-l)Hiw;p,p'-l,x,0) = -RiEiw);p)Eiw)^ [ ] ..^^l^^'^^ - (4.57) 

dw Je{c) E[wy - R{X;p) 

Note that the integral factor is real if E{w) is imaginary. Now, in the definition (4.6) of /(p,q,a;,i), the 
integration variable ^ lies on the unit circle in the upper half-plane, and therefore E{^) is imaginary, as is 
E'iOdS, and the boundary value taken by R{E{S^);p) on the branch cut. It then follows immediately that 
the integral condition 7 = holds. □ 

Note that the relation q = p^ — 1 is suggested by the fact that when t — the inverse-scattering problem 
reduces to that for the Zakharov-Shabat system (1.40) under the mapping A = E{w). The condition 
q = p^ — 1 maps the radical i?(w;p,q) into the radical i?(A;p) whose branching points are known (for 
Klaus-Shaw potentials, see [18]) to lie on the imaginary axis in the A-plane. 

Proposition 4.4. Suppose that t = 0, and that A = for a; > while V = for x < (for x = we may 
choose either A = or \/ — Assume also that q = q(a;) = p(a;)^ — 1 where p = p{x) is given by (4.48), 
and that the contour /3 U 7 coincides with the union of the unit circle \w\ — 1 and the real interval [a, b]. 
Then for each x G M., an analytic function g : C \ {(3 U M+) — >■ C is well-defined by Proposition 4-1 with 
associated functions : /3 U j ^ C and : /3 U 7 — s- C defined by (3.48 1, and the following hold: 

• $ = 0. 

• < for € "f if A = and > for ^ j if V — fb. Moreover, 4>{^) is bounded away from 
zero for ^ £ 7 except in a neighborhood of either of the two roots of i?(^; p, q)^ (which are endpoints 
ofl)- 

• 0(^) and 9oiS,) — ^(^) dre both real and nondecreasing (nonincreasing) with orientation if A = $ (if 
V = 0/ Moreover, for ^ e /3, 9'{^) is bounded away from zero except in neighborhoods of the two 
roots o/i?(^; p,q)^ (endpoints of j3) and, in case R, the point ^ = —1. 

• H{^) = H[S^;p{x), q(x),x,0) is bounded away from zero for ^ G /3 U 7 except in a neighborhood of 
^ = — 1 where H{^) has a simple zero. 
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Recalling the definition ( 1.59 ) of Xcrit, we note that ( 4.48 1 shows that for t = 0, g is in case R for |a;| < XcHt 
and in case L for |a;| > XcHt, while \x\ — Xcvit is the borderline case. 



Proof. According to Proposition 4.1 g{'w) — > as w — > cxi and the boundary values taken by g on (3 U 
are related as follows: 



0, 

2ig(C)+L(0Ti$, 



C e /3nM, 
e e /3nC±. 



(4.58) 



Writing g{w) = \/—wR{w)s{w), we see that s is a function analytic for w G C\/3 satisfying s{w) = o{w "^Z^) 



as w — )> cxD (and therefore since /? is bounded, s(w) — 0{w ^) as w 
boundary values of s on /3 are now determined from (4.58): 



s+(0-S-(0 = 



1 



It follows that s is necessarily given by a Cauchy integral: 



oo). Moreover, the differences of 

(4.59) 



C e /3nM, 
^ e /3nC±. 



s{w) 



1 

27h 



1 



-^R+iO^-w 



-^R+iOi~w 



(4.60) 



Now this formula provides apparent decay at the rate s{w) = 0{w ^) as w — > oo, so the true faster rate of 
decay s{'w) = ©(w"^) implies that 



1 

2W 



2iQ{o+m 



(4.61) 



where V is given by ( |4.33| ) in case L and by (4.34) in case R, and is obviously nonzero in both cases. But 
a t = and if p = p{x) and q = q(a;) while the contour /3 is as specified in the hypotheses, then it is easy 
to check that /3 is mapped onto itself preserving orientation by the involution ^ — while the integrand 
changes sign under this involution. This proves that $ = for all x at t — 0. 
Now Proposition |4. 1| also asserts that 

R{OH{0 = ^'iO - 2zQ'iO + L'iO - 2g'(0, ^ G 7, (4-62) 
where R{£,) — R{£,;p{x), q(x)). Holding ^ G 7 fixed, we differentiate with respect to x: 



lm)HiO-2.E'iO~2^^, 



(4.63) 



Recalling the function X{'w) defined in the proof of Proposition 4.2 by (4.36) and given in explicit form by 
(4.38), we have 

^Ri^)Hi0^2^E'iO-~ 



dx 



e e 7. 



But, d^/dx = 0, so we may evaluate X{£) in closed form using (4.38) and simple contour deformations: 

E{C) dQ 



-Ci?+(C)(C-0 WS.) 



(4.64) 
ns: 

(4.65) 



where we have used the fact that — q = 1 implies i?(0) — —1. Therefore, recalling the definition (1.39) of 
E(w) we find simply that 



d 1 d 



Using (4.49) and the identity 



dR 
5A 



(A;p(a^)) 



I6A OR 



then yields 



dx 



RiOHiO = - 



p'{x) dx 
8E{0E'{0 OR 



^ e 7. 



(A;p(a:)) 



p'(x) dx 
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(£;(0;P(x)), CG7- 



(4.66) 
(4.67) 
(4.68) 



Now with (T = 1 for A = and cr = — 1 for V = we integrate from xo(^) := crG ^(4ii?(^)) to x: 

RmiO = ~mOE'{0 r |^(i?(0;P(2/))-^, ^ex (4.69) 



where we observe that the contribution from the lower hmit of integration vanishes because by definition 
(see (4.5)) 7J(^; p(a;o), q(a;o), lo, 0) is finite and i?(^; p(a;o), q(a;o)) = 0. If Wk is either of the two roots of 

R{w\p{x),q{x)Y , then we have 4){wk) = ±i$ = 0, so 



m 



i?(C)ff(C)rfC = -8 



E{0 



A 



^\\-My))^d\ (4.70) 



iE(wk) JaG-i(4«A) dy p' [y) 

where we have made the substitution A — E{C). So, if A = (so a: > 0) then x > aG^^{AiX) so dy > and 
p'(y) > while dR{\;p{y)) / dy < and ArfA < 0, yielding 0(^) < for ^ G 7. On the other hand, if V = 
(so X < 0) then x < aG^^iAiX) so dy < and p'(y) < wh ile dR{X;p{y))/dy > and AdA < 0, yielding 
(j>{^) > for ^ e 7. In both cases we easily obtain from (4.69) the inequality 



|e|V2|i?(i5(e);p(^))| 



7^(^(0; P(2/))^ 

aG-^(ME(£,)) P (y) 



> 



^m)\\E'm 

|C|i/2sup,,R|G(y)G'(y)|' 



Proposition |4.1| asserts further that 
Differentiating with respect to a; for ^ e /3 fixed gives 



dx 



tR+iOHiO 



d^dx d^dx 



Substituting from the explicit formula (4.65) gives 

which can be equivalently written in the form 

-tR+{OH{0 = 



v/-ei?+(e)(x+(o+x_(o) 



- ■«^(^)^'(^)^^-(i.(e);p(.)), e./?, 



dx 



p'ix) 



dx 



^ e 7. 

(4.71) 
(4.72) 

(4.73) 
(4.74) 
(4.75) 



where _R+(i?(^); p(x)) refers to the boundary value from the left as the (vertical) branch cut is traversed by 
E{£^) when moves along an oriented arc of f3. Integrating from x' = o'G^^(4ii?(^)) to x' — x then gives 



df) 

--{0^^R+m{0^~8^E{0E\0 
d^ 



^^^-{E(O-My))-^, (4.76) 



LG-H4iEii)) dy p'{y) 

Of course d6/d^ = for ^ g 7. This shows that dd/dv is a monotone function of x, where v = —4iE{£^) is a 
real parameter for /3 U 7. The extreme value is attained at x = 0, at which point 7 vanishes and according 
to (1.45), for each £, E (3 the right-hand side of (4.76) becomes 0o(C). This proves the desired monotonicity 

of0o{o-m- 

Now 9{wk) ~ when Wk denotes either of the two roots of R{w] p{x), q(a;)), so 



^(0 



R+{C)H{C)dC = -^i 



A 



^^"-(A;P(y))-^dA, ^ g /3 . 



dy 



p'iy) 



(4.77) 



From this formula it follows that 0{^) > for ^ e /3 regardless of whether A = or V = 0. Also, from (4.76) 
we have the inequality 



8\E{0\\E'{0\ 
|e|V2|i?(i?(e);p(x))| 



aG-H4iEiO) 



dR+ I ^1^^ ^1 dy 
dy 



-{E{i)-p{y))-, 



P'(y) 



> 



^\EmE\i)\ 



|^|V2sup,,„|G(y)G'(y)r 



(4.78) 



Combining (4.711 and (4.781 shows that H is bounded away from zero along /3U7, uniformly with respect 
to a;, except possibly in neighborhoods of points ^ e /3U7 at which either i?(^) = (corresponding to ^ = 1) 
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or E'{^) = (corresponding to ^ = —1). While H{w) indeed has a simple zero aX w — —1 as is clear from 
(4.57), there is no zero at w = 1 as we will now show. Suppose that x € M is fixed, and ^6/3. Then by 

s that 

(4.79) 



collapsing the contour E{C) to the imaginary axis and extracting a residue at A = E{£), (4.571 shows that 

*'(A) d\ 



I 



A2 



i?(A;p) 



Taking the limit along /3 of ^ — > 1 (either from the upper or lower half-plane) shows that H has nonzero 
limiting values contributed by the first (residue) term in the above formula. 

This completes the proof of the assertion regarding the function H. The assertions regarding the functions 
(/) and e then follow from this, the formulae 0'(O = i?(C)i?(0 for ^ G 7 and 9'{C) = iR+i^)H{^) for ^ € (3, 
and the inequalities e{^) > for ^ e /3, (/)(0 < for ^ e 7 if A = 0, and (/)(^) > or ^ e 7 if V = 0. □ 

4.3. Continuation of g to nonzero t. We begin by establishing some differential identities. Here we are 
viewing M and / as functions of the roots of R{w;p, q) rather than as functions of p and q themselves. 

Proposition 4.5. Let Wk, fc = 0, 1, denote either of the two roots of the quadratic R{w;p, q)^ . Then 

DM 



dwh 



(4.80) 



owk 4 



R+iOd^ 



R-iOd^ 



(4.81) 



where in each case the partial derivative with respect to Wk is calculated holding the other root fixed, that is, 
by the chain rule in which u and v are expressed in terms of the two roots. 



Proof. Writing R{w)^ ~ {w — wo){w — wi) it follows easily that 



OR 
dwk 



1 Riw) 
2w — Wk 



(4.82) 



where both w and the other root are held fixed. The identity (4.80) follows by substituting -y/p^ — q — 
\/—wo^—wi into (4.4), differentiating under the integral sign using ( 4.82[ ) (the contours of C are independent 
of Wk), and comparing with the definition (4.5) of i?. In a similar way, one shows that 



- — R{OH[Q = -- H{wk), 

dWk 2 V-'f £.-Wk 



and using this in (4.6) proves (4.81). 

Proposition 4.6. Let wq and wi denote the two roots ofR{w)^. The Jacobian 

{wo, wi) :— det 



of the map {wq,wi) i— {M,I) is 

J{wo,wi) 



(4.83) 
□ 



-dM 


dM- 


dwo 


dwi 


dl 


dl 


-dwo 


dwi. 



(4.84) 



-V^/-wq^/—wiH{wq)H{wi){wi — wq) 



(4.85) 



where H is defined by (4.5) and T) is defined by (4.33) in case L and by (4.34) in case R. 



Proof. This is an immediate consequence of Proposition 4.5 the definition of R{w), and the formula (4.40) 
for T>, which applies in both cases L and R. □ 
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4.3.1. Continuation from x ^ {0,±a;crit}- 
Proposition 4.7. There exist disjoint open neighborhoods 

ff^ n M = ~a;crit) o.i^'d 



^+ n: 



and in the {x,t) -plane with 

(Xcrit,+00) 



and 



such that, with 



U 



n M = (-Xcrit, 0) and ^+ n M = (0, Xcrit), 
?P U G^ U , the following hold true. 



(4.86) 
(4.87) 



There are differentiable maps p : — )■ M and q : 



uniguely determined by the properties that 



p{x,0) = l-^G{x)' 



id q(a;,0) = p(.T,0)^ - 1, {x,0) G G, 



(4.88) 



M{p{x,t),c\{x,t),x,t) = I{p{x,t),q{x,t),x,t) = 0, {x,t) e 
where it is assumed that in the definition of M and /, A = for {x,t) G G^ U 



(4.89) 

?l while V = % for 



The guantity np{x,t) defined for (x , t) £ (? in terms ofp{x,t) andc{{x,t) by (4.351 satisfies np{x,0) = 
and 



and 



dup 

Ik 

drip 
dt 



(x, 0) < 0, (x, 0) e u 



(x, 0) > 0, {x, 0) e u 



R • 



(4.90) 



(4.91) 



• A connected contour /3 (J consisting of the union of (i) a Schwartz- symmetric closed curve passing 
through w = 1 and enclosing the origin with (ii) the interval [a, h], can be chosen so that for each 
{x,t) a 6 an analytic function g : C\(/3UM+) — ^ C is well-defined by Proposition 4-T with associated 
functions : U 7 — > C and <^ : /3 U 7 — > C defined by (3.481, and so that the following hold: 

— The function cj) satisfies K{(/)(.^)} < for ^ e ^ if {x,t) \J 0+ and 3?{?!)(^)} > for 
S, € "f if (xjt) € ff^ U ^1^. Moreover, 5R{(/i(f)} is bounded away from zero for ^ € 7 except in a 
neighborhood of either of the two roots of R{S,;p{x,t),q{x,t))^ (which are endpoints ofj). 

— The function 0(^) is real and nondecreasing (nonincreasing) with orientation for ^ € /3 if 
{x,t) e (if {x,t) e i^R U ^fj. Moreover, 9'{^) is bounded away from zero except 
in neighborhoods of the two roots of R(^;p(x,t),q{x,t))^ (endpoints of j3) and, in case R, a 
single point ^ = < that converges to ^ = ~l asi— >-0. 

— The function H(^) — H{(^; p{x, t), q{x, t), x, t) is bounded away from zero for ^ € /3 U 7 except in 
a neighborhood of ^ ^ where H{^) has a simple zero. 

The notation is meant to suggest ttie fact ttiat tiie configuration of tiie roots of R{w\ p{x, t),q{x, t)Y is of 
type R for {x,t) G ff^ U ff^ and is of type L for {x,t) € U . We have therefore defined the function 
np(x, t) appearing in the statements of Theorems 1 1 . 1 1 and 1 1 . 2 1 in terms of p{x,t) and q(a;,t) for (a;,t) e ^ by 
(4.35). We are now also in a position to define the accompanying function £{x,t) for (x,t) € G: 



£{x,t) 



p{x,t) 



{x,t) e 



q{x,t) 

We can now also identify the region S\_ involved in the statement of Theorem 



1.1 



(4.92) 

as the union ff^ [J . 



Proof. The existence of the maps p and q is a consequence of the Implicit Function Theorem. Indeed, 
since for x ^ the roots and wi of i?(w; p(a;, 0), q(x, 0))^ lie within the domain of analyticity of H{w), 
Proposition |4.5| shows that M and / are differentiable with respect to these roots. Furthermore, since under 
the additional hypothesis |a;| 7^ Xcrit the roots of R{w;p{x, 0), q(x, 0))^ are distinct and neither is equal to —1, 
it follows from Propositions 4.4 and 4.6 that the Jacobian (4.84) is nonzero when {x,0) G G. Therefore we 
may solve uniquely for wq and w\ in terms of (x,t) from the equations Af = / = 0, and since p (wo + wi)/2 
and q = (wq — wi)^/4 we also have p and q. 
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According to (4.35 1 in Proposition 4.2 
Also, 

drip 



the fact that \){x, 0)^ — q{x, 0) = 1 implies that Up = when t ^ 0. 

1 an 



n 



WqWi 



(4.93) 



so that dup/dt and dll/dt have opposite signs. By differentiation of the equations M = / = with respect 
to t one obtains the system of equations 

dM dwo dM dwi dM ^ , dl dwo dl dwi dl 



dwo dt dwi dt 
from which follows the identity 



and 



an 



1 



dwo dt dwi dt 
dl 



dt 



0, 



(4.94) 



I ^ dM dM \ dl 

dt /{wq.wi) 

Now by noting the explicit t dependence in M and (via the definition of H) in / we have 



wo^ wi- — 

OWq OWi 



dl \ dM 
dt~ 



dM 
~dt~ 



— — and 



di 



4vn 



^+(C;p,q) 



(4.95) 



(4.96) 



Therefore, as i = im plies that n = 1, we also have dM/dt = when t — Q, simplifying the identity (4.95 1. 
Using also Proposition 4.5 and the fact that when t = 0, R+{(,; p, q) = ^/—^R-^-{E{^);p{x, 0)) gives 



dn 

~dt 



1 



[wiy/—WiH{wi) — Wo\/—'WoH{wo)] 3? ■ 



/3nC+ ^ 



(4.97) 



Using (4.691 with <t = sgn(a::) and taking the limit as ^ approaches either root Wk of R{w;p,q)^ from 7, we 



have the formula 



^/^^H{wk) = -4- 



-U{x), t = 0, 



where we have also used the identity D{w) — 2wE'{w), and where U{x) is defined as 



U[x) 



E{wk) lim 

567 



dy 



R{E{i)-p{xM 7.g„(.)G-M4.B(?)) R{E[0-MyM 



t = 0. 



(4.98) 



(4.99) 



The quantity U{x) has the same value regardless of whether Wk — wq or Wk — wi because E{wo) = E^wi) 
at t = 0, which explains why we omit any notational dependence of U on k. Finally, using (4.98) together 
with formula (4.85) from Proposition 4.6 we write ( |4.97 ) in the form 



on 

~dt 



D{wi) - D{wo) 



8VD{wo)D{wi)U{x) wi - Wo 



3? ■ 



i?+(£;(O;p(x,0)) 

/3nc+ s 



(4.100) 



Here we have used the relationship n — wqWi = 1 which is valid at < = 0. 
We now determine the phases of the various factors in this formula. 

• Since E{wg) ~ E(wi) is a positive imaginary number, and since G~^{A:iE{S^)) — ;> as ^ — > with 
^ e 7, while^ e 7 implies that G'^^(4ii?(^)) < and in the range of integration E(i?(^); p(y, 0)) < 
for ^ € 7, it follows that U{x) is imaginary and has the same sign as x. 

• According to (4.331 and (4.34) from Proposition 4.2 I? is a positive quantity. 

• By explicit calculation, 

D{wi) - D{wo) i 1 1 \ i 1 



Wi - Wo 



4 y/-Wo + ^-Wi 



1 



2 ^/-wq + V-wi 



t = 0, 



a quantity that is negative imaginary. 
A similar direct calculation shows that 

D{wo)D{wi) -^(1 + wo){l + wi) 



1 

16 



{l + wo){l + wi), t = 0. 



i6Vn 

and this quantity is positive real for (x, 0) G but is negative real for {x,0) € ff^ 
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(4.101) 



(4.102) 



• Since /3 n C+ is, for t = 0, an arc of the unit circle that we may take (without loss of generality) to be 
oriented in the counterclockwise direction, we see that d^/S, = id9, a positive imaginary increment, 
while the boundary value R+{E{^); u{x, 0)) is negative imaginary, and hence 

i?+(ii;(a;p(x,0))^ eM+. (4.103) 

;3nc+ ^ 

Combining these phases then yields the sign structure of dli/dt at i = that produces the desired sign 
structure for dup/dt at t = 0. 

We now describe how to construct the contours /3 and 7 to guarantee all of the corresponding conditions 
in the statement of the proposition. Of course all of these conditions are generalizations for t oi 
corresponding conditions that hold true when t = according to Proposition |4.4| when /3 U 7 is taken to 
coincide with the union of the unit circle with the interval [a, b] so our argument will be a perturbative one, 
in which the unit circle is replaced by a suitable nearby curve. Firstly, since when t — 0, the function H 
is bounded away from zero on /3 U 7 except at w = — 1 where H has a simple root, the same holds (also at 
t = 0) throughout £7° if the latter is chosen without loss of generality to be close enough to PUj. Now since 
H{w;p{x,t), q{x,t), x,t) is an analytic function of w depending continuously on {x,t) near (x, 0) and that 
satisfies H{w*) = H{w) it will also be bounded away from zero in il° for t small except near some real point 
w = close to ui = — 1 where it has a simple zero. It is easy to check that H{w;p{x,t), q{x,t), x,t) has 
(two different) analytic continuations to a neighborhood of w = 1 from 51° from the upper and lower half 
planes, so the limiting values H{l±;p(x,t), q{x,t), x,t) are both finite and nonzero. The contour /? H C+ is 
then obtained by solving the well-posed autonomous initial-value problem 

^^^iR+iC,p{x,t),q{x,t))Hii;p{x,t),qix,t),x,t), r > 0, ^(0)^1 (4.104) 

where we interpret H(l) as iJ(l+). Clearly, r parameterizes a trajectory along which "^{0} is constant and 
^{9} is nonincreasing with parametrization r, since 

^^ir^ - f ■ f - [^R+mm HR+mm* = - \^R+mm' < o. (4.105) 



The vector field of (4.104) varies continuously with time t, and the only critical points in fl° are ^ — and 
the two distinct roots of R^. Since an integral curve of this vector field for t ^ Q connected = 1 with the 
root of R^ in C_|. (in case L) or with ^ = — 1 (in case R) and since the Melnikov-type integral condition 7 = 



continues to hold true for t ^ 0, the solution of the initial- value problem (4.104) is a curve terminating at 
the perturbed root of R (in case L, in finite t) or at the point ^ = < (in case R, in infinite r). This arc 
together with its Schwartz reflection in C_ and, in case R, the real interval connecting the roots of R^, is f3 
for t ^ 0. Upon assigning its orientation according to whether case A = or V = holds, we easily see 



that since no critical points of the vector field of (4.104) lie in /3, the desired strict inequality for d9/d^ holds 



along (3. To construct the contour 7 cruder methods suffice. In case L, 7 is the union of the real interval 
[o, b] with a Schwartz-symmetrical arc connecting the two complex-conjugate roots of R^; we define 7 H C+ 
for t ^ as the image of the same for t = under the linear mapping 

wo{x,t)-w+ wo{x,t)+w+wo{x,0) 1nR^ 
wo(x,0) + l ^+ w,{x,0) + l ^^■^^^> 

{wQ{x,t) is the root of R^ in C+) taking ^ = ^0(2;, 0) to WQ{x,t) and ^ = —1 to w+. In case R there is 
nothing to do since 7 has to be the union of real intervals [a, Wq{x, t)] U [wi{x, t), b] where ^0(2;, t) < Wi{x, t) 
are the two roots of R? . In both cases, an easy continuity argument together with the fact that 3?{(/)} = 
at the simple roots of R^ and a local analysis of d(j>/d£^ = i?(^)i?(^) near these roots shows that the desired 
inequality for 5ft{(/)(<^)} holds on 7. □ 



It should be noted that elements of this proof actually provide computationally feasible numerical methods 
for continuation of 5 as x and t vary. 
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4.3.2. Continuation from x = 0. When x — t = 0, the roots of i?(w;p,q)^ coincide with a and b and this 
imphes that M and / are not difFerentiable with respect to the roots at this point. To circumvent this 
difhculty, it now becomes necessary to exploit more than the two simplest configurations of A and V first 
introduced in §3.1[ we must now consider the possibility that neither A nor V is empty. 

We begin with several observations concerning the functions M(p, q, x, t) and /(p, q, x, t). Fix some small 
5 > and consider configurations of type R in which the roots — P ~ Vq and Wy. — P + >/q of R{w; p, q)^ 
(in this section this will be more suggestive notation than wq and wi) satisfy a < < a + 6 < —1 < 
b — S < < b, which bounds q > away from zero. If we further fix two real values and with 
a + S < < — 1 < < b — (5, then we may compare M and / for the various cases listed in 53.1 in which 
we use the transition point Tqo — when we have A = or V = P^^, and we use the transition point 

Too 



' OO 



Pj^^ or V 



The first observation is that in this situation the functions (Af , /) are the same in the case A — P^ as 
in the case V — P^^ , and are the same in the case V = P^^ as in the case A = P^^ ■ To see this, it is 
useful to note that by a simple contour deformation in which the components of \ T,^ and dV,^ \ S"^ 
are collapsed toward /3 U 7 we may write M, originally defined by (4.4 1, in the form 



M = 



t 



+ x + t + 



4 fd'oi^EM 

^(C;p,q) 



(4.107) 



Vp^-^ 

and we see that the only way that this formula depends on the choice of A is via the orientation of the 
contour arcs in 7, and so the desired equivalence for M follows because the transition points lie in the 
complementary contour /3. Exactly the same contour deformations, when applied to the definition of H^w), 
will result in the additional contribution of a residue at ^ = if w e /3 we have: 



H{w) 



1 



t 



'Vp' 



R{e,p,^M-w 



iR+{w;p,q) 



(4.108) 



Our analysis of the formula (4.107) applies to all but the last term. This last term does indeed distinguish 
between A = P^^ and V = P^^ and between A = P^^ and V = P^^ due to a change of orientation of 
/3nC+, which changes the sign of the boundary value i?+(w;p,q). However, this discrepancy contributes 
nothing to the integral /, since recalling the definition (4.6) we have 




fi+(C;p,q) 



iR+{£.;P,q) 



= ±2^{0o{w+) - *(0)} = 0, 



(4.109) 



and for the remaining terms we note that integrating over the (oriented) contour /?nC+ against the (oriented) 
boundary value i?+(^;p,q) is an orientation-invariant operation. 

The second observation is that since q is bounded away from zero, we may consider M and /, in any 
of the six choices of A listed in 5 3.1 as well-defined functions of the roots < w^. These are both 



analytic functions of and Wy in the intervals o < < a + 6 and b — d < Wy < b. Next, recall the 
part of Proposition 1.2 guaranteeing that 9o{w) has an analytic continuation from w > a and w < b to 
small neighborhoods of it; = a and of w = b respectively. Using this fact, we may construct the analytic 
continuation A4^{{AI, I)} of {M, I) with respect to (holding Wy, x, and t fixed) about a small positively- 
oriented loop about w = a beginning and ending on the real axis with a < < a-\-d. We may also construct 
a similar analytic continuation denoted Aiy{{M, I)} with respect to w-^ (holding w^, x, and t fixed) about 
a small positively-oriented loop about w — b beginning and ending on the real axis with b — S < Wy < b. 
The third and most important observation is that these two monodromy operations and My simply 



result in involutive permutations among the various cases of choice of A enumerated in 5 3.1 These relations 



are shown in Figure 4.3 These facts are elementary consequences of the formulae (4.107) and (4.108). Indeed, 



the analytic continuation operations both leave the non- integral terms invariant, and Ai^ (respectively Aiy) 
may be realized in these formulae by analytic continuation of the integrand from the real axis and the 
generalization of the real segment of 7 near w = a (respectively it; = b) to a straight-line contour connecting 
w = a with w — w-^ (respectively connecting w — b with w — Wy). This latter operation results in a 
change of sign of the square root in the integrand, which is equivalent to the reversal of orientation of the 
corresponding segment of 7, and hence in the desired permutation of formulae. 
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A = 




V = 



Figure 4.3. The effect of the two monodromy generators on the four distinct types of 
function pairs (M, /) . 



Let be the two-sheeted Riemann surface obtained from taking two copies of the disc |w — a| < 5 
sht along w < a and identified in the usual way. Similarly, let V-^ be the two-sheeted Riemann surface 
obtained from taking two copies of the disc |w — b| < 5 slit along w > b and identified in the usual way. The 
multivalued square roots := {w — o)^/^ and ky := (b — wY^^ are global analytic coordinates for and 



F^ respectively. We now define a function F^ x F^^ x 



by setting 



{M,I)ik^,ky,x,t) := <^ 



(M, I){p,q,x,t)\^^^ , 
(M, I){p,q,x,t)\^^p-<K , 
(M, I){p,q,x,t)\^^p^K , 
{M, I){p,q,x,t)\^^^, 



> 0, n{ky} > 

^{k^} < 0, ^{ky} > 

^{k^} > 0, ^{ky} < 

3?{fc^} < 0, ^{ky} < 0, 



(4.110) 



where on the right-hand side, p = {w^ + Wy)/2, q — {w^ — Wy)'^ /A, and := a + fc^ while Wy b — fc^. 
In the second and third lines we could have equivalently used V = and A = respectively. The 
monodromy arguments above show that (M, I) is a pair of single- valued analytic functions on F^ x F;^ for 
each {x,t) e M^, except possibly on the coordinate axes fc^ = or fc^ = 0. But it is easy to see that the 
pair (M, /) extends continuously to the axes, and hence any singularities are removable, so [M, I) is a pair 
of analytic functions defined on the whole product F^ x F^. 

In fact, given the above discussion, it is not hard to write down explicit formulae for the functions M and 



/. Indeed, starting from the formula (4.107) for M and taking into account the change of orientation of the 
two arcs of 7 corresponding to changes in sign of k^ and ky , we arrive at the formula 



M{k^,ky,x,t) 



t 



y^l - akl + bk. 



kz^ky 



+ x + t 



4fc_ 



^0(0 + fc^s)\/— a — k'^sds 



kl 



ki 



4/c. 



(4.111) 



k% 



kls 



Since the last term of the rewritten formula (4.108) for H{w) contributes nothing to /, we may omit it and 
apply the same process as used to arrive at the above formula for M to obtain 



I{k^, ky, 



x,t) 




R+{w; p, q)H{w) dw 



(4.112) 



where p — ^{a + b + k"^ — k'^) and q = i(b — a — fc^ — fc^)^ and where 



H{w) 



X — t 



w 



akl 



bk\-k\kl 



6'Q(a + fc^s)-\/— — k\s ds 



Jo (a + fcH,S - W)\/1 - Sy^ 

Ak^ 



(4.113) 



6*^1(5 - kls)^J~b + klsds 



(b — fc^s — w)-\/l — s^/b 



kls 



In particular, it is easy to see that the two functions M and / are difFerentiable with respect to k^ 
and fc^ at the origin k^ ~ k^ — Q where, essentially, all six cases of (M, /) coincide with the same value, 
corresponding to the configuration = a and = b that occurs when x = t — Q. Therefore, we may 
compute the Jacobian of {M,I) with respect to the global analytic coordinates {k^,ky) of the manifold 



/(fc^,fc^) :=det 



4Pn 



^(0,0)^--^(G(0)^-4) 



) we define 


dM dM' 




dk^ dky 




di di 




dk^ dky . 




nan is 


-f-*(w/4) 
av 


v=-G{0) 



(4.114) 



(4.115) 



where Vq denotes the positive quantity V defined by (4.34 1 in the case that the roots p ± of R{w) 
taken to be a and b. Hence, ^(0,0) < via Proposition 
independent of x and t. 



1.2 



and Assumption 1.6 and ^(0,0) is also 



Proof. We compute the Jacobian ^(0, 0) by calculating the partial derivatives of ^ (k^, ky) for A;^ > and 
fc^ > and then letting k^ I and fc^ J, 0. In this situation we can use the formula {M,I) = (M, /)|^^ 
with the right-hand side evaluated at — a + k'^ and i 

■ dM dM 



ky . But then, by Proposition 



4.6 



we have 



/(fc^,fc^) =det 



dw^ 
dl 



dWy 

dl 



dw^ dw^. 



■ det 



A=0 



'dw^ 


dw^ ' 


dk^ 


dky 


dwy 


dwy 


Idk^ 


dky\ 



= 4P 



(4.116) 



fc2 



kl{b-a-ki 



where the notation i?A=0('i') specifies that the formula (4.5) is to be interpreted in the case A = 
the limit fc^ ^ and ky 10 and recalling that ab = 1 gives 

/(0,0) =4I?o(b-a)- lim k^H^^^ia + k^_,) ■ lim fc^^A=0(b - fc' )• 
By applying elementary contour deformations to (|4.5|) we see that for any sufficiently small d > 0, 



lim k^H^^ijj{a + k 



and 



1 lini r mV^^d^ 



1 k 

lim k^H\—n,{b — kt) = lim — , , , ^ , „ , . 

k^MO ^ ^-"^ TT fc^,fc^;0 ^-b + kl Jb+d RiOi^ - b + fc2 ) 
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Taking 
(4.117) 

(4.118) 
(4.119) 



Now, 



and 



a - d < e < a implies R{0{(. ~ a - k%) = {a + k% - O^^^^Jb - kl - £, 
b<^<b + d implies i?(0(? ~ b + k'l) = -{^ - b + klf^^^^~a~k^ 



(4.120) 



(4.121) 



where in both cases we mean the positive 3/2 power on the right-hand side. Therefore, by the substitution 
= a + kl,( and a dominated convergence argument. 



lim 



fc2) 



Similarly, but now using the substitution ^ = b — fc^C, 



lim — -= 



kl 



Vb^' 
Vb^' 



Therefore, (4.1171 becomes 



Since E{a) ^ E{b) = -iG(0)/4, 



/(0,0) 



^(0,0) 



16Vo 



■^'{~tG{0)/AfE'{a)E'{b). 



Next, since E'{a)E'{b) = -(a + b + 2)/64 and A = iv/A, 

d 



^(0,0) 



^2 



vI/(it;/4) 



v=-G{0) 



(o+b + 2). 



Finally, recalling the definitions (3.2) of a and 6, we obtain (4.115). 



(4.122) 



(4.123) 



(4.124) 



(4.125) 



(4.126) 



□ 



Proposition 4.9. There exists an open neighborhood ff^ of the origin (0,0) in the (x,t)-plane such that the 
following hold true. 

• There are differentiable maps k^ : 6^ — > M and ky. : M uniquely characterized by the properties 

that k^ (0, 0) = ky (0, 0) = and 

M{k^{x,t),ky{x,t),x,t) = i{k^{x,t),ky{x,t),x,t) = 0, (x,i) G G^. (4.127) 

Via the identification = a+fc^ andwy. = b — k^ and the relations p = {w^+Wy)/2 and q — {w^ — 
Wy.Y /A, we obtain a solution of the equations M{p{x, t), q(a;, t),x,t) = and I(p{x, t), q(a;, i), x, t) = 
corresponding to a configuration with 
— A = when k^ > and ky > 0, 



P^^ or V 



V = P, 



when fc_ 



N 



< and ky > 0, 
A = P^^ when k^ > and ky < 0, 



— V = % when k^ < and ky < 0. 
Also, the functions k^{x,t) and ky{x,t) satisfy 

dk^ 



and 



dt 

dk^ 
dx 



(0,0) <0 and 
(0,0) >0 and 



dky 
~dt 

dky 
dx 



(0,0) > 
(0,0) > 0. 



(4.128) 



(4.129) 



Let A and V be chosen according to the signs of k^(x,t) and ky{x,t) as above. Then, there is a 
Schwartz-. symmetric closed curve transversely intersecting the real axis (only) at w — 1 and w = 
(see below) such that with /3 chosen as the union of this curve with the interval [w^{x,t),Wy{x,t)] 
and 7 chosen as the union of closed intervals [a, w^(a;,t)] and [wy{x,t),b] (recall that the local 
orientation of these contours depends on choice of A), there is for each {x,t) e G^ an analytic 
function g : C \ (/3UM+) — > C well-defined by Proposition 4-1 with associated functions : /?U7 — > C 
and (/) : /3 U 7 -^■ C defined by (3.48), so that the following hold: 
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The function (p satisfies 4>{£,) < for ^ e jD and (/)(<^) > /or ^ € 7 n . Moreover, </){(,) is 
bounded away from zero for ^67 except in neighborhoods of W-j and Wy ( which are endpoints 
ofl)- 

The function 6{S,) is real and nondecreasing (nonincreasing) with orientation for ^ G /3 H E^ 
(for e /3nE'^j. Moreover, 6'{£^) is bounded away from zero except in neighborhoods of — 
and — Wy (endpoints of j3) and the simple root ^ = of H converging to — —1 as 
{x,t)^{0,0). 

The function H(^) = H{(^; p{x, t), c\{x, t), x, t) is bounded away from zero for ^ g /3 U 7 except in 
a neighborhood 0/^ = , a point converging to £, = —1 as {x,t) — )■ (0,0), where H(£^) has a 
simple zero. 



Proof. By the Implicit Function Theorem, it foUows from Proposition 4.8 that the equations M = and 
I — may be solved uniquely near {x,t) — (0,0) and {k^,ky) = (0,0) for differentiable functions k^{x,t) 
and ky{x,t), which are both easily seen to be real- valued for real {x,t) € ff^. Sign changes in these two 



functions correspond to sheet exchanges on the Riemann surfaces and F.^, so from the definition (4.110) 



we confirm that we are in fact solving A/ = / = in various cases of choice of A. The inequalities (4.128) 



follow from implicit differentiation of the equations M = and / = with respect to t at {x, t) = (0, 0) 



and (k^,ky) — (0,0), and arguments similar to those used in the proof of Proposition 4.8 to compute the 



limiting values of various partial derivatives. The inequalities (4.129) can be obtained similarly, but also 
may be understood in the context of Proposition |4.7| since increasing (decreasing) x with t = leads to an 
overlap of with (with and in the latter we have A = (V = 0) corresponding in the present 
context to k^ and ky both positive (both negative). 

The perturbation theory of the simple root ^ = of H near = — 1 and the proof of existence of the 
contour /3 H connecting ^ = 1 with ^ = w+ along which 0(^) is real and monotone both work exactly the 



same as in the proof of Proposition 4.7 although we should point out that given {x,t) e the contour 
/?nC+ will generally be a different curve, and w"^ a different negative real value, for different allowed choices 
of A (this situation is relevant if and only if k^{x,t)ky{x,t) < 0). The desired monotonicity of d{^) along 
the intervals (w^, w^) and {w^ , Wy) also follows by simple perturbation arguments from t = 0. 

Of course when k^(x,t)ky{x,t) — one or both intervals of 7 have collapsed to points, so it remains 
to show that cj), necessarily real for ^ e 7 C M_, actually satisfies the desired inequalities in 7 when the 
degenerate configuration 7 = {o, b} present at {x,t) — (0,0) is unfolded. This will follow from an analysis 
of H(w) valid when w is near either a or b and (x, t) is near (0, 0). For some small d > Q fixed, 

r tJv^^+^W' (4.130) 

ttV-w Ja^d R{^\ P, q)(C - w) 

where cr^ = 1 if A = 0, A = P^^ , or V P^^ and cr^ = -1 if V 0, V = P^^ , or A = P^^, and where 



1.2 



the error term 0(1) is uniform for {x,t) G ff^. By choosing d small enough, we have from Proposition 
and E[a) -iG(0)/4 that 

e'^{0>\W>^^ a-d<(<a, (4.131) 
Also, since ^ lies to the left of both and Wy where R < 0, 



^ >__^> ^ >^^>o (4.132) 



^(C;p,q) -R(C;P,q) Vw^ - a + d^AI^ 



holds in the same interval o — d < ^ < a. So from (4.130) we have 



-a,HH> ''f]f~'\'^f-^K oil)>0, u^ia. (4.133) 
27r(D - a + d)^J-w 

H{w) = r' Jif^iJ'^^ ^ + ^(1)' ^ t b, (4.134) 
TrV^wJfa i?(^;p,q)(C- w) 

where = 1 if V = 0, V = P^^, or A = P!^^ and cr^ = -1 if A = 0, A P^^, or V = P^^ , and where 



Similarly, 



^ , Ul ^ J JY ailU. Uy _LI1^ — (U, ^ — ,U1 V 

the error term 0(1) is uniform for (x, t) G Over the interval of integration we have from Proposition 
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1.2 



and the fact that E{b) = -iG(0)/4 that 

- W > -Id'oW > 0, b<e<b + d, 
and since now ^ hes to the right of both and Wy in a region where again R < 0, 



> 



also holds for b < ^ < b + d, so 



> 



y=b" 



> 



V-b-d 



> 



cjyH{w) > 



mV^b^l log(b-«;)| 



+ 0(1), wtb. 



(4.135) 



(4.136) 



(4.137) 



27r(b - a + d)y/^ 

The inequality ( |4.133 ) shows that a^H{w) is large and negative for w to the right of w = a, while (4.137) 
shows that a^Hiw) is large and negative for w to the left of w = b, with both statements holding uniformly 
for {x,t) G Therefore, since i?(^;p,q) < both for a < ^ < and for Wy. < ^ < b, we learn that 
0'(^) = H{^)R{£^\ p, q) has the same sign as ct^ for a < ^ < and has the same sign as cr^ for Wy < ^ < b. 



Since (/)(^) = both when ^ 



and when ^ 



we obtain the desired inequalities on (f) for ^ G 7. □ 



Note that via the construction of g for each (x,t) G ff^ and Proposition 4.1 Proposition 4.9 effectively 
guarantees the existence of a real constant $ = $(x,t) that might at first sight appear to depend upon 
various artificial details of the choice of A. However, according to Proposition |4.2[ the partial derivatives 
d^/dx and d^/dt are necessarily given in terms of the functions u : — K and v : M., and this 

makes $(x,t) well-defined for {x,t) g up to a constant. The constant is then determined by the fact 
that $(0,0) = as guaranteed by Proposition 4.4 by taking limits from nonzero a; at t = 0. Therefore, 



$ : - 
defined in terms of g in 
The functions u : ff^ 



is indeed a well-defined differentiable function that also agrees with corresponding functions 



7^ where these domains overlap. 



and V 



agree with those previously defined in the overlap region 



?R n (€?f^ U €^p, ), and this allows us to extend the definitions (4.35) of np{x,t) and (4.92| of £{x,t) in a 



U 



The region Sr involved in the formulation of Theorem 



1.2 



consistent way to the full union U 
is exactly this union with two curves t — t±{x) omitted. The curve t = t^{x) is obtained by solving the 
equation k^{x,t) = for t, which is possible near the origin according to the inequalities (4.128 )-( 4.129[ ), 
and the signs of < — t^{x) and k^{x,t) are opposite. Similarly the curve t = t-{x) is obtained from the 
equation ky{x,t) = 0, and the signs of i — t_ (x) and ky {x, t) coincide. Therefore along the curve t — t+ {x) , 
we have — P — y/^ = n, while along the curve t = t^{x), we have Wy = p + -/q = b. The inequafities 
( 1.67 )-( 1.68) are a consequence of the inequalities — a + kt. > a and Wy 



k^ < b that hold for 



(2:,^) e Sr because k^{x,t) and ky{x,t) are real. 



4.4. Solution of the Whitham modulation equations. We are now in a position to describe how the 
functions p = p(a;, t) and q = q(a::, t) relate to the formally-derived Whitham modulation equations described 
in the introduction. 

Proposition 4.10. Let p — p{x,t) and q = c\ix,t) be the unique solutions of the equations M = I = Q 



matching the given initial data as described in Proposition ^.1 (or Proposition 4-9 for {x,t) near (0, 0)J. 
Let Up = np(x,t) and £ = £(x,t) be calculated explicitly from p and q via (4.35) and (4.92| respectively. 



Then these fields satisfy the Whitham modulation equations ( 1.32[ ) for superluminal rotational waves (when 
p and q correspond to case R, that is, for (x,<) (z Sr) and for superluminal librational waves (when p and q 
correspond to case L, that is, for {x,t) G Si). 

The rest of this section is devoted to the proof of this proposition. In fact, the roots wq and wi of the 
quadratic R{w;p, q)^, when expressed in terms of rip and £, will turn out to be Riemann invariants for the 



system ( 1.32 ) 



4.4.1. Diagonalization of the Whitham system. We begin by writing the Whitham system (1.32) in terms 



of standard complete elliptic integrals by evaluating J{£) given by (1.20) in the two cases of modulated 
superluminal wavetrains of librational and rotational types. The condition of superluminality (w^ > fc^) 



implies that for librational wavetrains we have J{£) 

Ir{-£). 



L\_{~'£) and for rotational wavetrains we have J{£) 
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For superluminal rotational waves, we recall the definition ( 1.17[ ) and obtain: 



J{8) = /r(-£) 



1 



7r\/2 



v/cos((/)) + £( 



V2 



\/cos{(j)) +£( 



for £ > 1. By the substitution cos((/)) = 1 — 2z^, this becomes simply 

4 /I +£ 



J{8) 



-E 



for superluminal rotational waves, 



where E(m) is the standard complete elliptic integral of the second kind: 



E{m) := / 

"'0 



dz, < m < 1. 



For superluminal librational waves, we recall instead the definition (1.16) and obtain: 
V2 



J{£) = I^{-£) 



\/cos(0) + £ I 



2^/2 



\/cos(0) + £ ( 



for — 1 < £ < 1. By the substitution cos(0) ~ 1 — (1 + £)z'^ ^ this becomes simply 



--(l-£)if 



1 + f 



for superluminal librational waves. 



(4.138) 



(4.139) 



(4.140) 



(4.141) 



(4.142) 



where K{m) is the standard complete elliptic integral of the first kind defined for < m < 1 by (1.60). 
Now recall the differential identities jl] 



Ejm) - (1 - m)K{m) Ejm) - K{m) 

K m) — ; r , E (m) = . 

^ ' 2to(1 - m) ^ ' 2m 



(4.143) 



Since in each case J(£) is a linear combination of K(m) and E{m) and m is a function of £, it is clear that 
K{m) and E{m) can be eliminated between J{£), J'{£), and J"{£) to yield a linear second-order differential 
equation satisfied by J{£). We shall now show that this equation is the same in both cases. In the case of 
superluminal rotational waves, the above formula for E' {m) yields 



J'{£) = 



1 



TT V 1 +£ 



K 



l + £ 



for superluminal rotational waves, 



and then differentiating again, eliminating K and E, and comparing with (4.139) we obtain 

1 



J"{£) 



J{£). 



4(1 -£-2) 

Similarly, in the case of superluminal librational waves, we have 



J'{£) = -K 



l + £ 



for superluminal librational waves, 



(4.144) 



(4.145) 



(4.146) 



so after one further differentiation we eliminate K and E and compare with (4.142) to arrive again at exactly 
the same second-order linear differential equation (4.145). 

Given two complex-conjugate or real and negative variables wo and wi, we express £ and n-p in terms of 
these by the relations 

Wo+Wi 1 - JWiiWi 



£ = 



(4.147) 



2y/woWi ' 1 + ^wqWi 

It is obvious that if wq and wi are real and negative, then £ > 1, and this corresponds to the case of 
modulated superluminal rotational waves. It is also clear that if wq and wi are complex conjugates, then 
— 1 < £ < 1, and this corresponds to the case of modulated superluminal librational waves. The Jacobian of 
the map {wo,wi) H> {np,£) is 

Wi Wq 



dup/dwa dup/dwi 
d£ /dwo d£ I dwi 



(1 + ^WqWi) V'^Ofi'l (1 + y^WoWi)'^y/WoWi 

Wi — Wq Wq — Wi 
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(4.148) 



Therefore, if we write the Whitham system in the form 



d_ 
Ft 



£ 



A(np,£; 



dx 



£ 



0, 



(4.149) 



where A(np,f) is the coefRcient matrix appearing in equation (1.32) (and of course we can equally well 
express this explicitly in terms of wq and wi using (4.147)) then we may equivalently write this as a system 
of equations for and wi as 



Ft 



Wo 
Wi 



+ S-1AS 



dx 



Wo 
Wi 



0. 



(4.150) 



It is now a direct calculation to show that, as a consequence of the differential equation (4.1451, S AS 



is a diagonal matrix, that is, the columns of the Jacobian J are independent eigenvectors of A. Moreover, 
once J"{£) is eliminated using (4.1451, the denominator J\f{np,£) defined by (1.331 factors as a difference 



of squares, and the numerator of each of the diagonal elements of S ^AS is divisible by exactly one of the 
two factors in J\f{np,£). It follows that the Whitham system can be written in the diagonal form 



9wj 
dt 



J = 0, 1 



where the characteristic velocities (eigenvalues of A) are given by 

_ y/woWi{l + ^woWi)J{£) + {wj - wi-j)(l - y/woWi)J'{£) 
^ ' ^woWi{l ~ ^woWi)J{£) + {wj - t«i-j)(l + ^woWi)J'{£) ' 
In other words, the variables wo and wi are Riemann invariants for the Whitham system (1.32) 



(4.151) 



(4.152) 



4.4.2. Diagonal quasilinear system satisfied by the roots of R{w)^ . We will now show that if wo{x,t) and 
wi{x,t) are obtained by solving the moment and integral conditions M = I = 0, then they also satisfy a 
system of partial differential equations in Riemann-invariant form. Implicitly differentiating M and / with 
respect to x and t gives 

MwoWo,(x,t) + M^iWi^{x,t) + M(^^t) = and IwoWo^x^t) + IwtWi^(^,t) + = 0- (4.153) 

Solving for the partial derivatives of Wj with respect to x and t assuming that the Jacobian determinant of 
{M,I) with respect to {wo,wi) is nonzero, we can easily confirm the identities 



dwi 



dwi 



= 0, J =0,1 



where 



dt •' dx 



J = 0, 1. 



(4.154) 



(4.155) 



We will now show that while the various partial derivatives appearing in (4.1551 contain explicit depen- 
dence on X and t, as well as dependence on initial data through the function 0o{w), the combinations Cj are 



in fact functions of wo and wi alone. Indeed, recalling (4.96), we have 



y/wowi ~ 1 



and It — 



1 



(4.156) 



and by similar explicit calculations using the definitions of M and / (the latter in terms of the function _ff ) 
we obtain 



\/woWy 



1 



and Ix = ^ 



1 



;3nC- 



(4.157) 



These partial derivatives are obviously functions of wo and wi (independent oix and and also not depending 
on the initial data). The partial derivatives of M and / with respect to wi_j are obtained from equations 



(4.80) and (4.81) respectively in the statement of Proposition 4.5 and while these depend also on a;, t, and 
^0 via a common factor of H{wi-j), this factor will cancel out of the expression for Cj. The result is that 

. _ ^ - (1 - ^woWi)Bj 



(4.158) 
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where 



and 



A ■= 



/3nC+ W— 



/3nC_ 



d^ 



(4.159) 



(4.160) 



It is now obvious that Cj — Cj{wa, wi) is an exphcit function of wq and wi only with no dependence on initial 
data. 

We may put Cj into a form that admits a comparison with the characteristic velocities Cj of the Whitham 
system (1.32) by noting that regardless of whether the radical R is in case R or in case L and regardless of 

1 in the 

dC+ I ^^TT^T^T^d^ (4.161) 



whether in the former case there may exist transition points, we may integrate by parts to write A in the 
form 



A^ 



2g - Wq - Wi 

R-iOV^ 



/3nC+ R+iOV^ JpnC- 
Therefore, some elementary algebraic manipulations show that cj may be written in the form 



y/woWi{l - ^wqWi)U + {wj - wi-j)(l + ^wqWi)V 



where 



and 



U 



A 



2g - Wq ~ Wi 

V^R+iO 



dC 



2£,-wo- wi 



/3^C-^ 



<R+iO 



<R-iO' 



<R-iO 



d^ 



(4.162) 



(4.163) 



(4.164) 



To complete the proof of Proposition 



4.10 



it therefore remains to show that the functions Cj{wo, wi) coincide 
with the Whitham characteristic velocities Cj(wo, wi), a calculation that is different in cases R and L. 



4.4.3. Equivalence of cj and Cj in case R. Comparing (4.162) with (4.152) we observe that to prove that 
Cj = Cj holds in case R it is sufficient to show that there is some nonvanishing quantity C{wq,wi) that is 
symmetrical in its two variables such that the following identities hold: 

2 



and 



CV = J'{£) = - 



TT V 1 



l + E 
2 

l + £ 



where £ — £{wq,wi) is defined by (4.1471. 



First, we evaluate V in case R. By elementary contour deformations, we have 

di „ dw 



V = 2 







= 2 



tui ^J-w{w - wo)(w - Wi) ' 



where < wi < 0. By the substitution (4.47) we then obtain (after some nontrivial algebra) 

2 



V 



[woWxY/^ V l + £ 



l + £ 



Therefore, the identity (4.166) holds if 



27r 



(4.165) 
(4.166) 

(4.167) 
(4.168) 
(4.169) 
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With C so determined, we carry out the same contour deformations and make exactly the same substitution 
\AA7\ to obtain 



1 2w ~ wq ~ wi 



7r(w;oWi)^/*^ Jwi y/-w(w - wo){w - wi) 



dw (4.170) 



4 jl+E fWl - ms^ - (1 - s^) ds 2 , , 

■/ ' ™— e(0,l). 



In the final integral, the numerator of the integrand cannot be broken up without introducing a nonintegrable 
singularity at s = 0. Now let Ri{s)'^ = 1 — and R2{s)'^ = 1 — ms^, take the branch cut of i?i(s) to 
be the real interval [—1,1] with Ri ^ is a.s s oo, and take the branch cut of R2{s) to be the union 
(— cx), — l/v^m] U [1/^m, +oo) with i?2(0) = 1- Then by elementary contour deformations. 



^/Y^^ - (1 - S^) ds _ \ ^/l - ms2 _ (l _ ^ ^ y- ^ „ ^2 „ ^^(-g) 

where L is a positively-oriented loop surrounding the branch cut of -Ri(s) but lying in the domain of ana- 
lyticity of R2{s). In this formulation the integral can indeed be broken into two separate integrals as s = 
is no longer on the path of integration. Therefore, 

^/l-ms-^-{l-s-') ds ^ \ f Ri{s) ds If ds 
Jo VT^Vl - ms^ s2 4/ii?2(s)s2 4/isi?i(s)- ' 

Since sRi{s) is analytic and nonvanishing outside of L, and since sRi{s) ~ zs^ as s — ^ cx), the second integral 
vanishes identically. By another contour deformation and the substitution x — l/(sy/rn), we therefore obtain 



^ ^ - ' ; -^^ / rfa; ^ £;(m). (4.173) 



%/l - s2^1 _ ms2 s2 Vtos2 - 1 s2 



a; 



Using this result in (4.1701 confirms the identity (4.165) and completes the proof that Cj = Cj in case R. 

AAA. Equivalence of Cj and Cj in case L. To prove that Cj = Cj in case L, it is sufficient to show that for 
some symmetrical nonvanishing function C{wo,wi) we have the identities 



CV = J'iS) = Ik (^-±^) . (4.175) 



and 

nv = .T'fF] = _ , 

TT V 2 

To prove these identities we consider, for p = 0, 1, the integral 

^ ^ f {2^ - Wq - Wi)P d^ r {2^ ~ wq - w,)P d^ {2C-wo^Wird^ 

where the second equality follows from a simple contour deformation, and i is a positively-oriented loop 
surrounding the branch cut of R, which we recall for case L is an arc connecting wq := re*^ and wi — re^*^ 
with r > and < 9 < it, passing through the real axis only at ^ = 1. We make the substitution 

C = (4.177) 

which maps the points ^ = re^*^ to z = ±itan(6'/2) and maps ^ = to z = — 1 and ^ = cx) to z = 1. Taking 
care with the sign of i?(^(z)) we find that Ip becomes 

J ^ 1 / 2r y [(l-cos(g)) + (l + cos(0))zFdz 

" zy^cos((?/2) i_,t,„(e/2) VI VT^^z^ +ta.nie/2)^ 
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In the case p = 0, set z = ix ta.n{d /2) to obtain 



dx 



^/f cos(6l/2) 7_i \/r^^^l + tan(6l/2)2a;2 Vrcos{0/2)Jo y/l^ y^TTtMdJW^ 



dx 



Then, set t = \/l — x'^ to obtain 



if 



1 + f 



To evaluate Ii , first muhiply the numerator and denominator of tlie integrand by 1 + z to find 



icos(6l/2) 
zcos(0/2) 



+»ta„(e/2) (J _ + 2z + (I - cos(6i))z2 

-»ta„(e/2) (I - z2)yr^^v/^2 _^ tan(6'/2)2 
+»ta„(e/2) „ ^ _ cos(6i))z2 

'itan(e/2 



(4.179) 



(4.180) 



(4.181) 



(1 - z2)^/^^^^z2 + tan(6l/2)2 

where we have used even/odd parity to simplify the result. Since (1 — cos(0)) + (1 
cos(6'))(l - z^), we get 



dz 



cos(6'))z2 



4VF 



+itan(0/2) 



(iz 



itan(0/2) 
+itan(e/2] 



icos{9/2) J_ 

icos{6/2) Jo 
Again using the substitution z 



/i = 8^^cos(6l/2)2 f 
Jo 



(1 - z2)3/2^22 +tan(6i/2)2 
dz 

{1- z^ f/^y/z^ +t&n{9/2f 
- ixta.n{9/2) followed by t = \/l 
dt 



- 2r(l + cos(6l))/o 
-A{wow,Y/Hl-£)I<(^- 



(4.182) 



a; 2 we arrive at 



\/r^(l-mt2)3/ 



--4(u;o«^i)^/^(l-f)i^(^) 



1+f 



(4.183) 



Finally, making the monotone decreasing substitution f i— > s given by (1 — m<2)(l — ms^) — 1 — m, we obtain 



(4.184) 



Since U = h/y/wowi and V = Iq, the formulae ( |4.180| and ( [4.184^ confirm the identities ( |4.174[ ) and ( |4.175[ ) 
upon making the choice of C{wq,wi) = {wqWi)^/'^ /tt . This completes the proof that Cj — cj in case L. 

5. Use of g{w) 



5.1. Opening a lens. Recall the matrix unknown N(z/;) defined by (3.43), and the matrix functions (w) 
and L^(w) defined by (3.49) and (3.50) respectively. The most crucial step in the Deift-Zhou steepest 
descent method is the opening of a "lens" about the contour /3, which in the current context takes the form 
of an explicit substitution leading to a new unknown 0{w) as follows. Let the lens A be an open subset 
of Q° containing /? as indicated in various configurations (depending on case L or case R and in the latter 
case whether or where a transition point exists) in Figures 5.1-5.4 The precise shape of the lens is not 
important, but rather that it lies close to /3 without coinciding except at the endpoints, that it fully abuts 
the real segments S^q or S>oi ^^'^ that if a transition point Tjy is present in case R then A fully abuts the 
contours E^^^ and E^^. 
Now set 

'n(u;)l^(w), wen'^nA 

N{w)L^{w), wen'^nA (5.1) 
^N(u;), w eC\(AUdQUR+). 

The reason for including the conto urs E ^'^ and E^^ in the interior of A should a transition point exist is 
now clear: according to Proposition 



0{w) 



3.2 



0{w) can be analytically continued to the contours E^^ and E^^, 
so 0{w) no longer has any jump discontinuity across any contour emanating into the complex plane from a 
transition point. The contour of discontinuity for 0(w) is therefore Eq ■— (E \ (E^^ U E'^^)) U dA, as is 
illustrated with solid curves in Figures |5.1H5.4[ 
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Figure 5.1. The contour Eq of discontinuity of the sectionally analytic function 0{w) in 
case L assuming either A = or V = 0. The lens A is shaded. 




Figure 5.2. The contour Eq of discontinuity of the sectionally analytic function 0(w) 
case R with either A = or V = 0. The lens A is shaded. 



Recall that on the arcs /3 n C± we have 4>{£,) = ±j<i> whereas on /3 n M we have </>{£_) = 0. Set 

$A = $A(.T,t):=$(2:,t)+^eAr#A, (5.2) 



where recall # A is defined in ( 3.4 ) . Then by direct calculation using the definition ( 3.49 ) , the jump condition 
(3.44) for N(u;), the definitions (3.48) of 9 and (j) in terms of boundary values of g, and the jump condition 
(3.21) for L, we see that 0{w) satisfies the jump conditions 



o+(o = o_(e) 







_jgT**A/«JV 





, ^ e /3 n n ( 



and 



O+(0 = O_(0 
Similarly, from (3.50) and ( 3.46[ ) we obtain 

o+(o = o_(e) 





-i 



, c G /3 n E^ n : 



ie=F**A/ejv 

jg±**A/ejv 
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, ^ e /3 n E^ n C± 



(5.3) 
(5.4) 

(5.5) 



Figure 5.3. The contour Eq of discontinuity of the sectionally analytic function 0{w) in 
case R with either A — Pji^ or V = Pjv^ ■ The lens A is shaded, and the dashed curves 



emanating from the transition point tn £ (3 do not belong to So according to Proposition 3.2 



Figure 5.4. The contour Eq of discontinuity of the sectionally analytic function 0(w) in 
case R with either A = P^^ or V = P^^ . The lens A is shaded, and the dashed curves 
emanating from the transition point tj^ £ j3 do not belong to Eq according to Proposition \3.^ 



and 



o+(0 = o_(e) 



(5.6) 



Therefore, the jump matrix for 0(w) along /3 is piecewise constant with respect to ^ (but has nontrivial 
dependence on the parameters x and t via $ = $(x, t) S M). Moreover, we observe that the jump matrices 
in (5.3) and (5.4) are respectively inverse to those in (|5.5|) and (5.6), while the contours /3n E"^ and (30 E^ 



are sub-arcs of (3 that are oppositely oriented. This shows that (5.3|-(5.4) are describing the same jump 



conditions as are (5.5)-(5.6) 



Next, we consider the jump conditions relating the boundary values taken by 0{w) from the upper and 



lower half-planes along the positive real axis, recalling Proposition 3.3 
that the boundary values of N(w) and of 0{w) agree for ^ S E>o, we have 



Clearly, from (3.58) and the fact 



O+(0 =a2O_(0fT2, eGE>o. 
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(5.7) 



To analyze the jump conditions for 0{w) on the contours E^g or S^q (only one or the other of which is present 
in any of the six choices of A under consideration), we need to calculate the exponent 2zQ+(^)+L+(^)— 2g_|_(^) 
appearing in ( 3.54[ ) and (3.56). Of course this exponent has an analytic continuation into il^ or as 
2iQ{w) + L{'w) — 2g{w), and taking the boundary value on /? near = 1 from the relevant domain we learn 
that 2i(5+(0 + -^+(C) ~ is the analytic continuation through the adjacent domain fi^ or fl^ of the 

function + i9o{^) — i0{S,) + zTre^r^A (mod 27riejv) defined on /3. Since (/>{£_) = ±z$ e iR along /?, we 
then have 



±[2iQ+(C) + L+(?)-2g+(C)]/cjv 



and 



±[2iQ+ (0-23+ 



\2 

O { CN—e 



O { eN — e 



-aA/e„g|sR{je(?)-jeo(5)}|/£jv 



A = £;+(0 >o, 



(5.8) 



A = £;_(C)>0, 



A 

>0' 



(5.9) 

where in b oth c ases the aforementioned analytic continuation of 0(^) — 0{){£,) from j3 is implied. Due to 
Proposition 1.2 this difference has an analytic continuation from /3 to a neighborhood of ^ = 1 (note that 
this will be a different analytic function depending upon which of the two arcs of /? meeting at ^ = 1 is 
involved), and furthermore, from the integral condition / = 0, we will have ^{i9{\) — i0o(l)} = 0. Using the 



formula (4.13), we can easily obtain that 



*^'(0-<(e) = 



^+(^;p,q) 



^i?(C;P,q)(C-0 



(5.10) 



so letting ^ 1 along /3 yields a purely imaginary quantity in the limit. It follows that di{i9{£^) — iOo{£^)} — 
O(A^) as A = -£'±(1^) — > for ^ real. Therefore, by choosing the width parameter Si of the rectangles D± 



defined after the statement of Proposition 1.2 to be sufficiently small we will have |3fi{j6'(^) — i6'o(C)}| ^ aA/2, 
and then from Proposition |3.3| we will have that 

e e E^n u E^ 



O+(e)=a2O_(0'T2 (l + 0(ejv)): 



Now suppose that ^ G 7. Then from (5.1) we have 0±(^) 



satisfied by 0{w ) we may recall the jump conditions for N(w) in the form ( 3.44 1-( 3.47) 
Proposition 3.1 the factors {£,) and T'^{£^) are uniformly bounded on 7 n E^ and 7 n E^ respectively. 



>o-->o- (5-11) 
N±(^), so to evaluate the jump conditions 

According to 



Furthermore, by Proposition |4.l[ we have 9{S,) = for ^ € 7. Finally, Proposition 4.7 or Proposition 4.9 
guarantees that 3?{0(^)} is strictly negative for ^£70 E^ and strictly positive for 1^ G 7 n E^ as long as is 
bounded away from the band endpoints. We conclude that 0+(^) = 0_(^)(I + exponentially small in eAr) 
holds uniformly for ^ G 7 bounded away from both band endpoints. 

Consider next the jump of 0{w) across the boundary of the lens A. We assume that the arcs of dA inherit 
orientation from the arcs of the band /? that they enclose. The matrix N(w) is continuous across 9A, so 
from (5.1 1, 

'o_(OLV(0^\ ^edAnnl 



o+(0 



(5.12) 



Referring to the definitions ( 3.49 )-( 3.50), we see that these may be written in the form 

0(e-[2iQ(?)+i(«)T«eo(5)-25(|)]/e„-)~j ^ 



O+(0 = 0_(e) (l + O(T^(0'/' - 1) + 0(T^(C)"'/' - 1) 

C £ dAnnl, 

and 

O+(0 = O_(0 (l + O(r^(0'/' - 1) + O{T^{0'^/^ - 1) + O(e[2'Q(«)+^«)T»''°(«)-29(0]A«)^ ^ 

^edAnn^, 

assuming that in each case the three error terms on the right-hand side are bounded. If as eat ^ the 
point ^ remains bounded away from the singular points a and b of T'^(-) and T'^{-), then according to 
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(5.13) 



(5.14) 



Proposition 3.1 the first two error terms in each case are ©(ejv)- The exponent 2iQ{^) + L{^)^ido{(_) — 2 g(^) 



4.1 



is, modulo iireN, the analytic continuation from /3 to SI J of (p{£,) T Since according to Proposition 

(/)(^) is an imaginary constant in (3, and since according to Proposition 4.7 or |4.9| ^(^) is analytic, real, and 
increasing along parts of f3 in with derivative bounded away from zero away from the band endpoints, it 
follows from the Cauchy-Riemann equations that ^{2iQ{(^) + i(^) =F *^o(0 ~ '^9(0} strictly positive for 
^ G (9Anf2j bounded away from the band endpoints. Therefore the final error term in (5.13) is exponentially 



small as Eat ^ uniformly for ^ bounded away from band endpoints. Completely analogous reasoning yields 
the same conclusion for the final error term in (5.14). We conclude that 0+(^) = 0_(^)(I + O(eAr)) holds 
for ^ e dA as long as ^ is bounded away from the two band endpoints (which, given the shape of the lens A, 
also implies that ^ is bounded away from a and b). 

Finally let us consider the jump conditions satisfied by 0{w) across the contours and S^, that 
is the contours that make up the boundary of the whole region Q. For ^ on any of these curves we have 
0± (^) = N± (^). Therefore, to calculate the jump matrices on these curves we need to combine the definition 
(3.43) of N(w) in terms of M(u') with the jump conditions ( 3.38 )-( 3.39). Since g{w) is analytic on and 



S^, we obtain the jump conditions 



o+(0 = o_(0 



1 



o+(o = o-(e) 



_jy(^)e[2«'3(«)+-f'(«)±*''«(«)-29(?)]/eiv 

_iy(-^)-lg-[2iQ(e)+L(C)Tieo(?)-25(4)]/e„ 
1 



(5.15) 



(5.16) 



The exponent 2iQ{^) + L{^) ± i9o{£,) — 2g{£^) occurring in (5.15) is, modulo iireN, the analytic continuation 
into the domain fij from the contour of the function 0(^) =F id{£^) ± 2i6'o(C). If i = 0, then (j>{^) and Oo{^) 
are real, and < 0, for ^ e S'^. According to Proposition 



Also, by Proposition 



1.2 



4.4 



)(C) ^ ^(C) is real and nondecreasing for 



Oq{^) is strictly increasing with derivative bounded below by a positive 
constant for ^ € S^, and hence the same is true of 26*0(5) ~ ^(0- It follows by a Cauchy-Riemann argument 
that if the width 2Si of the rectangle -D+ U D_ that is the closure of the image of SI under E{-) is sufficiently 
small, then when t = 0, di{2iQ{(,) + L{^) ± j6'o(C) — '^9(0} be strictly negative on all parts of with 
the possible exception of points near a and b where the contours meet E^. But if the band endpoints 



are bounded away from a and b, that is (according to Proposition 4.3) if x is bounded away from zero, then 
4>{a) and 0(b) will be strictly negative by Proposition 4.4 so the strict inequality on the exponent in (5.15) 



(b) will be strictly negative by Proposition 
persists right down to the real axis when t = 0. Since for x away from x = the inequality is strict uniformly 
on E^, it also persists uniformly on Ej for sufficiently small t ^ 0. (Assuming \t\ sufficiently small also 
ensures that the nonreal arcs of E^ or E^ are contained within SI.) Completely analogous arguments also 
show that the exponential factor occurring in (5.16) is exponentially small uniformly on E^ if x is bounded 



away from zero and t is small enough. If x approaches zero, then the exponential decay only fails near £, = a 
and 5 = b. Now, the factors Y{^) and Y{£_)~^ are uniformly bounded on Ej and E^ respectively, according 
to Proposition 3.1 We conclude that for t sufficiently small, 0+(5) = 0_(5)(I + exponentially small) holds 
uniformly for S E J and ^ G E^ except near o and b when x is also small. 

We formalize these observations concerning the jump conditions satisfied by O^w) in the following propo- 
sition. 



Proposition 5.1. Suppose that the point {x,t) lies in one of the domains G'^ , (see Proposition ^.l), 



or ff^ (see Proposition 4-9), and \t\ is sufficiently small. Let Uq and Ui be discs of small radius independent 
of X, t, and ejv centered at the band endpoints WQ(x,t) and Wi{x,t) (see Figures 5.8- 5.11). Then for ^ (z /3, 
0{w) satisfies exactly the piecewise constant jump conditions (5.3)-(5.6). For £ M-)., 0{w) satisfies 
exactly the .simple jump condition (5.7) except in a .small interval near ^ — 1 where the ©(eat) approximate 
version ( 5.11[ ) of this condition holds. Finally uniformly for ^ G Eq \ (/3 U U C/q U C/i), we have simply 
O+(0 = O_(e)(I + O(ejv)). 

5.2. Construction of a global parametrix. A global parametrix for 0{w) is a sectionally-analytic matrix 
function 0{w) designed to satisfy the jump conditions of 0{w) for ^ e Eq n (/3 U E+ U UqU Ui), with the 
only modification being that we use the jump condition (5.7) on all of rather than omitting a small 
interval near ^ = 1 where 0{w) satisfies the approximate relation (5.11). The global parametrix 0{w) will 
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be analytic (have no jump discontinuity) for ^ G So \ (/3 U M+ U C/q U Ui), that is, on the arcs of the jump 



contour Eq where Proposition 5.1 guarantees that 04.(^) — 0_(^)(I + 0{eN))- It is standard to construct 
0{w) by patching together (i) an outer parametrix denoted 0°"*(w) expected to be a valid approximation 
of 0{w) away from the two roots of i?(w;p,q)^ and (ii) two inner parametrices denoted □"(w) and 0']"(w) 
expected to be valid approximations of 0{w) for w Cz Uq and w ^ Ui respectively: 



OH (5.17) 
0°"*(w), otherwise. 



It is also an important part of the construction of 0{w) that the outer and inner parametrices are well- 
matched on the boundaries dUk of the discs. 

Our construction of 0{w) in this section rests essentially upon the following two facts that hold true for 
{x,t) e SlU Sr: 

• The roots of R{w; p, q)^ are distinct. Error terms will become uncontrolled if the roots are allowed 
to coalesce, which occurs only on the common boundary of Si and Sr consisting of the two points 
X = ibxcrit- In a forthcoming paper [5j we will give a complete analysis of the interesting dynamics 
that occur for x ~ ±Xcrit and |t| small. 

• No root of R{w;p,q)^ coincides with a or b. Again, error estimates will fail if one of the roots 
approaches either a or b, as occurs along the two curves t — t±(x) excluded from Sr along which the 
choice of A must be changed as explained in Proposition |4.9| Unlike the coalescence of two roots 
of p,q)^, which as mentioned above leads to interesting new dynamics, the exclusion of the 
curves t — t±{x) appears to us to be serving a merely technical purpose, allowing us to avoid more 
complicated inner parametrices. As mentioned in §1.3| there appears to be no exceptional behavior 
near the omitted curves in the (x, t)-plane. 

The outer parametrix must be considered separately in the two cases (L and R) because the contour 
geometry is qualitatively different. Taking into account just the jump conditions for 0{w) on the band /3 
and the positive real axis K+, the outer parametrix 0°"*(?ii) is to be determined as a solution of one or the 
other of the following two Riemann-Hilbert problems. The outer parametrix will depend parametrically on 
the fast phase ly := ^a/^n as well as on the geometry of the contour /3, although both of these dependencies 
are suppressed in our notation. 

Riemann-Hilbert Problem 5.1 (Outer parametrix, case L). Let v he a real number, and let wq E C+. 

Let (3^ denote an oriented arc in C+ from w — I to w — Wq and let /3_ denote the complex- conjugated arc 
from w = 1 to w ~ Wi = w^. Find a matrix 0°"'(w) with the following properties. 

Analyticity: 0°"*(?i;) is an analytic function of w for w £ C \ (/3+ U /?_ U M+) and Holder-^ 
continuous for any 7 < 1 with the exception of arbitrarily small neighborhoods of the points w = wq 
and w — wi — Wq. In the neighborhood Uk ofwk, the elements o/0°"*(?x;) are bounded by a multiple 
of \w - Wk\^^/'^. 

Jump condition: The boundary values taken by 0°"*(k;) along f3± and M.^ (the latter oriented 
from the origin to +00 ) satisfy the following jump conditions: 

6°"'(0 = O°"*(0«CTie^'""', Ce/3±, (5.18) 

and 

6r(e) = ^26°"*(0a2, ^eR+. (5.19) 

Normalization: The following normalization condition holds: 

lim 6™*(u;) =1. (5.20) 



5.5 



The jump conditions satisfied by 0°"*(w) in case L are summarized in Figure 

Riemann-Hilbert Problem 5.2 (Outer parametrix, case R). Let v be a real number, and let wq < < 
Wi < be given. Let j3+ denote an oriented arc in C+ from w ~ 1 to w ~ let /3_ denote the complex- 
conjugated arc in C_ from w — 1 to w — w'^ , and let [3^ and j3-^ denote real arcs oriented from w — to 
w = Wq and w — Wi respectively. Find a matrix 0°^^{w) with the following properties. 
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-co 



(0 = 0°^\^)iaie-"'''-' 



Figure 5.5. The jump conditions satisfied by the matrix 0°^^{w) that we will show is 
normalized in the stronger sense that 6°"*(m;) -I + Odwl^i/^) 

as w 00. The only other 
singularities admitted are 0{\w — near the points Wq and Wq. 



Analyticity: 0°"*(w) is an analytic function of w for w C \ (/3+ U /3_ U U U Holder-^ 
continuous for any 7 < 1 with the exception of arbitrarily small neighborhoods of the points w — wq 
and w = Wi. In the neighborhood Uk ofwk, the elements o/0°"*(?ii) are bounded by a multiple of 

\W - lU/cl"^/''. 

Jump condition: The boundary values taken by 0°"*(w) along j3±, j3^, (3-^, and M_|_ (the latter 
oriented from the origin to +00 ) satisfy the following jump conditions: 



and 



Normalization: The following normalization condition holds: 

lim 6°"*(u;) =1. 



(5.21) 
(5.22) 

(5.23) 
(5.24) 



The jump conditions satisfied by O {w) in case R are summarized in Figure 5.6 



07* (C) = 0°J^\i)iaie'''''-' 



f / \ 



Wq 



W 




+ \ Wi 




■00 



Figure 5.6. The jump conditions satisfied by the matrix 0°^^{w) that we will show is 
normalized in the stronger sense that 6°"*(w) = I + 0(1^1-1/2) 

as w ^ 00. The only other 
singularities admitted are 0{\w — w^]"^^'^) near the points Wq and Wq. 

A standard Liouville argument shows that these two Riemann-Hilbert problems have at most one solution. 
Existence of a unique solution may be accomplished by explicit construction involving Riemann 0-functions 
of genus one. For us to be able to continue the current line of argument it is sufficient to state the following 
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proposition, whose proof can be found in Appendix [B] Recall that in using the outer parametrix we will be 
setting V = ^a/^n- 



Proposition 5.2. Riemann-Hilbert Problems \5.1\ and \5.S\ each have a unique solution with the following 
properties: 

• 0°"*(z«) depends continuously on v and remains uniformly bounded (despite having no limit) as 
V oo for w &C\{UoU Ui). 

• \w — Wk\^^^0°^^*'{w) is uniformly bounded as v ^ oo for w G Uk, fc = 0, 1. 

• For all w where 0°"*(w) is defined, det(0°"*(u')) = 1. 

Whether it solves Riemann-Hilbert Problem \5.1\ ( case L) or Riemann-Hilbert Problem 5.2 ( case R /, the matrix 
0°"'(w) has the following asymptotic forms: 

6°"*H = 6"-" + 6"'V^ + 0(w), w->0 (5.25) 

and 

Qooa 

0°"*(w) =1+ ^= +0(^-1), w^oo, (5.26) 



and the matrix coefficients O"'*^, O*^'^, and 0°°'^ are all uniformly bounded as v ^ oo. 

In terms of the matrix elements of the coefficients O"'", 0°'^, and 0°°'^ we now define the following 
quantities: 

C := (-l)#^0?i°, (5.27) 

S := (-1)#^02T, (5-28) 

and 

•,o,o^-lQO,l 



~^12 ^22 ^12 ^12 ^22 ' 



(5.29) 



where the second line follows from the first because according to Proposition 5.2 det(0°''') = det(O°"'(0)) 
1. 

In Appendix |b] the following simple formulae for C, S, and G are established. 

Proposition 5.3. Let v ~ ^aI^n — ^{x, t) /e jy + irf^A, and let the contour /3 depend on {x,t) G M.' 
described in Proposition J^.l or Proposition 4-9 //0°"*(w) is the solution of Riemann-Hilbert Problem 
(case L), the derived quantities defined by ( |5.27[ ) ~(5.29) are given by 



5.1 



■ ■ f 2^K{m) 
C — CN(x,t) = an \ — ;m 



S — SN{x,t) 



(5.30) 



. • , , 4X(to)9$ 
G = GN[X,t) = —\/mcn 



TT dt 



where the elliptic parameter is 



1 



m = TOL := sin(C) , Q < C, := - arg(u;o) < 



(5.31) 



which coincides with the function of E given by (1.62). On the other hand, if 0°"'(w) is the solution of 
Riemann-Hilbert Problem 5.2. the derived quantities defined by (5.27 1 -(5.29) are given by 

{2^K{m) 



C = Cn{x, t) = cn 
5 — Sn{x, t) - 



/2^K{m) 



(5.32) 



G = GN{x,t) = - 



AK{m)d<i>, /2$i^(m) 

TT Ot \ TTEtv 
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whe 



m = TTiR := 



(5.33) 



which coincides with the function of £ given by (1.69). Here K{-) is the complete elliptic integral of the first 



kind as defined by (I.6OI, and both elliptic parameters correspond to the so-called normal case of < m < 1. 



In both cases, the quantities Cn, Sn, and Gn are all periodic in the fast phase variable ^/^n with period 
27r, and the differential relations 



\^Cf,{x,t)GN(x,t)^0{tN) and eN^{x,t) = -^SNix,t)GN{x,t) + 0{eN) (5.34) 



(where m is a function of x and t through the roots Wq ~ Wo{x,t) and Wi 
uniformly for bounded {x,t). 



wi{x,t) of R{w;p,c\)'^) hold 



Now we describe the inner parametrices 0™{w), which are constructed in a fairly standard way from Airy 
functions. The use of such "Airy" parametrices has been a hnchpin of many papers using the Deift-Zhou 
methodology going back to the original reference . A reference that describes a similar construction as in 
the present case where the functions T^(u') and T'^(u') appear in the jump conditions is H]. The basic Airy 
parametrix is the solution of the following Riemann-Hilbert Problem. 




Figure 5.7. The jump matrix VAiry(C) defined on the contour Sai 



Riemann-Hilbert Problem 5.3 (Airy parametrix). Consider the contour SAiry illustrated in Figure 5.7 
consisting of four rays with angles arg((') = 0, arg(C) — ±2tt/3, and arg(— C) ~ 0. Find a 2 x 2 matrix 
function Z(^) with the following properties: 

Analyticity: Z(C) is an analytic function 0/ G C \ SAiry and Holder-^ continuous for any 7 < 1 
in each sector of analyticity. 

Jump condition: The boundary values taken by Z(C) on the rays 0/ SAiry are related by the jump 
condition Z-|_(^) = Z_(^)VAiry(0 /"^ C ^ ^Airy; whcrc the jump matrix VAiry(C) '■s as defined in 
Figure \5. 7\ 
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Normalization: Z(C) satisfies the normalization condition 



lim Z(C)UC 



-T3/4 _ 



uniformly with respect to direction, where U is the unitary matrix 

gi7r/4 g-i7r/4 ■ 



u 



1 
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(5.35) 



(5.36) 



It is well-documented [H [2] that this problem has a unique solution given by the following formulae. Let 
r := (|)2/3C. Then 

-3/6 ^g-3.^/4Ai'(^) gll-Vl2Ai'(re-2"/3)] 2r3/V,/3 

e-"/4Ai(r) e"/i2Ai(re"2'^*/3) 



Z(C) 
Z(C) :=y2^ 
Z(C) := \/2^ 



< arg(C) < 



27r 



(5.37) 



(T3/6 



o"3/6 r 



e-7"/12Ai(Te2^'/3) g^</12Ai(^g-2^»/3) 

jll7rj/12Ai'(T-g-27ri/3^ g77ri/12Ai'(T-g2Tri/3y 
g7r»/12 Ai(^e-27rV3) g57rz/12 Ai(^g27rz/3) ^ 



.2x3/^3/3 



V3/3 



27r 



<arg(C)<7r, (5.38) 



Z(C) 



0-3/6 



-3.rV4Ai'(^) g7^z/12Ai'(^e27r»/3y 

-"/4Ai(T) e5'^^/i2Ai(Te2"/3) 



2r3/2<T3/3 



27r 

T 



-TT < arg(C) < 



< arg(C) < 0. 



2tt 



(5.39) 
(5.40) 



Moreover, it is easy to show from standard asymptotic formulae for the Airy function Ai(-) that the normal- 

o. (5.41) 



ization condition (5.35) holds in the stronger sense that 

-0(^-3/2) 0(^-1) 



Z(C)UC 



-0-3/4 _ 



c 



To construct 0^^{w) from Z(C), we must consider two cases, depending upon whether the point w = Wk 
lies in or E^. We assume that the corresponding disc Uk is small enough that all of Uk n (E^ U E'^) lies 
also in E^ or E^ respectively. 

If Wk & E"^, then in the adjacent contour 7 the function satisfies 3?{(/)(^)} < and 0'(^) — R{£,)H{£_) 
where the quadratic R{^)^ has a simple zero at Wk and H{£,) is analytic and bounded away from zero near 
Wk- It follows that the function defined by taking the principal branch in the formula 

Wv{w):={(f,{wk)-(f>{wjf/\ we 7 (5.42) 

(recall that (j>{wk) = if Wfc € M, while (j>{wk) — ±i<I' if Wk € C±) can be analytically continued from -fOUk 
to the full neighborhood Uk, and (by taking Uk sufficiently small, but independent of ejv) Wi^{w) 7^ for 
w G Uk- Thus, W = Wy{w) defines a conformal map taking Uk to a neighborhood of the origin in the 
ly-plane. At this point, we choose the parts of the four contour arcs meeting at Wk within Uk so that their 
images under Wy are straight segments with angles arg(Vl^) = (for Wx/i.'j)), arg(— W^) — (for W-^{l3)), 
and arg(iy) — ±27r/3 (for W^{dA n i^Y))- To tie the independent variable C of Z(C) to w, we set 



Wyjw) 

2/3 ' 



Wk e E^. 



(5.43) 



d^{w) :-- 



Thus, the disc Uk is mapped under it; 1— C to a neighborhood of the origin in the C-plane whose outer 
boundary is expanding at the uniform rate of e^^^'^ as eat i 0. Set cv := e"/''*'*^/^^'^"^ if Wk € C± and 
c-v '-— if Wk S K, and let a nonzero piecewise analytic function d^{w) be defined in Uk as follows: 

'1, \aTg{-Wv{w))\<7T/3 
T^{wY/^, |arg(iyv(u;))| <2V3. 

It is then straightforward to verify that if H^(?i;) is any matrix that is holomorphic for w £ Uk, with C, 
defined in terms of w by (5.43), 

6i"(u;) H^(u;)Z(C)(-^al)4=^d^(«;)^^ w G Uk, Wk G E^, 
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(5.44) 



(5.45) 



is analytic exactly where 0{w) is and satisfies exactly the same jump conditions as 0(u;) does within the 
neighborhood Uk- It remains to determine the holomorphic prefactor H^(w), and this is done to achieve 
accurate matching onto the outer parametrix ©""'(w) on the disc boundary dUk- To do this, we now 
observe that the unimodular matrices 0°"'(u;)cy'^^ (iai) and W^{wy'^^'^\J^ satisfy the same analyticity and 
jump conditions for w ^ U and grow at the same rate as w Wk', therefore, their matrix ratio B^(i/7) 
is unimodular and analytic in Uk, and moreover, it is a consequence of Proposition |5.2| that this ratio is 
uniformly bounded as ejy i 0. We then set 

-1 



■■{lai) 



(5.46) 



U^iw) ■.= B^{w)e%''\ where (w) 
Then, it is a direct matter to check that 

6i"H6™*(w)-i = c^Hr"^/4(2(^)uc-"^/4)c"^/4c^H-iD^(u;), w e Uk, Wk e S^, (5.47) 

where 

C^(u.) 6°"*(w)cv"'(iai)U and D^(u;) := 6°"*(w)d^(w)'"^'6°"*(w)-i. (5.48) 

guarantees that 0°"'(^) and its inverse are bounded on dUk 



5.2 



Now if w = ^ G dUk, then Proposition 

uniformly as e^v i 0, while Proposition 3.1 guarantees that d^{w)'^^ = I + C'(ejv) holds uniformly for w £ Uk, 
so it follows that C^(0 = C(l), C^'TO"^ = and D^(^) = I + O(eAr) are uniform estimates for 

^ e dUk- Furthermore, since w — ^ G dUk is equivalent to — C'(e^'^), we may use the large-C asymptotic 
formula ( [Kilj ) to obtain (-'^3/4.(^2{0lJC^/^)C^^'^ = 1 + 0{C^/^) = 1 + 0(6^) for w = ^ e dUk- The result 
of these calculations is the uniform estimate 

Ojr(C)O°"*(0"' =II + 0(eiv), ^edUk, fc = 0,l. (5.49) 

This relation shows that our choice of the holomorphic prefactor H^(ii;) yields an accurate match between 
the inner and outer parametrices on the disc boundary dUk- 

If instead Wk & S^, then in the adjacent contour 7 the function 0(^) satisfies 5R{(/)(^)} > so the correct 
conformal mapping is defined by taking the principal branch in 

Wa{w) -.^ {^{w) ~ <j){wk))^^^ we-/ (5.50) 

and analytically continuing from ^ DUk to Uk (using the fact that 4>{w) — 4>{wk) behaves as {w — Wk)^^"^ 
near w = Wk)- Choosing the contours within Uk so that their images under W = Wa{w) lie along rays with 
angles arg(W^) = (for I^a(7)), arg(-t^) = (for WAil3)), and arg(W^) = ±27r/3 (for WaIOA n n^)), we 
choose the independent variable ^ in the Airy parametrix Z(^) to be given by 

WAiw) 



C 



2/3 



Setting CA := e-"/4±**^/(2^") if Wk e 

d'^iw) := < 



Z± and Ca 



T^iwy 



-1/2 



Wk e S^. (5.51) 
-^'74 jf g and defining d^{w) for w e Uk hj 

(5.52) 



\argi-WA{w))\<Tr/3 
\aTg{WA{w))\ < 27r/3, 



we define the inner parametrix as 

6i^{w) :=H^(«;)Z(C)c^='d^(w)^ 



w e Uk, Wk € T,^, 

where C is a function of w and e^v by (5.51), and the holomorphic prefactor H'^(w) is given by 



U^iw) := B^{w)e 



CT3/6 

N ' 



where B'^( w] 



'Hca"' 



WAiwy/^V^ 



(5.53) 



(5.54) 



It is again straightforward to confirm that for w e Uk, d^^{w) is analytic where 0{w) is, and satisfies exactly 
the same jump conditions as does 0{w). Also, 

6i"(w)6™'(w)-i =C'^(u;)r"^/4(Z(C)Ur"3/4)C"^/4c^(w)-iD^(w;), w£Uk, Wfe G S^, (5.55) 
where 

6™*(w)c;^''^U and n^{w) ■.^6°'-'\w)d^{w)'''6°''\w)-\ (5.56) 



C^{w) 
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Completely analogous reasoning as in the case that Wk E then shows that the estimate (5.49) holds also 
when Wk & S^. 

We formalize these results in the following proposition. 

Proposition 5.4. Suppose that the root Wk of the quadratic i?(w;p,q)^ is hounded away from the other 
root and also from a and b, and let the inner parametrix 0™(w) he defined for w in a suitahly small (hut 
independent ofeiq) neighborhood Uk ofwk by either (5.45)-( 5.46[ ) (ifwk G Tj^ ) or (5.53)-(5.54) (ifwk € S'^J. 
Then det{0^^{w)) = 1 where defined, 0{w)0™{w)~^ is analytic for w G Uk, and the mismatch with the outer 



parametrix on dUk is characterized by the estimate (5.49) 



We emphasize that it is a consequence of the differential identities ( |4.13 ) and (4.18) along with the fact 
that H is bounded away from zero near each of the distinct roots of R{w;p, q)^ that the matrix Z(C) can 
be used to construct the correct inner parametrix for 0{w) in Uk- Indeed, were the function (/'(w) — (/>(wfc) 
to vanish to higher order due to coalescence of roots of R{w;p, q)^ or the presence of a zero of the analytic 
function H, a more exotic inner parametrix would be required because W^{w) or Wa{w) would fail to be 
a proper conformal map. A more fruitful approach in such a situation is to investigate the double-scaling 
limit where </> degenerates due to allowing {x, t) to converge to some critical point at an appropriate rate 
as eN i 0. We will carry out such analysis in a forthcoming paper [5. for the case when the two roots of 
R{w;p, q)^ coalesce on the real axis when t = and x — ±a;crit- Modified inner parametrices will als o be 
required if one or the other roots of R{w;p, q)^ "bounces off" of a or b as explained in Proposition 
this is a far less interesting special case. 



4.9 



but 



5.3. The efTect of conjugation. Estimation of the error. Consider the matrix E{w) (the error in 
approximating 0{w) with the global parametrix 0{w)) defined by 

E(w) := 0{w)6{w)-^ (5.57) 

for all w E C where both matrices on the right-hand side are well-defined. Since according to Propositions |5.2| 
and 5.4 the outer parametrix 0{w) defined by ( |5.17[ ) is a unimodular sectionally analytic matrix function, it 
follows that the same is true of E(i(;). Also, since (i) the outer parametrix 0°^^{w) satisfies exactly the same 
jump condition as does 0{w) on the arcs of the contour /? and (u) the inner parametrices 0™(?«) satisfy 
exactly the same jump conditions as does 0{w) on all contours within the open discs Uk, the error E(k;) 
can be analytically continued to all of these contour arcs. Thus, the jump contour for E(i«), denoted Ee, is 



as illustrated in Figures 5.8 5.11 




Figure 5.8. The contour Ee of discontinuity of the sectionally analytic function E{w) in 
case L. The circles are the boundaries of the discs Ui and U2- 
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Figure 5.9. The contour Ee of discontinuity of the sectionally analytic function E(w) in 
case R with either A — $ or V — The circles are the boundaries of the discs Ui and U2- 




Figure 5.10. The contour Eg of discontinuity of the sectionally analytic function 'Ei[w) in 

The circles are the boundaries of the discs Ui 



case R with either A 
and U2- 



P^R or V = P, 



N 



Now while the global parametrix 0(w) is known, the matrix 0{w) is only characterized by being related 
via explicit transformations to H(?i'), which in turn is specified only as the (unknown) solution of Riemann- 
Hilbert Problem |2.1[ However, the conditions of Riemann-Hilbert Problem |2 . 1 1 imply an equivalent Riemann- 
Hilbert problem whose solution must give E(i(;). As both 0{w) and 0{w) tend to the identity matrix as 
w — ^ oo (in the case of 0{w) this follows from the sequence of explicit transformations relating it back to 



VL{w) and the normalization condition (2.6) for the latter, and in the case of 0{w) this follows from the 
0°^^{w) for large \w\ and from the normalization condition on the outer parametrix as 

I as 



fact that d{w) 



specified from the conditions of Riemann-Hilbert Problem 5.1 
w — >■ oo. To formulate the Riemann-Hilbert problem for the error it remains to analyze the jump conditions 



5.2), we must require that 'Ei{w) 



satisfied by 'E,{w) along the contour Eg pictured in Figures 5.8 -5.11 



C Eg. On eit her sid e of we have 
we therefore 



5.2 



First consider the jump of £(?«) across the positive real axis ^ e 
O-i-(^) — 05-"'(i^), and according to the jump conditions of Riemam 
have the exact relation 0-)_(^) — cr20-(C)o'2 for ^ G According to Proposition l5jj the corresponding 
boundary values of 0(w) are related by 04.(^) ~ (720- (C)'''2(I + C'(eAr)), where the error term is identically 
zero except in the interval J where the lens A abuts the positive real axis from above and below. Using the 
fact (see Proposition 5.2) that 0-(C) = 0°"*(^) is, along with its inverse, uniformly bounded for ^ > 0, we 
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Figure 5.11. The contour "E-e of discontinuity of the sectionally analytic function 'E{w) in 
case R with either A = or V = . The circles are the boundaries of the discs Ui 

and U2- 



then see that E+(0 = cr2E_(0CT2(I + X(^)) holds for ^ > 0, where X(^) = Oie^) for ^ € J and otherwise 

x(C)^o. 

Across the disc boundaries dUk C Se we have no discontinuity of 0{w), but E(u;) is discontinuous 
because the global parametrix 0{w) is discontinuous due to the mismatch between the outer and inner 
parametrices. If we take the disc boundaries to be oriented in the cl ockw ise direction, then we have E+(^) = 



E_(^)[0)f(^)0°"*(^) ^] for ^ G dllki and according to Proposition 5.4, this jump condition can be written 
in the form E+(C) = E„(0(I + ©(ew))- 



Finally, consider C e Se \ {^+ U dU^ U dUi). According to Proposition [5T] we have 0+(^) = O_(0(I + 
©(eAf)) holding uniformly for all such ^. But the global parametrix has no jump, and we have O(^) = 0°"'(^) 
on these contour arcs, and again recalling Proposition |5.2| as ejv i the global parametrix is uniformly 
bounded here along with its inverse. Therefore E+(^) = E_(06°"'(0(I + 0(eAr))6°"*(^)^i = E_(0(I + 
0{eM)) holds uniformly for ^ G Se \ (K+ U dU^ U dUi). 

It follows from these considerations that 'Ei{w) may be characterized as the solution of the following 
Riemann-Hilbert problem. 

Riemann-Hilbert Problem 5.4 (Error). Seek a 2 x 2 matrix E(w) with the following properties: 

Analyticity: E{w) is analytic for w G C\I]e where Se is the contour (independent ofej^j) pictured 
in various cases in Figures \5. 8\]5. 1 1\ and in each component of the domain of analyticity is uniformly 
Holder continuous with any exponent 7 < 1 . 

Jump conditions: The boundary values taken by E(w) on Ee o,re related as follows. For G M_|_, 

E+(e) = rT2E_(0fT2(I + X(0) (5.58) 

where X(^) — ©(ejv) for ^ G J and X(^) — for ^ G \ J. For all remaining ^ G Se we have the 
uniform estimate 

E+(0 = E_(0(I + O(ejv)) (5.59) 

as Eat \. 0. 

Normalization: The matrix E(z«) satisfies the condition 

lim E(w) = I. (5.60) 

I arg(-M))|<7r 



This Riemann-Hilbert problem closely resembles a problem of "small-norm" type, except for the form of 
the jump condition along the positive real axis (the branch cut of ^/—w). But if we consider unfolding the 
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branch cut by setting w — z'^ and then definmg a matrix function F(z) in terms of by 



F(z) 




> 
< 0, 



(5.61) 



then the jump contour Ef for F(z) in the z-plane includes two disjoint images of Se \ one in the upper 
half-plane and one in the lower half-plane; moreover from ( |5.59 ) the jump of F{z) across the arcs of either 
of these two images is of the form F_|_(2;) = F_(z)(I + ©(e^v))- If z € K, then it follows from (5.581 that 
F+(z) — F_(z)(I -I- 0{eN)) where the error term vanishes identically if ^ J. Thus, F(z) is the solution 
of a "small-norm" Riemann-Hilbert problem of standard form. It is a metatheorem in this subject that 
such problems have unique solutions that are uniformly close to the identity matrix on compact sets that 
avoid the jump contour and in a full neighborhood of the point at infinity. Indeed, solving the singular 
integral equations corresponding to such a problem involves inverting an operator (on, say, L^(Ef)) that is 
a perturbation of the identity of size 0{e]\[) in operator norm. Such a problem can of course be solved by 
iteration, and the resulting Neumann series also functions as an asymptotic series as ejv i 0. This yields a 
representation of F(z) in terms of a Cauchy integral: 



F(z) 



1 

27ri 



(5.62) 



where Y e L^j-j]^) ^j^h ||Y(- 
Y e LI(Sf) with ||Y(-)| 



©(e^r). Note that since Ep is compact, by Cauchy-Schwarz we also have 



©(eat). It follows that if ii' is a compact subset of ( 



sup 



\F{z) 



0{eN) 



disjoint from Ep, then 

(5.63) 



because the Cauchy kernel (^ — z)^^ is uniformly bounded for ^ e Ef and z E K. Also, if z lies outside of a 
sufficiently large disc containing Ef, then the geometric series (^ — z)~^ = — (z~^ + (.z^"^ + ^'^z^^ + 
uniformly convergent for ^ G Ef, and so we obtain the convergent series expansion for F(z) as z — > o 



IS 



F(z) 



n=l 



■ f Yi^^-U^, (5.64) 

and the coefficient of each negative power of z is ©(e^r). Corresponding results hold for E(it;) by restricting 
z to the upper half-plane and using (5.611. 

We have thus shown that the significance of the global parametrix 0{w) deffiied by (5.171 is the following 
approximation result. 



Proposition 5.5. 



Suppose {x,t) is a point in one of the domains , 0"^ 



(see Proposition or 

(see Proposition 4-9), and that \t\ is sufficiently small. If also the roots of the quadratic R{w; p, q)^ do not 
coincide, nor does either root equal a or b, then Riemann-Hilbert Problem 2.1 has a unique solution H(w), 
and the matrix 0{w) obtained therefrom by means of the systematic substitutions H i— )■ J (see (2.32)J, 



J H- !■ M (see ( |3.6[ ) ), M M- N (see (3.43)j, and N M- O (see (5.1 )j has expansions for large and small w of 
the form 

(5.65) 



0{w) = 0^jf{x,t) + 0'^'(x,t)V-w + 0(u;), w^O 



,0,1 



and 



0{w) 



0^'^{x,t) 



+ C'(w"^), w 



oo, 



and the coefficients satisfy the estimates (see Proposition 5.2) 

50,0 

00,1 



0°M°{x,t) 



'N 
)0,1 

'n 

\oo,l 



O^T {x,t) 



O(ejv) 
O(ejv) 



(5.66) 



(5.67) 



where the dependence on N and {x,t) on the right-hand side enters through v — ^{x,t)/e^ + tt^A and the 
motion of the contour /3. The ©(ejv) error terms are also uniform with respect to {x,t) as long as the roots 
of the quadratic R{w; p, q)^ are bounded away from each other and from a and b. 
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If w lies in a small neighborhood of the origin, or alternatively if |w| is sufficiently large, then according 
to (5.1) 0{w) coincides with N(w), and according to (3.6) M(ii;) coincides with J{w), the latter matrix 
function being the solution of Riemann-Hilbert Problem |2.2[ Recalling the relations ( |2.32 ) and (3.43), we 
see that for such w, 

-w + \/—y 



H(u;) = 0{w)e 



Since .g(0) — g{oo) — 0, the diagonal factor relating H(u') and 0{w) has the expansions 




o9(-"')o-3/e]V 






oo 



(5.68) 



(5.69) 



for some diagonal matrices (x, t) and C^' {x, t). The matrices Aj^{x, t), B^(a;, t), and B^(a;, t) obtained 



from H(u') via ( 2.9 )-( 2.11 1 are then written in terms of the expansion coefficients of 0(w) from ( 5.65 )-( 5.66) 



as 



AN{x,t)^{-l)*^0%"{x,t) 

B%(x,t) - {-l)*^C'^\x,t) + 0%"{x,t)-^0%\x,t) 
B^ix,t) = 0^'\x,t) 



(5.70) 



Then using Proposition 5.5 these are expressed asymptotically in terms of the corresponding expansion 
coefRcients of the outer parametrix ©""'(w) as follows: 

ANix,t)^i^l)*^6"-° + 0{eN) 

B%{x,t) = {~i)*^c%\x,t) + (6"-°)-i6"-i 

B^{x, t) = 6°°^i + C^'\x, t) + 0{€n) 
(recall that according to Proposition 5.2 the coefficients 0°'*^, 



0{eM) 



(5.71) 



0°'^, and 0°°'i are bounded as cm i 0, 
along with (0"'°)~^). Now we recall the definitions of the quantities cos(\uN{x^t)) and sin(iwjv(a;, t)) (see 
(2.15)), and of eNUN^t{x,t) (see (2.16l) characterizing the fluxon condensate {wAr(a;, t)}^^jVo according to 



Definition 2.1 and compare with the definitions (5.27|-(5.29p to obtain the asymptotic formulae 
cos ^UN{x,t)^ — C]y{x,t) + 0{e]y) and 

and 



sin ( ^UN{x,t) 



du 

eN—Q^{x,t) = GN{x,t) +0[eN), 



(5.72) 



(5.73) 



where the asymptotics are valid for the same ranges of [x, t) and with the same nature of convergence as 
in the statement of Proposition 5.5 The asymptotic formulae (5.72 1 are differentiable with respect to t 



(yielding (5.73l) according to (5.34) from Proposition |5.3| 



Appendix A. Proofs of Propositions Concerning Initial Data 

A.l. Proof of Proposition Given a real-valued differentiable function f{v) defined on < u < T^, 
define 



f{v)dv 



-V <w <Q 



(A.l) 



as the right-hand side of (1.46) with f'{v) = (p{v). If / is in the range of (1.45) for some G satisfying 
Assumptions 1.2 and 1.3 then we know that V = —G{0) and 



ds 



< w < y = -G(o). 



(A.2) 
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Therefore in this case we have 

J Jq 



2v 



dsdv, G(0)<u><0. 



Exchanging the order of integration and setting t = yields 

mM = - / , , 

The inner integral evaluates to tt (independent of s and w) so 

rG-Hw) 

/[/]H= / ds = G~\w) 



■ ds. 



(A.3) 



(A.4) 



(A.5) 



yielding the identity (1.46) as desired 



A. 2. Proof of Proposition 1.2, First note that analyticity and strict monotonicity of G{x) automatically 



ensures the analyticity and positivity of in the open interval of its definition. Now assume that < f < 
— G(0). Introducing m = G(s)^ as a change of variables in (1.45) yields 



1 



G(0)^ 



G(0)2 - m 



(A.6) 



where ^ is defined in terms of G by (1.481 and satisfies the conditions of Assumption 1.4 Let S{m]v) be 
the analytic function of m for m G C \ [v^, G'(O)^] that satisfies 

9 

Sim; vf = ^, ,„ and lim S(m + iS; v) > for < to < G(Of. (A. 7) 

G(0)2 - TO sia ^ ' ^ ' 



Then, with Li being the contour loop shown in Figure pV.l[ 



*(A) ^ \ i Sim;v)^{m) — . 
4 /li m 



Since S is also analytic as a function of v when to lies on Li and v lies in the open region bounded by 







Figure A.l. The loop contour Li surrounds the branch cut [w^,G(0)^] ofS{m;v) while the 
loop contour L2 also encloses the origin. Both contours lie within the presumed domain of 
analyticity of'i^{m) (shaded). 



Li, and since S is continuous in ni for such v and to, this formula immediately shows that ^'(A) is analytic 
for each v E (0,— G(0)). To analyze the behavior near the endpoints, we proceed as follows. By a simple 
contour deformation, 



*(A) = — 5(0;i;)^(0) + - / S{m;v)^im) — 
2 4/^3 to 



(A.9) 



where we used the fact that S{0,v) = —iv/G{0). Now S(m;v) has the following convergent expansions: 

(2fc- l)!!w2fc 



5'(to; v) = S{m; 0) 



21 1/2 

TO 



S'(to,0) 



1-E 



(2fc)!! TO^ 



|m| > lul"^. 



(A.IO) 
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Sim;v) = S'(m;-G(0)) 



S{m;-G{0)) 



-2G(0)(v + G(0)) + (w + G(0))2 



{2k - 1) 



fc=i 



G(0)2 
-2G(0)(w 



1/2 



G(0)) + (« + G(0))^ 



S'(to;-G(0)) 



(2fc)!! 
G(0)(i; + G(0)) 



G(0)2 



G(0)2 
G(Of| 



> 



2G(0)(w 



G(0)) + (^; + G(0))2|. 

(A.ll) 

Note also that S'(to; — G(0)) = i. Both of these expansions are uniformly convergent on the contour L2 if v 
is confined to a sufficiently small neighborhood of = or u = — G(0) respectively, and hence the integral 



in (A.9) may be calculated term-by-term. In the case of the expansion for v small, the result is a convergent 



series in even nonnegative powers of v with purely imaginary coefficients. In the case of the expansion for 

V near — G(0), the result is a convergent series in (generally) all nonnegative integer powers oi v + G(0), 
again with purely imaginary coefhcients. Since power series always converge in disks, we now see that we 
have constructed the analytic continuation of \E'(A) valid for v in full complex neighborhoods of w = and 

V = — G(0) respectively. The form of the Taylor series (1.52) is now clear, but it remains to confirm the 
positivity of a and the value of ^'(0). But from (A.9), 



27r^(0) 
' G(0) 



> 



(A.12) 



since G(0) < and J#(0) > by hypothesis, and the constant term ^'(0) 
directly by passing to the limit A 



confirm the relations (1.51), one may pass to the limit A 



G||i/4 is easier to evaluate 
in the formula (1.45) than by working with the series. Finally, to 
iG(0)/4 in (fOsl) to obtain «'(-iG(0)/4) = 0, 



and then also from (A.9) and (A.ll) 
d 



dv 



*(A) 



7r^(0) zG(0) 
2G(0) ^ 4 



^(m) dm 7r^(G(0)2 



G(0)2 



2G(0) 



< 



(A.13) 



since ^{G{0)^) > by hypothesis, where the integral over L2 is evaluated explicitly by residues. 



Appendix B. Details of the Outer Parametrix in Cases L and R 

Since many of the important details of the asymptotic behavior of the sine-Gordon equation are de- 
rived from the appropriate outer parametrix, we here provide all details of the solution of Riemann-Hilbert 
Problems |5.1| and |5.2| in terms of Riemann 0-functions of genus one. We also explain how the extracted 
potentials (approximate solutions of sine-Gordon) can be reduced to a very simple form in terms of Jacobi 
elliptic functions. This appendix contains all details of the proofs of Propositions |5.2| and |5.3| 



B.l. The outer parametrix in case L. Proof of Proposition 5.2 in this case. 



B.1.1. Solution of Riemann-Hilbe rt P roblem \5A\ in terms of Baker- Akhiezer functions. The first step in 
solving Riemann-Hilbert Problem 5.1 for ©""'(w) is to introduce a new, equivalent, unknown P(w), given 
in terms of 0°"*(tx;) by 



P{w) 



0°"*(u;), 

(726™'(w)a2, 



Q{w} > 
^{w} < 0. 



(B.l) 



The equivalent Riemann-Hilbert problem satisfied by P(w) is described by the scheme shown in Figure B.l 
As was the case with 0°"* (w), the boundary values taken on the two disjoint components of the jump contour 
are continuous and bounded with the exception of the endpoints wq and Wq where inverse fourth roots are 
tolerated. 

Next, we remove the real parameter v from the jump conditions by defining the scalar function h{w) as 
follows: 

s{w) r" ds 

2TTi J^, S+{s){s-w)' 
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h{w) 



(B.2) 



P+(0 = P-(0^aie--3 




P+(0 = P-(0^aie^' 



Figure B.l. The jump conditions satisfied by the matrix P{w) normalized as P{w) = 
I + 0(|w|~i/2) asw-^oo. 



where the integration is along the jump contour shown in Figure B.l 
S{w) is analytic in the complement of the jump contours, and S{w) - 
branch of w^/^). Note that 

h{w) = pw^^'^ + 0{\w 

where 

1 

P 



and where S{w)'^ := w{w — Wo){w — Wq), 
= w^^^{l+0{w~^)) as ui — >■ CO (principal 



The function defined by (B.2| satisfies 



2ni 



ds 



w — ?> oo, 



S+is) 



(B.3) 
(B.4) 



for ^ on the negative real axis, and + 



— 1 for ^ on the contour connecting Wq and wq- It takes continuous and bounded boundary values on the 
entire jump contour. The new unknown we define in place of P{w) is then 

Qiw) ■.= P{w)e"'^^'"^'''. (B.5) 



Direct calculations then show that the conditions determining Q(io) are as indicated in Figure B.2 Inverse 



Q+(e) = Q-(0^^i 




Q+(0 = Q-(0^^i 



Figure B.2. The jump conditions satisfied by the matrix Q(it;) normalized as 
Q(u;)e-*P''"'''''"3 = I + 0(1^1-1/2) as u; ^ oo. 



fourth root singularities are again admitted at wq and Wq. 

The next transformation is undertaken to diagonalize the prefactor of a2 in the jump condition on the 
negative real axis. So, since 

giTr/4 gi7r/4 



t72 = VCT3V 

the matrix defined by 



V 



V2 



det(V) = l, V-i=Vt, 



(B.6) 



R(w) := V^Q(u;) (B.7) 

satisfies the conditions indicated in Figure [B. 3 [ Inverse fourth root singularities are once again admitted at 
Wo and Wq. 
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R+(0=R-(0^ai 




R+(0 = R-(0^^i 



Figure B.3. The jump conditions satisfied by the matrix R(?i;) normalized as 
R(u;)e-^P'^"''^'^^' = Vt + 0(1^1-1/2) as ix; ^ oo. 



The next transformations aim to convert the post-multiphcative jump matrices both into the permutation 
matrix ai. To accomphsh this on the arc connecting Wq to wo, we introduce the function q{w) satisfying 

= (B.8) 

W — Wq 

that is uniquely specified as being analytic except on the contour arc connecting Wq to Wq and satisfying 
q{w) = 1 + 0{w^^) as w — >■ oo. This function satisfies q+{w) = iq^{w) for w on the contour arc of 
discontinuity. Setting 

S(w) g(w)-iR(w), (B.9) 



we find that the new unknown S{w) corresponds to the conditions in Figure B.4 Note that multiplication 



s+(0 = s-(0^i 



S+(0 = ^3S-(0^2 




S+(e) = S_(Oai 



Figure B.4. The jump conditions satisfied by the matrix S{w) normalized 



as 



+ Oilwl-^/"^) as w ^ oo. Unlike 6°^'(u)), P{w), Cl{w), and R(w), 



the matrix S{w) is required to be bounded in a neighborhood of w ~ Wq while we allow a 
stronger singularity at w — wq: S{w) = 0{\w — wo\^^^^)- 



by q{w) changes the nature of the admissible singularities at wq and Wq: S{w) is now required to be bounded 
at Wq, and an inverse square root singularity is admitted at wq. 

To convert the jump condition on the negative real axis to the same form, we first separate the two rows 
of S{w) by writing 



(B.IO) 



and we also write V — [vi,V2] to give notation for the normalized eigenvectors of (T2 (see (B.6)). As is 
diagonal, the jump conditions for the individual rows decouple and in each case these conditions only involve 
matrix multiplication on the right. Thus we have Sj+(^)^ = Sj_(^)^(Ti for j = 1,2 and ^ on the contour 
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arc connecting Wq and wo, while we have Si+(^)^ = Si_(^)^CT2 and S2+(0^ = §2- (^)^(— (72) for ^ on the 
negative real axis. Now set 

ti{wf %/2si(w)'^e-*^'=("')'"^/^ and t2(w)'^ := iV2s2{wf e'"''''^"'''''''/^ , where k{w) := ^ + h{w). 

(B.ll) 

The function k is analytic on the complement of the contours and satisfies fc+(C) + ^-(0 = on the arc 
connecting w'q with wq while k^{^) + ^-(0 — 1 fo^' negative real ^. Using (B.3I one sees that k{w) has the 
asymptotic behavior 

k{w) ^pw^^^ + ^+0{\w\-^^^), w-^00. (B.12) 
From this information it is easy to see that on all jump contours we have 

t,+(e)^ = t,_(oVi, J = 1,2, (B.13) 

and we also have the normalization conditions 

t^iivfe-'P'^^""'^'"' = [1, 1] + Oi\w\^^/^), w^oo (B.14) 

where 

TT TT 

ipi := u - - and ip2:=i^+-. (B.15) 

Both ti(w)^ and t2(w)^ may become unbounded in the finite w-plane only as w — > wq, where all four scalar 
components must be 0{\w — wqI"^^^)- 



Finally, we implement the involutive jump conditions (B.13) by viewing the elements of the row vectors 
tj{w)^ as single- valued scalar functions on an appropriate Riemann surface. Let X be the Riemann surface 
of the equation j/^ = S{w)'^ = 'w{w — 'Wo){w — Wq) compactified a,t y = w = 00. We view the finite part of 
X as two copies of the w-plane (sheets) cut along the branch cuts of the function S{'w) and glued together 
in the usual way. We define on X \ {00, wq} two scalar functions ti{P) and t2(-P) as follows: 

' \[t,(u;(F))^]2, Pe sheet 2, ^ ' 

where for row vectors [ui,M2]j := Uji and w{P) denotes the "sheet projection" function. These definitions 
are consistent along the cuts where the sheets are identified precisely because the jump matrices for tj(w)^ 
have all been reduced to the simple permutation (sheet exchange) matrix ai. The Baker-Akhiezer functions 
tj '■ X \ {wq, 00} — !■ C are analytic in their domain of definition, which omits just two points of X. Since 

fp^ f{w{P)-woy/\ Pe sheet 1 

I -(w(P) - Wo) ' , P e sheet 2 

is a holomorphic local coordinate for X near the branch point P = wq, we see that tj{P) admits a simple 
pole at this point. Near the point P = 00, tj{P) has exponential behavior: 

tj{P)e''P^^y-^P^'" =l + 0{y^{P)), P^oo (B.18) 

where 

, , \w{P)-^/'^, P e sheet 1 

2/00 P := ^ / , 1 2 , _ B.19 

I -w[P) ^1^, P e sheet 2 

is a holomorphic local coordinate for X near the branch point P ^ 00. 

B.1.2. Construction of the Baker-Akhiezer functions. To express tj{P) in terms of special functions requires 
a few ingredients. Firstly, one chooses a basis of homology cycles on X consisting of two noncontractible 
oriented closed paths, o and &, such that b intersects a exactly once, from the right of a. For later convenience, 
we suppose that neither of these paths passes through the point P = 00. The homology cycles we choose are 



illustrated in Figure B.5 As X is a genus one Riemann surface (elliptic curve), there is a one-dimensional 



space of holomorphic differentials; as a basis we choose the unique holomorphic differential on X of the form 

.{P) = c*^^ (B.20) 
^ ' S{P) 
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Figure B.5. The homology basis on the sheets of X. Thick curves indicate branch cuts 
where the two sheets are identified, solid curves are on sheet 1 and dashed curves are on 
sheet 2. 



where S{P) denotes the hft to X of the function S{w): 

^(P)^ fWP)), sheet 1 

^ ^ \-S{w{P)), Pe sheet 2, ^ ^ 

and where the constant c is chosen so that 

i uj{P) = 2TTi. (B.22) 

J a 

Denote by % the other loop integral of uj{P): 

r dw{P) 

H:^^.(P)^2..|^. (B.23) 

fa S{P) 

It is easy to see that although H is complex, ^{H} < 0. The Riemann Q-function corresponding to % is 
the entire function of z G C given by the Fourier series 

°° 1 2 

Q{z:n) := 62^" e"^ (B.24) 



n— — CO 



A simple relabeling of the sum by n i— > — n shows that 

e(-z;-H) = e(z;-H). (B.25) 
The Riemann 0-function satisfies the automorphic identities: 

e{z + 2Tri;H) ^eiz;H) (B.26) 

and 

e{z + H;H)^e-i'"e-''Q{z;H). (B.27) 
The function Q{z]'H) vanishes to first order at all points z of the form z — K, + 27rim + T-Ln and nowhere 
else, where m and n are integers and where 

/C = /CCH) iTT + i-H (B.28) 

is the Riemann constant. We choose as a base point on X the branch point Pq = w;o, and then define the 
Abel map by 

A{P)-~ / uj{P') (mod 27rirn + -Hn), m,neZ. (B.29) 

Jpo 

The value of A{P) is not completely determined only because the path is only determined modulo the cycles 
a and b. Finally, let ^{P) denote the abelian differential of the second kind with double pole at P — oo: 

niP):^''^^^^dwiP) (B.30) 
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where the constant C is chosen so that 



ln{P) = 0. 

J a 



Note that by asymptotic expansion of (B.301, we have 

Jpo Voo{P) 

Let n denote the other loop integral of Vl{P): 



K := <j>n{P). 



P 



(B.31) 



(B.32) 



(B.33) 



Lemma B.l (see [3). k — 2c. 

Proof. Let X denote the canonical dissection of X obtained by cutting X along the cycles a and b. X is 
a simply-connected complex manifold with boundary illustrated in Figure B.6 Consider the meromorphic 




Figure B.6. The canonical dissection X of the compact Riemann surface X. Note that if 
we omit the upper and right-hand boundaries, then the points of X and those of X are in 
one-to-one correspondence. 



differential defined on X given by the following expression: 

r(p) A{p)niP). 



(B.34) 



Here, the Abel mapping is well defined since X is simply connected. Let dX denote the positively-oriented 
boundary of X, and consider the integral 



/ r(p). 

JdX 



(B.35) 



We will evaluate / two different ways. On the one hand, we may evaluate / by residues. The only singularity 
of r(P) is a double pole at the point P — oo E X , and while ft{P) has no residue there, r(P) indeed has 
one: 



r(p) 



A{oo) 



u^(P') 



\v^{P)+0{yUPf) 



= [A[^) - 2cy^{P) + 0{yoo{Pf 
A{oo) 2c 



-yUP) + 0{yUPf 



-2 



dyoo{P) 

yooiPr 



(B.36) 



VooiPr yUP) 



o{i) 



dyoo{P), P^oo, 
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so by the Residue Theorem, 



J = / r(P) = inic. 
Jax 



(B.37) 



On the other hand, we may evaluate / directly: 



/ = / r(P) = / [Aight(^) - A,,i,{P)] n{P) + f [AtopiP) ~ Aottom(i')] ^{P), (B.38) 

JdX Ja Jb 

where the subscripts on A indicate where the corresponding points P live on the diagram of X shown in 
Figure [BTS] But clearly. 



bright 

Therefore, 



uj{P') = -n and Atop(P) - ^bottom(P) = f oj{P') = 2m. 



I = -n (f n{p) + 2TTi (f n{p), 

Ja Jb 



and due to the choice of the constant C, we simply find 

/ = 2TriK. 



Comparing (B.37) with (B.41), we see that n = 2c, as desired. 



It then follows from (B.4), (B.20), and the normalization condition (|B.22|) that 

1 



pK 



ids 



S+{s) 



1 

'27ri 



OJ = -1. 



(B.39) 
(B.40) 

(B.41) 

□ 

(B.42) 



We now give Krichever's formula [16] for the Baker- Akhiezer functions tj{P): 



t,{P) -.^N, 



e{A{p) + jc + ipKipf,n) 

Q{A{P) + IC;H) 



cxp ipifj 



Q{A{P) + JC-iipj-'H) 
Q{A{P) + K-n) 



exp ip^Pj 



n{p') 



n{p') 



(B.43) 



where Nj is a normalizing constant. Here, the path in the exponent is intended to be the same as that in the 
Abel map A{P). This expression is well-defined in spite of the indeterminacy of the path of integration due 
to this identification of paths and the two automorphic identities (B.26) and (B.27) satisfied by 6. Indeed: 

• Adding an a-cycle to the path does not change the exponent due to (B.31 1, and the 0-functions are 
also unchanged due to (B.22) and (B.26). 



• Adding a &-cycle to the path adds —iipj to the exponent according to (B.33I and (B.42), but this is 
compensated by the ratio of 0-functions according to (B.23) and (B.27). 
Moreover, since by our choice of base point A{P) vanishes to first order in the holomorphic local parameter 
yoiP) when P = wq, it is clear that tj{P) may have a simple pole at this point. Also, from the asymptotic 
behavior (B.32) of fl{P) it is clear that if the constant Nj is chosen correctly tj{P) will have the desired 
asymptotic behavior (B.18) as P — > oo. 

To compute Nj, let us select a path from P = to P = oo that lies entirely on sheet 1 of X and along 
which 5{i(;(P)} > 0. This unambiguously determines A(oo) as well as the constant term in the asymptotic 
expansion of the exponent. It is easy to verify for such a path that 

^M = -^j(^(P) = -^H. (B.44) 

Also, as P tends to P = oo along such a path (which we take to lie along the negative real ax;is for large w), 



r n{p') = \i> n{p') - I In{p') + H^'^^ 

Jwn ^ J a Jb Jo 



c 



25+ (w) 



dw 



— K - 
2 



w(P) 



-C 
2S+{w) 
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(B.45) 



dw. 



and 



w(P) 



C 



2S+{w) 



dw — 



w{P) r 







w{P) r 



w + C 


1 


2S+{w) 




w + C 


1 


2S+{w) 





1 



(B.46) 



dw 



vupy 



The remaining integrand is integrable at infinity, so by doubling the contour of integration along the branch 
cut on the negative real axis and then closing the contour in the right half-plane we may write 



v{P) r 



w + C 
2S+{w) 



-1/2 



dw 



C 



2S+{w) 
'w{P) + C 



-1/2 



2S{P) 
1 



dw{P) 



Therefore, 



and so we find 



fn{p') = 

J wn 



Voo{P) 



K 

2 



a VHP) 



0{y^{P)), 



dw + 0{yoo{P)) 
dw{P) + 0{Voo{.P)) 
0{yoo{P)) 



(B.47) 



lim t-fp)e-*'^^P^~(-P)"' - jv. Q(^^~^y^J'^) e»y,/2 



P^oo 



It follows that the normalization constant is simply 



N. 



1 



-iipj /2 



(B.4 



(B.49) 



(B.50) 



This completes the construction of tj{P), and hence of 0°^*(w). To obtain an especially useful formula for 
0°^*^{w), we need another intermediate result. The function 



l(w) 



P 



C 



2Siw') 



dw' 



(B.51) 



is well-defined for w in the cut plane (the first sheet of X as illustrated in Figure B.5) due to the condition 
(B.31 1, which makes the integral independent of path. 



Lemma B.2. l{w) = k{w) := h{w) + 

Proof. Clearly, I and k are analytic in the same domain, and are both uniformly bounded on bounded subsets 
of this domain. Since S changes sign across the cut connecting wq and Wq it is obvious that along this cut, 
l+{w) + l-{w) = k+{w) + k-{w) = 0. Now because the integrand of I is the restriction of the differential 



to the first sheet, it is easy to see that 1{Q) = —\pk, so by (B.42), in fact — ^. Again using the 
fact that S changes sign across the negative real axis, we see that along this cut l+{w) + l-{w) = 21(0) = 
k^{w) + k_{w) = 1. Also, both l{w) and k{w) have the same leading asymptotic behavior as w — > oo, 
namely, pw^^^ + 0{w~^^^). Setting m{w) := {l{w) — k{w))/S{w) we see that m{w) is a function that is 
analytic where S is, with at worst inverse square-root singularities at w = 0,wa,WQ. Moreover, along the 
open arcs of discontinuity of S", we have m+(w) = m_(w), so m extends continuously to these arcs. It follows 
that m(w) is necessarily an entire function of w. Then, since m{w) — 0{w^'^) as w oo, m{w) = by 
Liouville's Theorem, and the result follows. □ 

We now present a formula for 0°^^{w). For w in the cut plane, let Pk{w) denote the preimage of w under 
w{P) on sheet k of X. Then, by composing the transformations leading from 0°"'(w) to the Baker- Akhiezer 
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functions tj{P), we find that for > 0, 

1 (t2(Pi)e-''^^''-ii(Pi)e-'"^^'') f (ti(P2)e*'^^'' + i2(i'2)e*'^^''). 



0™'(w;) 
while for < 0, 

6°"*(u;) 



I (i2(P2)e'^^'' - ti(P2)e*^^") I (ti(Pi)e-^'^^'' + i2(Pi)e 



(B.52) 



(B.53) 



The dependence on w enters these formulae via q — q{w), Pj = Pj{'w), and h — h{w). 

To simplify tj{Pi{w))e~^'^^^'^^\ we evaluate the formula (B.43I by selecting in both the Abel map and 
the integral in the exponent a path from wq to Pi{w) lying entirely in the finite w-plane on sheet 1 of X. 
Therefore, using Lemma B.2 and (B.50I, 

e(A(Pi(w)) + /c - iipj^n) ,^^iM-i^^h(w) 



tj{Pi{w))e-'^^''^'"^ = K 



' e{A{Pi{w)) + ic;n) 

e{iTr;H)e{A{Pi{w)) + IC - iipj;H) 
e(i7r - iipj;n)e{A{Pi (w)) +IC;n)' 



(B.54) 



To evaluate tj{P2{w))e 



we proceed similarly, choosing a path from wq to P2(w) lying in the finite 



(B.55) 



w-plane on sheet 2 of X to obtain 

^ ei^n;H)eiA{P2{w))+JC-ivf,H) 
QiiTT ~ iifj;n)e{A{P2{w)) + IC;'Hf ' 

This latter formula can be further simplified with the observation that A{P2{w)) — ~A{Pi{w)) because the 
base point Pq has been chosen as the branch point wq, so with the use of the relations (B.25)-(B.27| and 
the definition (B.28I of /C we obtain ultimately 

tj(P2(u;))e"^^''('") = 



e{iTr;n)eiA{Pi{w))+IC + iipj;n) 
e{iTr + iiPj;n)e{A{Pi{w))+IC;n)' 

In the formulae (B.54) and ( |B.56 ), A{Pi{w)) represents any value of the integral 



(B.56) 



5.2 



This essentially completes the proof of Proposition 
follow from the formulae (B.52), (B.53), (B.54), and (|B.56|) upon noting that the only dependence on i/ 



in case L. Indeed, the uniform bounds on 0°"*(k;) 



occurs via the angles (pj defined by (B.15), and that 8 satisfies the periodicity relation (B.26). That 
det(0°"'(i«)) = 1 is a consequence of the fact that the jump matrices have determinant equal to 1, via a 
standard Liouville argument. 



B.2. The outer parametrix in case R. Proof of Proposition 5.2 in this case. 

B.2.1. Solution of Riemann-Hilbert Prohle m\5.S\ in terms of Baker- Akhiezer functions. The contour of dis- 
continuity of C)°"*(w) as illustrated in Figure p . 6| divides the complementary region into a bounded component 
To and an unbounded component Too- The first step in solving Riemann-Hilbert Problem 5.2 is to make an 
explicit substitution to simplify the contour. We therefore define a new unknown P(u') in terms of 0°"*(i/7) 

by 

'6°"'(w), > and w e Too 

6°"'(^^;)wle''""^ > and w G Tq 

a2d°'^\'w)a2, < and w G Too 

^cr20™'(«^)c^3e''"'', < and w G Tq. 
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P(w) 



(B.58) 



This substitution has the effect of collapsing the contour to the real axis and removing the jump discontinuity 
along The matrix P{w) may be analytically continued to the domain C\{—oo, 0], and its boundary values 
on the negative real axis are necessarily continuous except at the points wq and wi at which negative one- 



fourth power singularities are admitted. The jump conditions satisfied by F{w) are illustrated in Figure B.7 



P+(0 = a2P-(0(-^3) 
. \ 



P+ (0 = ^2?- (0^2 



Wo 



Wi 



Figure B.7. The jump conditions satisfied by the matrix P{w) normalized as P(w) 
I + as w OO. 



Next, we remove the real parameter v from the jump conditions by defining the scalar function h{'w) by 

. , S(w) f° ds , , 

Hw):^ ^/ r, B.59 

Jw^ S+{s){s~w) 

where S{'w)'^ :— w{w — wo){w — wi), S{w) is analytic in the complement of the real intervals — oo < w < wq 
and wi < w < 0, and S{'w) = w^^'^{l + 0{w~^)) as w — t- oo (principal branch of w'^/^). The asymptotic 
behavior of h{w) as w — > oo is given by 

h{w) = pw^/^ + 0{\w\-'^^^), w^oo, (B.60) 

where 



1 



ds 



(B.61) 



TTZ J^^ S+{s)' 

The function defined by ( B.59| ) is analytic exactly where S{w) is, and it satisfies + ^-(C) ~ for 

— oo < ^ < Wo and h+{^) + h-{£^) = —2 for wi < ^ < 0, taking continuous and bounded boundary values. 
In place of P{w) we now take the new unknown defined by 

Qiw) ■.= P{w)e"'''^'"'^'''. (B.62) 

By direct calculation we obtain the jump conditions determining Q(i(j) as shown in Figure B.8 The boundary 



OO ■ 



Q+(0 = ^2Q-(0(-^3) 

■ — \ ^ 



Q+(0 = fT2Q_(0^2 



Wo 



Wi 



Q+(0 = ^2Q-(e)(-^2) 



Figure B.8. The jump conditions satisfied by the matrix Q(w) normalized as 
Q{w)e-'P'""'"''''' = I + 0{\w\-^/^) asw~^oo. 

values taken by Q(w) are continuous except at the points wq and wi at which inverse fourth-root singularities 
are admitted. 

We next diagonalize the prefactor of a2 with the use of the eigenvector matrix V given in (B.6|. Setting 

Il{w) = V^Qiw), (B.63) 
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one finds that the jump conditions on the negative real axis are reduced to those shown in Figure [B.9[ and 



R+(0=^3R-(0(-^3) 



oo ■ 



R+(0=tX3R-(0^2 



1^0 



Wi 



R+(0 = a3R-(0(-^2) 



Figure B.9. The jump conditions satisfied by the matrix R(ti;) normalized as 



as w ^ oo. 



the only discontinuities of the boundary values are once again inverse fourth roots admitted at wq and wi. 
Now, let q{'w) be the function satisfying 



4 w-wo 
q{w) = 



(B.64) 



w ~ Wi 

and the normalization condition lunw^oc qIw) — 1, and taken to be analytic for w € C \ [wo,wi]. Its 
boundary values taken on the branch cut [wo,wi] are related by q+{S,) = —iq-{S,). Note also that 

/ \ 1/4 

g(0)=f^j >0. (B.65) 

Now define a new unknown S{w) by 

S{w) -.^qiwy^Riw). (B.66) 

Taking into account the type of singularities that R(z/;) may have near w = wq and w — wi, we see that 
S(w) = 0{\w — Wq\~^^'^) for w near wq, while S{w) is bounded near w — wi. The jump conditions satisfied 



by S{w) on the negative real axis are illustrated in Figure B.IO 



s+(0 = ^3S-(0H^3) 

. \ 



oo ■ 



S+ (0=^38.(0^2 



Wq 



Wi 



S+(0 = a3S_(0(-^2) 



Figure B.IO. The jump conditions satisfied by the matrix S(w) normalized as 
Siw)e-'P'''"'^^''^ = Vt + 0(|w|-i/2) as w ^ oo. Unlike 6°^'(w), P{w), Q{w), and K{w), 
the matrix S(w) is required to be bounded in a neighborhood of w = wi while we admit a 
stronger singularity at w = wq: S{w) = 0{\w — wo\~^^^). 

We now separate the rows of the matrix S{w) by writing 



Siw) 



and introduce two new row vectors ti(i/7)^ and t2(w)^ by setting 



(B.67) 



(B.68) 
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where k{w) is the function analytic for w € C \ IR_ given by 

kiw) := +hiw)-^ 
2 2m 



ds 



'wo S{s){s-w)' 

Note that k{w) is bounded on compact sets in the w-plane and has the asymptotic behavior 



k{w) = 



P 



1 

27ri 



ds 



,1/2 



+ -+0{\w\-^/^), W^OO. 



(B.69) 



(B.70) 



Also, the jump conditions satisfied by k on the negative real axis are as follows: + fc-(C) = 1 for 

^ < Wo, k+{^) - fc_(^) = -1 for Wo < C < wi, and k+{^) + = -1 for wi < ^ < 0. 

From this information it follows that tj{w)^ are analytic for u; G C \ ((— oo, wq] U [wi, 0]), and they satisfy 
the involutive jump conditions 

t,+(c)^ = t,-(e)Vi, J = 1,2, 

for ^ in either of the two intervals of discontinuity, and we also have the normalization conditions 
where (fix and tp2 are given by 

vr 1 ds 
ipi := v — / and Lp2 



''^■^ = [1,1] +0(1^1-1/2), 



(B.71) 
(B.72) 



4ip 



S{s) 



n 
2 



4«p 



ds 



(B.73) 



Both ti(u;)^ and t2{w)^ may become unbounded in the finite w-plane only as w — > wq, where all four scalar 
components must be 0{\w — wo\^^^^). 

We are now in a position to identify the components of the row vectors tj(w)^ as sheet projections 
of scalar Baker- Akhiezer functions tj{P) onto the two sheets of the Riemann surface X of the equation 

— S{w)'^ = w{w — wq){w — wi) compactified at ?; = w = oo. Viewing X as two copies (sheets) of the 
w-plane cut along the intervals (— oo,wo] and [wi,0] and appropriately glued together, the Baker- Akhiezer 
functions are then defined just as in case L: 

'[tj(w(F))"^]i, Pe sheet 1 
Jtj(w(P))"^]2, Pe sheet 2. 

These functions are analytic on X except at exactly two points: the branch point w — wq at which tj{P) 



(B.74) 



admits a simple pole (in the holomorphic local coordinate yo{P) given by (B.17l), and the branch point 
w = oo &i which tj{P) has exponential behavior in terms of the holomorphic local coordinate yoo(-P) defined 



by (B.19): 



t,(P)e-'f'^^-^-(^)" - 1 + 0{y^{P)), P ^ oo. 



(B.75) 



B.2.2. Construction of the Baker- Akhiezer functions. To write down Krichever's formula for tj{P), we define 
a homology basis on X as shown in Figure[B.ll[ With this basis selected, we define a holomorphic differential 



OO 




Figure B.ll. A homology basis for the Riemann surface X. Solid curves lie on sheet 1 
and dashed curves lie on sheet 2. 



a;(P), corresponding constants % and /C, Riemann 
point Po — meromorphic differential ^1{P), and 



0-function Q{z;T-L), the Abel mapping A{P) with base 
corresponding constant k by exactly the same sequence 
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of formulae as in case L, namely (B.20)-(B.33|. Lemma B.l also holds in the current context, with the 
implication that 



Jw^ S+{s) m J a 



Also, we see that ipj can be expressed in terms of H as follows: 

TT iH TT iH 

</?i = J^---— and '/'2 = 2 + 
The Krichever formula for the Baker- Akhiezer functions tj (P) is then 

tj{P) -.^Nj r^/./m , r-.-iv^ ^P^i / ) 



(B.76) 



(B.77) 



e(A(P)+/C;H) 

e{A{P)+IC-2iipj;H) 
Q{A{P) + IC;n) 



exp ipipj / ^{P') 



wa 
P 



(B.78) 



where Nj is a constant chosen to enforce the normalization condition (B.75I. 

To compute the normalization constants, choose a path from P = to P = oo on sheet 1 with Sjui} > 0, 
and obtain 



A[oo) = /w(P) = -TTZ. 

2 J a 

Also, as P tends to P = cx) along such a path (we may take P on sheet 1 with w{P) > 0), 



(B.79) 



[\{P') = \i>n{P')-\i>n{P') 



to(P) 



C 



2S{w) 



dw 



(B.80) 



2S{w) 



and 



w + C 
2S{w) 



dw — 



0{P) 



w{P) 



W + C 1 
2S{w) ' 2^" 
w + C 1 



1/2 
1/2 



dw + w{Pf/^ 
dw 



(B.81) 



2S{w) 2 J ^ ' y^{Py 

The remaining integrand is integrable at infinity, and doubling the contour of integration and rotating the 
contours through the upper and lower half-planes respectively to lie along the negative real axis gives 

v[P) 



Jo 



w + C 1 
2S{w) ^ 2' 



-1/2 



dw 



w + C 1 



-w 



-1/2 



and therefore 



It follows that 



n{p') = 



2S{w) 2 
\jn{P') + 0{y^{P)) 

\n + 0{yUP)l 



dw + 0{yoo{P)) 



(B.82) 



yoo{P) 



lim tj{P)e 



-ipiPjVoaiP) 



ei^n - 2iipj;n) 



P^oo 



so the normalization constants are given by 



^' Q{\n~2iipj-ny 
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(B.83) 
(B.84) 

(B.85) 



completing the construction of the Baker- Akhiezer functions tj{P), and hence of 0°"*(w). To obtain a useful 
expression for 0°"*(w) we recall the function defined by the formula (B.51) for w E C \ ((— 00,^0] U 
[wi, 0]); the analogue of Lemma B.2 in the current context is then the following. 

Lemma B.3. l{w) = h{w). 

Proof. The proof is virtually the same as that of Lemma [B.2| one uses the identity (B.76) to show that l(w) 
and h(w) satisfy the same additive jump conditions on the intervals (— oOjWq) and (wi,0), and then uses 
this information to establish that the ratio (/(w) — h{w))/S{'w) is an entire function of w that vanishes as 
w — >■ 00. □ 

To write formulae for 0°"*(u;), it is convenient to introduce the scalar function 

H S{w) 
n{w) := —h{w) — 



S{s){s — w) 



(B.86) 



0° 



Then, with Pj{w) denoting the point on sheet j of X corresponding to w, for G{-u;} > and w G Too: 

J (<2(Pi)e-"^^''-" - ii(Pi)e-'''^i''+") 
for > and w £ Tq: 

i (ti(P2)e''^i''-" - ^2(^2)6*'^"''+") e"" 



I (ii(P2)e*'^i''-" + i2(^'2)e*'^^'^+") 



(B.87) 



0°"*(u;) 



(ti(Pi)e-^'^i''+" + t2{Pi)e-^^'''-") 



_2_ 
2i 



6°"*(w) 



(ii(P2)e^'^i''-" + t2(^2)e''^'''+") e"" | (t2(Pi)e"^'^^''"" - ti(Pi)e-^'^i''+") 

for Sjui} < and w £ Tq: 

" f (i2(Pi)e"*^^'*~" - ti(Pi)e-*^i'*+") e-'' 

^ (ti(Pi)e-*^i'*+" + t2(^'i)e-''^^''-") e-' 
and for Syjui} < and w £ Too- 

I (ti(P2)e^'^i''-" + t2(^2)e''^=''+") -§i {t2{Pi)e-^'^^^-'^ - ti(Pi)e-*'^i^+") 

-§l (ti(P2)e^'^i''-" - t2(P2)e^'^2''+") I (ti(Pi)e-''^i''+" -|-t2(-Pi)e"^'^^''~") 



I (ti(P2)e*^^'*-" +t2(P2)e*'^^'*+") 
I (ii(P2)e*'^i'*-" - i2(^2)e''^'''+") e" 



(Bi 



(B.89) 



6°"'(u;) 



(B.90) 



It is then easy to check from these formulae that the exponentials e^^ViHw) ^^.g exactly cancelled out by 
corresponding exponentials in the Baker- Akhiezer functions; indeed assuming the path of integration to be 
confined to the half-plane containing w for 5{ti;} 7^ 0, 



and 



e(i-H; 'H)Q{A{Pi{w)) + JC - 2i(pj;n) 
" e(i-H - 2iipj;H)QiA{Pi{w)) + JC;H) 

e{^H;H)e{A{P2{w)) + IC - 2iipj;H) 
Qi^H - 2iifj]'H)Q{A{P2{w)) +IC;H)' 



(B.91) 



(B.92) 



due to Lemma B.3 Therefore, all remaining dependence on v is either within the arguments of the 0- 
functions or in the exponential factors e**", leading in both cases to bounded oscillations as — > 00. This 
essentially completes the proof of Proposition |5.2| in case R. 



B.3. Recovery of the potentials. Proof of Proposition 5.3 Here we obtain explicit formulae for C, 
S, and G as defined from 6°"*(w) by ((5^-(|5^. 



B.3.1. Formulae for C , S, and G in case L. To obtain such formulae in case L, it is enough to evaluate the 
asymptotic behavior of 0°"'(i(;) for Syjui} > using (B.52). The value of A{Pi{w)) ioi w k, Q and w ~ oo 
may be determined by choosing the path in ( B.57[ ) to lie in the upper half- plane, avoiding (except of course 
at the endpoints w = and w = or w = oo) the branch cuts of S. Then it is easy to see that 



^(^i(o)) - Ij^^iP) ^ = 



(B.93) 



and from (B.44) we have ^(Pi(oo)) = — H/2. Also, using the fact that w^/^ = i^—w (principal branches) 
for > 0, we have 



A{P,{w)) - A{P,m = 



cd£, 

W) 



Kiei/2(i + o(0) 



2c 

l^^ol 
2ic 

\wq\ 



rW 



+ 0(u,3/2) 



(B.94) 



> 



and 



A(Pi(w))-A(Fi(oo)) = 



cdi 



cd^ 

C3/2(1 + 0(C-1)) 

= -2c«;-i/2 + 0(^-3/2) 
The final ingredients needed are the asymptotic formulae for q{w): 



(B.95) 



oo. 



> 0. 



q{w)^e'^ +0{w), 0<C 



1 TT 

-arg(«;o) < -, 



w 



and 



(B.96) 
(B.97) 



q{w) — l + 0{w 1), ^ oo. 

Using the Taylor expansion of the entire function Q{z;'H) with respect to z and the identities (B.25) and 
(B.26), it then follows from substituting (B.54) and ( |B.56 ) into (B.52) and using ip2 — Vi 



TT that 



and 



A0,0 



A0,0 
^12 



■ 0,0^ e< e(z^;H) 
22 2 9(0; 

.o,o_ e'^e(i7r;H) 

^^21 



(B.98) 



(B.99) 



2i 6(0; -H) 

Note that the matrix 0°'° = O°"*(0) is invariant under conjugation by a2, as must be the case since by the 
conditions of Riemann-Hilbert Problem 



5.1 



are sufficient to calculate C and S: substitution into (5.271 gives 



C=i-1) 



0°"-*{w) is to be Holder continuous up to M+. These formulae 

(B.lOO) 



2 6(0; -H) 



and substitution into (5.28) gives 

^ = ( -|-i#Ae"' Q(^^;"^) 

^ ^ 2i 6(0; 



(B.lOl) 
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The higher-order coefBcients that we need to calculate GN{x,t) are: 



AO,! 



o 



12 



' 0(O;H) 

"K^ 9(0; 



e(^(p2;H) 
e'(i^i;H) 



e(zv?i;H) 



(B.102) 
(B.103) 



and 



Aoo,l 



(B.104) 



Here by Q'{z]'H) we mean the partial derivative with respect to z holding % fixed. Substituting into (5.29) 
we have 

e(i7r;-H)21 \Q'{iip2]H) Q'ivfi^'H)' 



G 



1 



hoi e(0;H)= 



6(1^2; e(z(pi;H) 



(B.105) 



Evaluation ofH. Transformations of Q with respect to %. We now seek to simplify the formulae for C, 5, 
and G. Of course one would like to introduce Jacobi elliptic functions at this stage, but it turns out that 
since T-L is complex the resulting formulae involve elliptic parameters m that are not in the so-called normal 
case of < m < 1. To arrive at the simplest formulae in terms of Jacobi elliptic functions with a parameter 
m G (0, 1) we will therefore first make some transformations of %. Elementary contour definitions in the two 
integrals involved in the definition (B.23) show that T-L can be written in the form 

1 



n=-{'H^ + 2m), 



where "Hq is given by the ratio of integrals 



"Ho 



-27r 



-t-C30 



dw 



^Jw{w - wo)(w - uf^ 
dw 



(B.106) 



(B.107) 



'-00 ^ — wiw — Wo)iw — Wq) 

where the square roots are both positive. Using the substitution 

, .z- 1 



1 



followed by 



and recalling the complete elliptic integral of the first kind defined by (1.60) we obtain 



dw 



2K{cos{Cr) 



^/w{w - Wo){w - Wq) 

Similarly by means of the substitution (4.45) we find 

dw 



2ir(sin(C)2) 



Therefore, 



^y-w{w - Wo){w - W^) 

K{m') 



(B.108) 



(B.109) 



(B.llO) 



Ho 



-2tt- 



(B.lll) 

1 — m €E (0, 1) are given by 
(B.112) 



K{m) 

where the elliptic parameter m € (0,1) and its complementary parameter m' :- 

m := sin(<;^)^, m' := cos(C)^. 

Note that Hq is a negative real number. According to the theory of elliptic functions (see P^) it is Riemann 0- 
functions with parameter Hq that are associated with Jacobi elliptic functions of elliptic parameter m € (0, 1). 

Adding 27ri to Hq and then dividing by two amounts to the composition of two classical transformations 
of 0-functions [1 . Indeed, the so-called first principal first- degree transformation implies that 

e{z;H-^2Tri) ^Q{z + i7T;H), z e C, ^{H} < 0, (B.113) 
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an identity that allows us to add 27ri to Hq. Also, the so-called Gauss transformation (a second-degree 
transformation) implies that 



e(i7r; ^n)e{z; \n) = e(z + m- nf -f e-^e^/^Q{z + iTT- 1%; Uf, z e 



and 



6(0; \n)Q{z + in; \n) = Q{z + ivr; Uf - e-'e^'^Q{z + in - 1%; nf, z e 



n{H} < 0, (B.114) 
^{n} < 0, (B.115) 



identities that allow us to divide Hq + 2ni by two. Combining these together and using the identity (B.26) 
yields 

e(^^;H)e(z;H) ^ e(z; Hp)^ + ^e--e^o/^e(z + »^ - ^Hp; Hp)^ -u^li-u +o ^ milfi^ 
&{z + in;H)Q{Q;n) e(z; Hp)^ - ie-e«o/4e(^ + _ i^^^; T^o)^ ' 2^"°^ ^''^ ^ ' ^ 

The elliptic parameter m can be expressed directly in terms of special values of 0-functions with parameter 
"Hp as follows [T]: 



9(^7r;Hp)^ 
e(0;Ho)4 



1 — m and 



p«o/2 



6(— ^Hp; "Hp)** 
e(0;Ho)4 



(B.117) 



Since Hp < we see from (B.24) that all four 0-function values appearing in this identity are real, so we 
may take the positive square root of both sides. 



We may now introduce the Jacobi elliptic functions with parameter m related to "Hp by (B.lll) or 
equivalently by (B.117): 



cn 



K{m)z 



K{m)z 



-zjl 



and 



Setting 



dn 



K{m)z 



9(0;Hp)e(z + z7r~- |-Hp;-Ho) 
e(-iHo;Ho)e(z + ^7r;Ho) ' 

9(»7r;7^o)9(z- iHo;Ho) 
e(-iHp;Ho)e(z + i^;Ho)' 

e(i^;Ho)e(z;Ho) 



z e 



z e 



e(0;-Hp)e(z + i7r;-Hp)^ 



K(m) 

Zj := 



(B.118) 
(B.119) 

(B.120) 

(B.121) 



and substituting into ( |B.100[ ) and ( |B.101[ ) firstly from ( |B.116| -( |bH71 ) and then from ( |B.118[ )-( |rI201 ) we 
obtain: 



^=(-l)#^- 



2 
27 



dn(Zi; rn)^ — i\J ram' sn{Z\ \ mf dn(Z2; mf — i\/mm'sn[Z2', mf 



dn(Zi; to)2 -|- i\J ram! sn{Z\ \ mf dn{Z2', mf + iy/mm'sn[Z2', mf 

dn(Z2; mf — i^mm'sn(Z2', raf dn(Zi; mf — i^mm'sn{Zi; mf 
dn(Z2; to)^ -I- i^/mm'sn(Z2', ra)^ dn(Zi; to)^ -|- i\/mm'sn(Zi; m)^ 



(B.122) 



(B.123) 



To express G in terms of Jacobi elliptic functions, first we note that using z — in in (B.116I, taking into 
account the periodicity relation (B.26), and then using the positive square roots of the identities (B.117) 

■e{in;'Hf 



together with (B.112) gives 



1. 



(B.124) 



Now, by simple contour deformations and the use of the normalization condition (B.22) defining c, it follows 
that 2cl? = TT, where V is the denominator defined by (4.33). Therefore, if we allow the branch points wq 
and Wq to depend on {x,t) via the moment and integral conditions M = / = 0, it follows from (4.32) of 
Proposition 



4.2 



that we may write G in the form 



G = 2 



~dt 



e{iip2-H) e{i^i;H) 
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(B.125) 



Since <P2 = <;5i + tt, this can be written as a logarithmic derivative: 

'ld{i(pi + in; Ti) 



G 



^.9$ d 

dt dipi ^ 



(B.126) 



Applying (B.116) and the positive square roots of (B.117), we may then substitute from (B.118)-(B.120) to 
obtain 



^ 2 '^^ Y ( dn (Zi; m)^ + iVmm'sn (Zi; m)^ 
dt dipi \ dn (Zi; m)^ — i\/ mm'sn {Zi; m)^ 



(B.127) 



Use of elliptic function identities for a fixed elliptic parameter. Now, we simplify C, S", and G further by 
recalling some identities relating elliptic functions at a fixed value of the elliptic parameter m. Firstly, using 
the Pythagorean identities [T] 



sn(-; m)^ + cn(-; m)^ — dn(-; m)^ + msn(-; m)^ = 1 
to eliminate dn(Zj;r7i)^ and recalling that e"> = \/?7i' + i^/m, we have 



and 



G 



2 

2^ 



1 — iy'roe *^sn(Zi;m)^ 1 — zy^e '^sn(Z2;m)^ 



1 + iy'me*^ sn(Zi; m)^ 1 + i\/rrie'?sn(Z2; m)^ 
1 — iY^e~*''Sn(Z2; m)^ 1 — iy/me~'^'^sn(Zi; m) 



1 + iY^e'^sn(Z2; m)^ 1 + iyme*^sn(Zi; to)^ 



2i——\o 
dt dipi ^ 



1 + iy/me'^'^sn{Zi; m)^ 



1 — i-^/roe '^sn(Zi; m)^ 



Next, using the double-angle identity 

sn(Z; m) 



these can be written in the form 



2 77icn(2Z; m) — dn(2Z; m) + to' 
mcn(2Z; to) — TOdn(2Z; to) ' 

•y/TOe^""cn(2Zi; to) + idn(2Zi; to 



iVto'c 



■\/ TOe*? cn(2Zi ; to) — idn(2Zi ; to) + i\fm! e'? 
•\/TOe^*^cn(2Z2; to) + idn{2Z2', to) — i\fm'e~^'^ 



S={-\) 



2i 



TOe*^cn(2Z2; to) — idn(2Z2; to) + iy/rn'e^'' 
TOe^*^cn(2Z2; to) +idn(2Z2;m) — i^fm'e^^'^ 



TOe*^cn(2Z2; to) — idn{2Z2; to) + iy/m/e^i 
TOe^'^cn(2Zi; to) + idn(2Zi; to) — i\/rn/ e~^^ 



and 



Then, since 



G 



TOe'^cn(2Zi; to) — idn(2Zi; to) + m'e*^ 

2.9$ d ^ / -y/TOe"'Cn(2Zi; to) — idn(2Zi; m) + iVTO'e*^ 
* (9t diyJi i v/TOe-^^cn(2Zi; to) + idn(2Zi; m) - i\fm! e-^< 



2Zi^W - K{m) and 2Z2 + if (m) where 
the use of the identities 



2vK{m) 



cn(W ± K(m);m) — ^Vm' — - — ' — - and dn{W ± K{m);m) ^ — -. r 

dn(W^; to) dn(T/F; m) 



(B.128) 

(B.129) 
(B.130) 

(B.131) 
(B.132) 



(B.133) 



(B.134) 



(B.135) 



(B.136) 



(B.137) 
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yields 



me *^sn(W; m) + i — ie '^dn(PF; m) 



me*'^sn(Vl^; to) — i + ie*?dn(M^; m) 

TOe~*^sn(VF; to) + « — ie"*^dn(T4/'; to) 



2^ 



^me^'^sn{W; m) — i + ie'-'^dn{W; to) 

me~*^sn(VF; to) + i — ie~*^dn(VF; m) 
— ■\/me*fsn(Vl^; to) — i + ie*^dn(W^; m) 
•\/TOe~*''Sn(M^; m) + i — ie^"^^ dn{W \ to) 



(B.138) 



and 



^ AK(m)d^ d , 

G = — : — — log 



TOe*^sn(Vr; to) — i + ie"^dn(VF; to) 
msniW; m) — ze~*^ + idniW; to) 



JTT dt dW °\ Y^sn(VF; to) + ze*^ — idn(VF; to) 
Again applying (B.128I and e'^ = \/m' + i^/ro, and the differential identities jl] 

d 



(B.139) 



(B.140) 



dW 
d 

dW 
d 

dW 



sn(T/F; m) = cn(VF; m)dn(W^; to) 
cn(T/F; m) = — sn(T4/'; TO)dn(M^; to) 
dn(T/F; to) = — TOCn(Vl^; TO)sn(W^; to), 



(B.141) 



these become simply 



C = (-l)#'^dn(VF;m) 

5* = -(-l)#'^VTOsn(W^;TO) 

G = ^ — ^ — VTOcnfW^iTO). 

TT Ot 



(B.142) 



We have already introduced (a:, t)-dependence via the conditions M = I = 0. If we also recall that 
v = $/eAr + 7r#A, then W takes the form 



(B.143) 



But, since is even in case L, and since sn(-; to), cn(-; m), and dn(-; m) are all periodic with period 4ii'(TO), 
the formulae (B.142) reduce to the expressions (5.30), as desired. 



This nearly completes the proof of Proposition 5.3 in case L. It only remains to confirm the differential 

(B.144) 



= u, 



relations (5.34). Note, however, that by partial differentiation with respect to to of the relation 

.sn(«;m) 

Jo VT^Vl - mt^ 

the use of the Pythagorean identities (B.128) shows that 

d 



-sn(u;m) = —cn(u;m)dn(u;m) 
dm jQ 



a(ti:m) g 

dm 



1 



VF^Vl - mf^ 



dt. 



(B.145) 



Now the integral increases by 4:K'{m) when u increases by AK{m), the fundamental real period of sn(-; to) 
Therefore, 



_d_ 
dm 



sn(u; m) ~ — cn(u; m)dn(u; to) 



K'{m) 
K{m) 



u + /(m; to) 



(B.146) 
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where f{u + AK{m)]m) — f(u;m), making / periodic and hence bounded with respect to u. By partial 
differentiation of 5jv(a^,i) as given by the formula (5.30) we obtain 

dSN 



eN dm ^dm d 

eN^—(x,t) = -— -=— -sn(u;m) - eNVm—--—sn{u;m) 
at 2y/m ot at am 



TO— — sn( w; to) 
ou 



2K(m) a$ 2<i>K'{m) dm 



dt 



dt 



2$X(to) 



(B.147) 



Using (B.141 ) and (B.146I to evaluate the partial derivatives of sn(u; to) (note that the derivatives in (B.141 ) 
are in fact partial derivatives with to held fixed) this becomes 



dSN , . 



2K{m) 5$ 
^ ~dt 



TOCn(u; TO)dn(u; to) 



dm 
dt 



TOCn(u; m)dn(u; m)f(u; to) ^sn(u; to) 

2VTO 



(B.148) 



and the first of the relations (5.34 1 then follows upon comparing with the formulae (5.30). The second 



relation is established in a completely analogous manner. This finally completes the proof of Proposition 5.3 
in case L. 

B.3.2. Formulae for C , S, and G in case R. Once again to obtain the asymptotic behavior of 0°"*(w) as 
w ^ Q and w — > oo, it is sufficient to assume that Sjui} > 0, and therefore to use the formula (B.87) to 
analyze the limit ui — > oo and the formula (B.88) to analyze the limit w — > oo. By choosing a path in the 

(B.149) 



open upper half-plane, we easily obtain the following asymptotic formulae for the Abel mapping: 



A{Pi[w)) = -H -ii: + V-w + 0(w^/^), w^O, > 0, 

2 ^wqWi 



and 



A{Pi{w)) ^ -iTT + 



2ic 



0{w 



-3/2N 



w oo, Syjui} > 0, 



and of course A{P2{w)) — -~A{Pi{w)). Also, we have the following asymptotic formulae for q{w): 



q{w) 



1/4 



+ 0{w), w^O 



and 



q{w) ^ I + 0{w^^), w ^ oo. 



(B.150) 

(B.151) 
(B.152) 



We will also require asymptotic expansions of n{w) defined by (B.86) valid for large and small w. Beginning 
with the exact formula 



n{w) 



Siw) 



ds 



n 



27ri J^, S+{s){s-w) 



ds 



S{s){s — w) 



(B.153) 



we may use the expansion S{w) = w ^^^{l + 0{w ^)) as w — > oo and expand {s — w) ^ in a geometric series 
for large w to obtain 



,1/2 



nlw) 



n 



n 



s ds 



2T:i J^^ S+{s) 
But, by definition of we may write this in the form 



ds 




ds 


S+{s) ^ 


J Wo 


S{s) 




s ds 




) Jw„ 


S{s) 



(B.154) 



n{w) 



1 



2^1/2 



n 



c 



2ni J^^ 2S+{s) 



ds 



C 



2S{s) 



0{w- 



ds + 0{w-^) 



w — )■ oo. 



w — )■ oo. 



(B.155) 
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Identifying the integrals as cycles of the meromorphic differential ^1{P) then yields 



n{w) 



1 

2^1/2 
c 

2i\/—w 



(B.156) 



where in the second line we have used Lemma B.l On the other hand, we may write n(w) in the form 



n{w) 



Siw) 



n 



ds 



ds 



(B.157) 



Ani S{s){s~w) J^^ S{s)is-w)_ 

where it is understood that w lies outside of the loop contour a pictured in Figure [B.ll[ If we let w cross 
the contour a to approach the origin, then we obtain a residue contribution: 



, , H S{w) 

n[w) = 

^ ' 8 4 



n 



ds 



ds 



(B.158) 



4ni S{s){s-w) J^^ S{s){s~w) 

where it is now understood that w lies inside of a. In particular, the expression in brackets has a well-defined 
value at u; = (it is in fact analytic in a neighborhood of u; = 0). Since S{'w) — ^/wowiw^^^^l + 0{'w)) as 
w — >■ 0, we therefore see that 



n{w) = — 
The differential identity 



n 



ds 



d{w) 



2TTi S{s)s 
(w — wo){w — Wi) 



"'1 ds 

,0 '5'(S)S 



3/2^ 



d\w) = 



WqWi 



S{w) ' ' 2S{w) 2wS{w) 

allows integration by parts, leading to the equivalent expansion 



(B.159) 



(B.160) 



1 



8 8^woWi 



n 



2^1 J, S{s) 



s ds 

W) 



3/2N 



(B.161) 



Using the definitions of % and the meromorphic differential ^1{P), this may be written in the form 



n{w) — — 



I 



8 4^woWi 

n K 

V 

8 A^WqWi 

n c 

8 2i^WQWi 



n 

'2Tii 



Ja Jb 



(B.162) 



where Lemma B.I has again been used. 

Applying these expansions for w small with Syjiu} > to the formula (B.88I with the help of (B.9I |-(B.92) 
then yields 



" 22 2 \wi 



and 



AO.O _ A0,0 
"-^12 — *^21 



6(0; -H) 
9(0; -H) 



e(2z(^i;H)^,, ^ ei2iif2;n) ^_ 



e{2iip2;n) 



e{2iipi;n) 



2i \wi 



/8 



e{2iip2;n) 



Q{2iipi-U) 



(B.I63) 



(B.I64) 



As in case L, we observe t hat t he ma trix °'° = O°^'(0) is invariant under conjugation by (T2- Substitution 
from these formulae into (5.271 and (5.28) gives, respectively. 



C 



(-i)#^ (w,\'/^Q{\n-.n)^^/ 



Wi 



9(0; 



and 



S 



(-i)#^ ( w^\"^ Qiln^n) 



2i 



Wi 



6(0; 



9(2»yi;H) 6(2z^2;H) _ 
Q{2iip2;nf Q{2iipi-nf 



6(2z^i;H)^,, Q{2i^2;'H) ^_ 



6(2z(^2;H) 



e(2i^i;H) 



(B.I65) 
(B.I66) 
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The higher-order coefficients we need to calculate G are obtained by continuing the expansion for w near 
the origin to higher order, and also by expanding (B.87I as w — > oo with Qjw} > with the help of 
(|R9Tl)-(|R92l): 



qO.I 



9(|'H;'H) ^/g 



6(0; -H) 



ie(2i^i;H) i e(2»^2;H) _ 

-e — - „ . —re 



4e(2^^2;H) 



4e(2iv?i;-H) 



(B.167) 



'^12 — 



Wo 



1/4 o/l 



e(i-H;-H) 

eio-n) ' 



e{2iip2;H) 



e(2z^i;H) 



ie(2zvPi;H) „ , 1 9(2*^2; 



and 



Aoo,l 
"-^12 



4e(2i(/j2;'H) 4e(2i^i;-H) 

e'{2iip2;n) _ Q'{2iipi;H) _ 1" 
e{2i(p2;H) Q{2iipi;H) 2 



Substituting into (5.29) we have 



G 



wi 9(0; H)2 



Q' {2iip2;'H) _ Q'{2iipi;H) _ 1 
e{2i(p2;H) e{2iipi;H) 2 



(B.168) 



(B.169) 



(B.170) 



Evalua tion o f % . Transformations of 9 with respect to % . By simple contour deformations (referring to 
Figure B.ll I and the definitions of S{P) and of S{w), we see that 

'""1 dw 

y^w{w — Wo){'W — Wl) 



n 



^27r 



dw 



(B.171) 



1 ^y-w{w - Wo)iw - Wl) 



where in each case the positive square root is meant. Although this is a negative real quantity and hence is 
associated with an elliptic parameter that is in the normal case, it will be in fact convenient to rewrite H in 
the form H = 2'Ho and associate Ho with an elliptic parameter m. To this end, we recall the substitution 
(|4.47|) and obtain 



dw 



2K{m) 



where 



1 ^—w{w — wo){w — Wl) y/^wo + y/--wi 
4,/woWi 

V <' i ^ (Q^-^)^ 



iV-Wo + y/-Wi)^ 

For the denominator of we use the globally bijective fractional-linear substitution 

- {y/-Wo + y/-Wi) - {y/-Wo - y/-Wi)s 



W 



-^^WqWi 



to obtain 



iV-Wo + y/-Wi) + (V-WO - y/-Wi)s 

dw AK{m') 



^Jw{w — Wo){w — Wl) \/~Wo + \/—Wi 



where 



:= 1 — m. Therefore, it follows that 



T~L — 2'Hoj 'Ho 



-2tt 



K{m') 
K{m) 



(B.172) 
(B.173) 

(B.174) 
(B.175) 

(B.176) 



where the elliptic modulus m is given by (B.173). 
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Multiplying by two corresponds to a Gauss transformation (in reverse). If we try directly to apply 
(B.114|-(B.115) to express 0-functions with parameter H in terms of those with parameter Ho we will need 
to choose branches of square roots. This difficulty can, however, be circumvented by using the formula [26 



e{2z;2n) 



29(0; 2-H) 



2 e 



^{H} < 0. 



(B.177) 



We may now rewrite (B.165) and (B.166) in terms of 9-functions with parameter Hq with the help of 
(B.177 1, and then substitute into these formulae the definitions of the Jacobi elliptic functions (B.118)- 
(B.120I. Using also the positive square roots of the identities (B.117I, the standard 0-function properties 
(B.25 |-(B.27) along with the fact that 8(i7r + ^^Hq'jHq) = 0, and also the identity 

1/4 



Wo 
Wi 



1 



following from (|B.173|), we see that 

dn(Zi; m)^ 



C 



and 



2i 



dn(Z2; m) 



Vmm'sn(Zi; to)^ — 
dn(Zi;m)2 



dn(Z2; m)^ + Vm' 



\/mm,'sn(Zi; m)^ — ^ymcn{Zl^, m)^ •\/mm'sn(Z2; m)-^ — ^ymcn{Z2', m)^ 



(B.178) 



(B.179) 



(B.180) 



where Zj are defined in terms of ipj by (B.121I (although in the current case ipj and m have different 
meanings than they did in case L). 

To express G given by (B.170) in terms of elliptic functions, we first note that since H = 27io, (B.177) 
together with the positive square roots of the identities (B.117), the definition (B.173I, and the fact that 
6(i'Ho + itt; "Ho) = imply that 



1 9(^H;H)^e«/4 
e(0;-H)2 



^/wowi 



(B.181) 



Elementary contour deformations show that cD = tt, where V is defined for case R by (4.34 1. If we now let 
Wo and wi depend on {x,t) such that the moment and integral conditions M — I — are satisfied, then 
Proposition (4.2 1 applies, and from (4.321 one finds that 

1 



1 e(iH;H)2e«/4 



dt dt ■ 



Also, since ip2 = 



wi 9(0; 7^)2 

iH/A, we may write G in the form 

,^^9(2»(pi^iH;H) 



G = 



2i——\o 
dt dipi ^ 



9(2i^i;H) 



(B.182) 



(B.183) 



Substituting from (B.118)-(B.120) after first using (B.177I and the positive square roots of (B.117I we arrive 
at the formula 

cn(Zi;m)2 - \/m'sn(Zi; m)^ \ 



2i——\o 



dn{Zi ; m)2 + \/ m' 



(B.184) 



Use of elliptic function identities for a fixed elliptic parameter. Using the Pythagorean identities (B.128) 
allows us to write G and S in terms of sn(Zi; to)^ and sn(Z2; m)^ only, and G as a logarithmic derivative of 
a quantity involving sn(Zi;m)2 only. Next, we apply the double-angle identity (B.132) and note that 

K{m) + iK{m) and 2Z2 ^ W + K{m) - iK{m') where W : 



2Zi = W 
allowing the use of the identities 

cn(M^ ± K{m) =F iK{m'); m) = 



2vK{m) 



(B.185) 



i\l m' 



and dnlW ± K{m) ^ iK{'m');m) — ^iVm 



jSn{W; to) 

cn{W] to) 
(B.186) 
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Finally, we again apply (B.128), and in the case of G the differential identities (B.141), to express C, S', and 
G simply as: 

C ={-!)* ^cxv{W]m) 
AK{m) 5$ 



G 



-dn(W^; m) 



TT dt 

Finally, inserting the value ofz/asz^ = $/ejv + tt^A, the phase W becomes 

2^K{m) 



W = 



(B. 



and since sii{u + 2K{m)]m) = — sn(u;m), cn{u + 2K(m);m) = — cn(u;m), while dn{u + 2K{m);m) — 
dn{u;m), we see that even though is not necessarily even in case R, the formulae (B.187I indeed have 
the desired form (5.32). 

To prove that (5.34) holds also in case R, we differentiate SNix,t) given by (5.32), keeping in mind that 
m — m(x^ t): 



dm d , , 



d_ 
du 



sn(u; m) 



2'^K'{m) dm 2K{m) 9$ 



dt 



dt 



Computing the partial derivatives of sn(u;TO) using ( |B.141 ) and (B.146) this becomes 



dSN I . 



2K{m)d<^ dm 

— cnfit; TOjdnlM, m) + ejv-;:^cn(u; mjdnlu; mjjlu; m) 

TT dt dt 



2^K{m) 

Trejv 

(B.189) 



(B.190) 



from which the first of the formulae (5.34) follows by comparison with (5.32) (the second is proved analo- 
gously). This completes the proof of Proposition 5.3 in case R. 
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